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ducting astronomical research. Possessing MSLAP in whole or part indicates that the user has 
accepted all of the terms of the licensing agreement, which serves to insure the integrity of the 
software. Basically, the terms of the license agreement are: 

1) You are not permitted to distribute the standard source code to any other site. There are several 

distribution centers that supply standard source code. 

2) You are permitted to distribute additional modules (subroutines) written for use with MSLAP to 

other sites, provided the new code meets documentation requirements. 

3) You are permitted to alter the standard source code as you desire, but you must maintain the ori- 

ginal copyright notice in the source file and you may not distribute the altered code without 
written authorization. Also, any alterations to the standard source code must meet documenta- 
tion requirements. 

4) If you represent a guest user facility, you are not permitted to incorporate modifications that 

remove the modularity of MSLAP, thus making it difficult for guest users to customize 
MSLAP. 
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I. Introduction to MSLAP 


The Modular Spectral Line Analysis Program (MSLAP), as its name implies, forms a back- 
bone of programs, designed so that customized subroutines can be inserted and implemented 
with minimal difficulty. MSLAP, a third generation package of software, also is a complete and 
powerful stand-alone program for analyzing spectra, providing the basic structure to identify 
spectral features, to make quantitative measurements of these features, and to store the mea- 
surements for convenient access. MSLAP can be used to measure not only the zeroth moment 
(equivalent width) of a profile, but also the first and second moments. Optical depths and 
the corresponding column densities across the profile can be measured as well for sufficiently 
high resolution data. 

The software was developed for an interactive, graphical analysis where the computer carries 
most of the computational and data organizational burden and the investigator is responsible 
for all judgment decisions. Cursors are used not only to provide graphical input, but also for 
logical control branching. It employs sophisticated statistical techniques for determining the 
best polynomial fit to the continuum and for calculating the uncertainties. MSLAP, making 
use of data structures, provides substantially more capabilities in the handling, presentation, 
editing, and manipulation of intermediate results than do its predecessors. 

MSLAP is written in the Interactive Data Language (IDL) and issues some commands to the 
UNIX operating system. While MSLAP has not been transported to a VMS environment, 
there is no apparent reason why it could not be. MSLAP is specifically designed for worksta- 
tions, running in either an X Window or Sunview environment. 

The basic structure of this document is as follows. A method for getting started quickly is 
given in §11, where it is assumed that MSLAP has already been installed in some convenient 
library directory and the reader is anxious to try running MSLAP for the first time. Section 
III provides a detailed, step-by-step description of the six major portions of the program. 
Some common mistakes in data handling are listed as cautions in §IV. The documentation for 
modifying and customizing MSLAP is given in §V and for installing MSLAP on system is in 
§VI. Finally, for convenient access, Appendix A contains a listing of the source code, including 
a table of contents, while the three data structures used by MSLAP provided in Appendix B. 


II. Quick Start 

It has been said that a fundamental difference between men and women is that men do not 
like to read directions, preferring to get started quickly, while women do read the directions, 
opting to do the job only once. Whatever your persuasion, this section is designed to get the 
novice going with the minimal amount of reading. Once program has been started, the novice 
can rely on the instructions provided by MSLAP as it runs or rely on §111, which provides a 
more detailed description. 


a. Getting Ready 

MSLAP runs in either an X Window or Sunview window environment. Sunview is perhaps 
the more sophisticated system, but is specific only to Sun Workstations, while X Windows is 
a industry standard that is available on most machines. X Windows has the added advantage 
that a user can run software on a remote machine, displaying the graphics and printed matter 
on the local console. It is not possible, however, to run MSLAP on a remote machine without 
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a local workstation or terminal that supports windows. 


Both Sunview and X Windows can be configured in many different ways. If you are at a Guest 
Facility, the chances are that your machine automatically starts one of these systems when you 
log onto the workstation. MSLAP needs a large window in the lower left corner for printing 
instructions. If you do not have a window in the lower left corner and are not an experienced 
user of one of these environments, it is recommended that you copy the appropriate files from 
the MSLAP library to your login directory. Your system manager has placed the MSLAP 

library in the directory: For X Windows, the files are .xlldefaults, 

.twmrc, and .xinitrc. For Sunview, only .sunview is needed. If necessary, make sure that your 

path includes: (Your path is found either in your .login or 

.cshrc file and can be changed using your favorite editor.) Then type: "source .login” to reset 
your path name. Start the window system by typing either: 


"sunview” 


or 


"xinit ; kbd_mode -a ; clear” 

for Sunview or X windows, respectively. The latter could be aliased to something simpler or 
either could be place in your .login file. Now each window that appears, effectively represents 
a separate terminal. Each time the cursor, associated with the mouse, is placed into a window, 
that window becomes the "active terminal”. (In some cases, it is also necessary to click one 
of the mouse buttons to activate the new window.) All other windows, however, can still 
continue processes that were previously started, including output printing. 

Place the mouse cursor into the large window in the lower left of the screen. Change to the di- 
rectory where you want to run MSLAP and copy from the MSLAP library the files: mslap.pro 
and dgets.pro You are now ready to proceed to run MSLAP. 


b. Running MSLAP 

MSLAP is written in the Interactive Data Language (IDL). Use the mouse to place the cursor 
in the large window in the lower left corner of the screen and start IDL by typing: "idl”. Once 
you see the prompt ”IDL>”, type: ".ran mslap.pro". A menu of options will then appear. 
The first two options allow the researcher to scan quickly through the spectra, searching for 
features, identifying the species in look up tables, making measurements of the profiles, and 
storing all of this information in a data base for latter use. The second option also produces 
optical depths as a function of velocity across the profiles. 

In options 1: and 2: of the FIRST MENU, a series of questions appear to configure MSLAP. 
The default values for all Yes/No questions is assumed to be NO except for answers of either 
a ”Y” or a ”y”. 

The main body of MSLAP begins after these questions are answered. In this portion of the 
software, the mouse is used to input graphical information to the program as well as to con- 
trol the data reduction. The basic philosophy of the mouse operation is as follows. The left 
button should be used in response to the general flow of the program while the middle button 
is used to signal the computer that the researcher wishes to move on to the next task. For 
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example when the main spectrum is showing, the left button causes the computer to attempt 
to identify features of interest (the first task MSLAP expects to perform) and the middle 
button causes the computer to go fetch the next portion of the spectrum. Another example 
would be in the continuum- fitting portion of MSLAP. There, the left mouse is used to select 
the portions of the spectrum to be fit, while the middle mouse button indicates all portions 
are selected - now go on and do the fits. 

Finally, all other branching is accomplished by using the right mouse button, which causes a 
Menu of additional options to appear. The researcher is invited to examine this Menu from 
time to time. There is no harm in striking the right mouse button since "Take No Action” 
is always one of the options. Note: once the Menu appears, the left mouse is used to make 
the selection. This Menu can be modified easily and customized routines inserted into it by 
editing the mslap.pro file that was transferred to your working directory. (See §V for details 
on modifying this file.) The researcher is encouraged to examine this file once he or she has 
become familiar with the operation of MSLAP. 

At this point put down this manual and concentrate on running MSLAP. The program pro- 
vides self-contained directions that appear in the lower left window. However, if you run into 
trouble, §111 of this manual has a detailed, step-by-step set of instructions. 
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III. Step-by-step Description 


This section is organized into 6 subheading that represent each of the 6 options from the first 
menu that MSLAP displays when started. This menu appear as follows: 


Hay 17, 1990 

9 

9 

J *********************************************************************** 

; Modular Spectral Line Analysis Program 

; M.S.L.A.P. version 1.0 

9 

; copyright (c) 1991 by Charles L. Joseph and Edvard B. Jenkins 
; All rights reserved. A license may be obtained from the 

; first author or from an authorized distribution center. 

» 

; This software distributed through: Princeton University 

• *********************************************************************** 


THE FOLLOWING OPTIONS ARE AVAILABLE: 


1: Measure Profile Moments Only 

(Equivalent Width, Profile Centroid, etc.) 
2: Measure Moments and Optical Depths 

3: Manipulate Tabled Data 

4: Edit Tabled Data 

5: Analyze Column Densities vs Velocity 

on Option #2 Data 

6: Compare Data to Curve of Growth 

10: TO EXIT FROM THIS PROGRAM 


Which option would you like ? 


1: Measure Profile Moments Only 

Option 1: of the first menu is for measuring the zeroth (equivalent width), first (velocity), 
and second moments of the spectral profiles. Both this and and the second option ask the 
same set of questions, which will be discussed below. The computer I/O will be displayed in 
small style type like this, while the descriptive narration describing these questions will be display 
in ordinary style type. Answers to the questions will be provided as an example. 

The first two questions deal with establishing an output data file. In the case below, the 
output file name is called "test” and a file called "test.DTL” will be created or opened, 
depending whether the file already exists. If a file named test.DTL is found, the program 
examines the contents to see how many measurement have been made previously and reports 
this information to the user. The investigator then has the option of starting with a number 
of measurements already made or initializing (erasing) the contents of the file. Note: if you 
wish to save the previous work, but not combine new data to the file, answer the "Append 
new data to the old?” question with a ”y” and exit the program later. For the purposes of 
this example, a ”n" was used. 


Vhat is ths output file ntmt ?t«st 

Fils Already Exists 

4 ■sasursJisnts hays already been aads. 

Append nev data to the old? 

(Vote: a KO will erase old data.)n 


Next the program seeks to determine the radial velocity expected for the profiles. This in- 
formation is used only to help select the correct entries in the Look Up Tables in order to 
identify the various spectral features. The information can be input in the form of Delta- 
Lambda/Lambda or as a Radial Velocity in km/s. If you are uncertain of the real value, enter 
a ”0”. The computer then mirrors the value that has been entered. 


Tb* radial velocity is used only to help identify speciea. 
Enter 0 if you are uncertain of the real Talue. 

Enter Delta-Laabda/Laabda or Radial Velocity for the source: 0 
The Delta-Laabda/Laabda is: 0.00000 


There are several Look Up Tables which contain the identification of various species along 
with important information such as the oscillator strength and rest wavelength. (See §V for 
details of the use and implementation of these tables.) MSLAP has the capability to search 
more than one of these tables automatically. In addition, a personal User Look Up Table, 
which is capable of holding up to 100 entries, is scanned. 

Below is a list of Look Up Table Options. In the example, the investigator has requested both 
the Morton and Smith table of interstellar lines as well as a compilation of molecular bands 
be searched each time an identification is requested. 
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Awailabla lookup tablaa ara: 


1) Intarstallar Linas (Horton and Saith 1973 plus updatas) 

2) Holacular Hydrogan, HD, and CO 

3) Options 13 2 Coabinad 

4) Hot-Star Linas (not iaplaaantad) 


Which ona would you lika ?3 


Next the programs queries as to which dget (i.e. data-getting) routine is to be used and the 
name of the file holding the spectra. There are a number of data formats used by different 
observatories so it is appropriate for many users to keep several dget routines on hand. While 
only 5 dget routines are supported at a time, each researcher gets to incorporate the 5 that 
best meets his needs. (See §V for more details on the implementation of the dgets.pro file.) In 
the example, the investigator has requested a very specialized dget routine which reads data 
taken with an echelle spectrograph on a sounding rocket. This particular dget differs from 
the dget loaded into slot 3) because it contains information about the background level and 
associated errors. The file to be used is in another subdirectory: ../rocket and is named: imaps. 


DATAGET Options: 

1) IUE Standard GO Filaa for High-Ras. Spactra (DISXGET) 

2) HOT Baing Uaad 

3) ASCII Fonat of VsTalangth-Spactra-Quality-BG 

4) 1024-Elanant Stand-Alona Data 
6) XHAPS psaudo standard (ASCII) 


Which ona would you lika ? 5 


Entar conplats IIFUT Data Filanaaa 

Including tha path if nacassary . ,/rockat/iaaps 


The final question before the researcher is off and running is the issue of coherence length. 
The message is pretty much self explanatory, but its importance should not be slighted (see 
§IV). Note that the coherence length need not be an integer value. For example, if the value 
of a pixel is influenced only slightly by its adjacent neighbors, the coherence length could be 
1.2 for instance. A non-integer value might also be expected if the data have been smoothed 
by a running-gaussian smoothing routine. 


DfPORTAIT: Tha coharanca langth is us ad to calculata ALL uncartaintias . 
It is tha nuabar of pixals influancing tha yalua in a giwan pixal. 

Thus, tha coharanca langth wust bs 1 or grastsr. For axanpla, data 
swoothad by a 3-point Running Box Car has a coharanca langth of 3. 

Vhat is tha coharanca langth of tha prasant data? 3 
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The main program then begins and the screen resembles Figure 1. Generally speaking, in- 
structions appear in the text (lower left) window, indicating how the graphical input is to be 
made in the plot window. As this portion of the software starts, the following information 
is displayed in the text window. Particularly important pieces of information are always boxed. 


>•* Data Hara Baas Raad — Data Catting Option: S 


| Tha background uncertainty frequently has a major impact on the I 
| uncertainty of the various measurements . i 
1 Current background is: 12.0000 with an error: 5.00000 I 
| Use the RIGHT Mouse Button if these are unsat is fact 0 x 7 . I 
| lots: all graphical input is performed by placing the cross-hairs I 
| at the point of interest and pressing a mouse button. I 


Left Mouse to Locate Feature to be Measured Heed B0T be centered 

Middle Mouse to GO 01, get new data 

Right Mouse to bring up UeerMenu of other options. 


Notice that the last 3 lines contain the instructions for the use of the Mouse. These instructions 
are constantly being updated as you proceed through the program. Once you become familiar 
with MSLAP, your focus should move almost exclusively to the plot in the upper right window. 



Fig. 1 showing the typical screen layout. 


The plot in figure 1 has several features worth highlighting. In addition to the spectrum, there 
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is a data quality vector plotted. Not all data sets have such a vector, in which case it is set to 
zero as shown in figure 1. The three horizontal lines plotted between about 1125.4 and 1125.5 
A indicate the assumed background level and its uncertainty. This information is also printed 
in the text window, surrounded by a box, every time a new data is read. The graphical cross 
hairs or cursors are also depicted and as the Mouse is moved across the pad, the cursor moves 
accordingly. A small window, labeled Cursor Coordinates, provides continuous readout of the 
cursor location. 


At this point, the branching possibilities become large and it is nearly impossible to provide 
a detailed description of each path on a step-by-step basis. The general philosophy, however, 
is to scam through the spectra, looking for spectral features of interest. The Center Mouse 
Button is reserved to indicate that the user is finished with the present task. Since the present 
task is to locate features to be measured, the Center Button indicates a request for new spec- 
tra. If your dget is configured for multiple portions of the spectrum in a single data file, 
striking the Center Button repeatedly displays pieces of the spectrum sequentially. 

The Left Mouse Button is used to identify species by comparing the observed wavelength 
obtained from the cursor location at the time the button was struck to a set of laboratory 
wavelengths located in various Look Up Tables. Immediately after the Left Button is de- 
pressed, a small menu of options is provided as shown in Figure 2. The three closest matches 
from the requested Look Up Table(s) are always provided plus any close matches from the 
researcher’s personal User Look Up table. There were no entries found from the User’s Ta- 
ble in the example shown in Figure 2. The investigator must make a selection by moving 
the Mouse up or down and striking the Left Button. The current position is noted with an 
arrow and the potential selection is highlighted. Notice that "NONE - Return to Spectra” 
and "Input Identity” are also options. If the latter is chosen, the program will present a 
number of questions to obtain information, including species name, its rest wavelength, and 
its oscillator strength. This complete species identification can be used as a temporary set of 
variables or can be permanently stored in the researcher’s personal table for subsequent access. 



Fig. 2 showing the entries from the Look Up Tables 
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The Right Mouse Button provides all of the remaining flexibility. Striking it brings up the 
UBerMenu of options shown in Figure 3. The novice is strongly encouraged to press this but- 
ton, if for no other reason than to see the selection. This is a relatively safe operation since one 
of the options is to "Take No Action”, just in case the Mouse Button was used inadvertantly. 


"Tot# No Action 
Qu& frettim itS £***♦ 

Delete lost Meoture men t 
input Remork on Lost Meot. 
Moke HARD-COPY 
Gat Next Spectra 
Liter Prog 1 — Dummy Pretent 
Liter Prog 2 - Dummy Pretent 
Stop for keyboard Input 
Click on Ret. Wavelength 
Adjust Plot Parameters 
Adjust Background 4 Error 


v * i • ' 1 1 * ■ 

T 

0 

-50 



1 1 . 1 • ■ < 1 -4 * ‘ 


1125.0 1125.2 1125.4 1125.6 1125.8 

Wove length (A) 


Fig. 3 showing the entries from the UserMenu 


The UserMenu shown in Figure 3 works like all menus appearing above the plot window. A 
selection is made by using the mouse to place the arrow over the desired entry and depressing 
the Left Mouse Button. While the functionality of most of the entries are straight forward, 
the "Click on Ref. Wavelength" deserves a few comments. If this option is invoked, the user 
iB solicited to move the cursor to any wavelength which he wishes "to declare” to be the 
rest wavelength and to click the Left Mouse Button. Now the abscissa coordinate, which is 
continually printed, reads in terms of velocity until the next graphical input. This capability 
is particularly useful when searching for a weak (1-2 sigma) feature after its doppler velocity 
has been determined from a strong line. 

Once a species has been identified using the Left Mouse Button and making a selection as in 
Figure 2, the program enters the continuum fitting subroutine. MSLAP provides an expanded 
plot with the individual data points highlighted. In this portion of the software, the Left Mouse 
Button is used to isolate regions of the spectrum that are believed to be featureless, regions 
that will be used to perform a polynomial fit. The user may specify 15 or less regions, which 
are marked on the plot with numbers from 1 to 2, 3 to 4, 5 to 6, and so forth (see Fig. 4). As 
before, the Right Mouse brings up various MENUs, while the Center Mouse Button signifies 
all regions of interest have been identified (go on to the next task). Menu options that are 
available in the continuum fitting routine include: 1) defining discrete continuum points, 2) 
identifying additional profiles for measurment and corresponding species identifications, 3) ad- 
justing the markers of the profile centers, and 4) bringing up the previous UserMenu of options. 

Once the featureless portions have been identified, the program calculates polynomials fits of 
order 1 through 7. Then, MSLAP, starting with a polynomial of order 1, sequentially tests 
the polynomials of increasing order, searching for the case where no statistically significant 
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improvement of the fit is realized by polynomials of higher orders. Specifically, the next higher 
and next, next higher orders are tested for polynomials up to order 5. The program uses F 
Distribution Tests with a 5% significance to make this choice. 


Then, a menu of options appears over the plot as shown in Figure 4. This menu is presented 
with MSLAP’s choice of polynomial being indicated. The user can over ride the automatic 
selection of the order number, 4 in this particular case. He can in fact over plot various con- 
tinuum fits until he is satisfied. Simultaneously, the following instructions plus a reminder of 
MSLAP’s choice appear in the text window: 


The recommended Order is 4 (highlighted), but you may select another. 

To assist in other choices: X is the observational difference in the 
reduced chi squares divided by the reduced chi square, which if larger 
than the theoretical F(l,n) indicates that going to the next higher order 
polynomial is justified statistically. T is similar to X except it is 
for comparison to F(2,n), an order that is 2 higher. The F Distributions 
are at the EX confidence level for: 180 points. 

Select the Order of the Polynomial 

- Other Polynomials may be examined before deciding. 

- Selecting order 4 implies use that polynomial. 


1125.0 1125.2 1125.4 1125.6 1125.8 1125.0 1125.2 1125.4 1125.6 1125.8 








Figure 5 shows several polynomial fits. For comparison, a polynomial with an order number 
that is one less is also plotted as a dashed line. (These were made using one of the options 
shown in Figure 4.) In this example, there is a significant difference between a polynomial 
of order 1 and 2, but not much difference between orders 2 and 3. There is also significant 
changes between orders 3 and 4, but not between 4 and 5 or between 5 and 6 (not pictured). 

A polynomial is finally determined by selecting the order number that is currently plotted. 
Notice that there are two entries that always indicate the order number currently plotted, but 
request either the use of the alternate integrating routine be used or a fixed sized integration 
window. The latter allows the user to specify some predetermined velocity interval over which 
the measurements are to be made. If this window is undefined, MSLAP will prompt the user 
for input. 

Finally, the cursors and Left Mouse Button are used to set the end points for the integrating 
region. (See §IV for cautions on setting this range.) The program then calculates the various 
moments and produces a plot such as the ones shown in Figures 6 and 7. The latter depicts 
a case where two profiles were identified as CN from the User’s customized table and were 
simultaneously measured. Up to 4 profiles can be measured at one time. The solid, horizontal 
lines near the bottom of the plot indicate those regions between the dotted vertical lines that 
were used in the continuum fit. The area of integration is shadded. It is possible, for example 
to integrate only part of the profile. 


Rel. Flux 
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./rocket/imaps Thu ^P r * 13:45:18 1991 



1125.10 1125.20 1125.30 1125.40 1125.50 1125.60 1125.70 1125.80 

Wavelength (A) 


Measurements in file: test.DTL , ^ eal S/N * 1 

Apparent S/N Ratio: 23.0 based on 185 points with a Noise Coherence Length of 3.0 

Continuum was fit with o polynomial of order: 4 0.25 0.53 

BG (bac^ r °un d ) was token to he: .! 2 -°. + /~« „ 


For F e II 1125.440 error contributions from BG: 0.0018 0.002 ^ 0.024 

Errors below are the Addition in Quodroture of the Background ond RMS*~Noise Errors 


Species 
Fe II 


Lab. Wave, f Obs. Wave 

1125.440 0.0110 1125.445 

Errors: 


EQW (A) 1st (km/s) 2nd (km/s/s) 
0.0323 1.245 72.094 

0.0023 0.433 4.905 


Fig. 6 showing the results for a single profile 




3873.0 3873.5 3874.0 3874.5 3875.0 3875.5 3876.0 

Wavelength (A) 


Measurements in file: test.DTL ^ t Real S/N — 548.6 

Apparent S/N Ratio: 548.6 based on 47 points with a Noise Coherence Length of 1.0 
Continuum was fit with o polynomial of order: 1 1.42 2.11 

BG (background) was taken to be: 0.0 +/ — 0-0 

For CN R(0) ....... 3874.608 error contributions from BG: 0.0000 0.000 . 0.000 

Errors below are the Addition in Quadrature of the Background and RMS— Noise Errors 


Species 
CN R(0) 

CN R(1) 


Lab. Wove, f Obs. Wave 
3874.608 0.0338 3874.628 

Errors: 

3874.000 0.0228 3874.025 

Errors: 


EQW (A) 1st (km/s) 
0.0077 1.586 

0.0002 0.550 

0.0020 1.915 

0.0002 1.174 


2nd (km/s/s) 
15.469 
16.462 
9.835 
22.954 


Fig. 7 showing the results for multiple profiles 


2: Measure Moments and Optical Depths 


This option is identical to the previous one, except in that example it creates an additional file 
test.TAU as well as test.DTL. Also, it plots the optical depths immediately after it displays 
the results of the moments (i.e. Figures 6 or 7). The optical depths associated with Figure 6 
are shown in Figure 8. 



15 



-20 0 20 

velocity 


Fig. 8 showing the optical depths across the profile 

The uncertainties, pictured in Figure 8 by open circles, only represent the systematic errors 
due to the uncertainties in the background and continuum placements. These data can be 
accessed later in Option 5: as profiles of the column density. 

3: Manipulate Tabled Data 

This portion of MSLAP enables the user to make customized tables of the measurements that 
were made in options 1: and 2: of the FIRST MENU. The investigator loads the dtl data 
structure by reading the .DTL file that was defined in option 1: or 2: of the FIRST MENU. 
He then can print an abreviated form of this dtl structure, create a ASCII format table of 
only the entries of interest, or has several sorting options, including reordering by Laboratory 
Wavelengths, by Observed Wavelengths, or by ION. 

The working text window should be expanded to include the entire screen in this portion 
of MSLAP. The table creation is designed to handle large tables and the full-screen window 
prevents text from rapping around, making it difficult to read. 

The routine is very generalized. Values supplied by standard MSLAP as well as those in the 
UserParameter (up) can be accessed and displayed in any order. Data from more than one 
file c an be combined into one master table. Measurements that are absent in one file, but 
present in another, are supplied with epsillons in the master table. However, all files that are 
to be combined must be sorted in the same fashion. 

Option 5 (to create a customized table in ASCII format) is particularly powerful. It allows 
the user to build the table, taking any number of pieces in any order and add these to an 
existing 2-dimensional character string. Parts of several data files can be included as well. 
Figure 9 shows the initial table at the start of Option 5. A small (the first 6 lines) portion of 
the table is shown at the top, while a list of options for adding to this 2-dimensional character 
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field is provided at the bottom. Initially only the species identifications, the wavelengths, and 
the oscillator strengths are included. A scale is provided immediately below the partial table 
to assist with the addition and deletion of columns of individual characters. The researcher 
then adds to, or subtracts from this 2-dimensional character field, using options 1 through 11. 


»»»»»»»»»»»» Partial Display of the TABLE <<<<<<<<<<<<<<<<<<<<<<< 
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Figure 9 showing the screen at the start of the creation of a customized table 
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Figure 10 same as Figure 9, except additional columns of data have 
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"figure 10 show# the development as the result of sequentially invoking options 4, 5, 6, 7, 8, 
and 9. Notice that an underline starting with the file name "test” extends to the right over 
all of this data. Once data from a new file is added the name changes and a new under line 
will continue from that point. Two data sets that have been combined from files: test2 and 
test3 as shown below. Notice that some minor editing of the file (such as the creation of true 
ellipses) is required. 


test2. 


test3 


Species 

Wavelength 

f 

EQW (A) 

ME (A) 

EQW (A) 

ME (A) 

Fe II 

... 1125.44 

0.0110 

0.0275 

0.0013 

0.0275 

0.0013 

Fe II 

... 1125.44 

0.0110 

0.0281 

0.0012 

0.0281 

0.0012 

Fa TT 

1125.44 

0.0110 



0.0262 

0.0012 

CM R(0) . . . 

... 3874.61 

0.0338 

0.0076 

0.0002 

0.0076 

0.0002 

CH R(0) . . . 

... 3874.61 

0.0338 

0.0169 

0.0009 

0.0076 

0.0002 

CM RfO} 

. . 3874.61 

0.0338 



0.0077 

0.0002 

CM RfO'i 

3874.61 

0.0338 



0.0077 

0.0002 

CM 

... 3874.61 

0.0338 



0.0169 

0.0009 

CM R(0} 

3874.61 

0.0338 



0.0879 

0.0227 

CM RfO’i 

... 3874.61 

0.0338 



0.1160 

0.0413 

vn ft v v/ • • • 

CM RfO} 

... 3874.61 

0.0338 



0.0076 

0.0002 

vn xv \ v/ / «9i 

CM R(l) ... 

... 3874.00 

0.0228 

0.0021 

0.0002 

0.0020 

0.0002 

CH R(l) . . . 

... 3874.00 

0.0228 

0.0020 

0.0003 

0.0020 

0.0002 

CH RCl) 

... 3874.00 

0.0228 

0.0020 

0.0002 

0.0020 

0.0002 


4: Edit Tabled Data 

This portion of MSLAP allows the user to perform basic editing functions on the dtl data 
structure. The routine is modeled after primitive line editors that were common in the days of 
the PDP-11 computers. Single key strokes control the editing functions. The user can move 
up or down by one line, move to the top or bottom of the data structure, insert or delete a 
line, or change a line. A line is defined currently as a complete set of measurements for a 
single profile, including any comments and any values in the user parameter. A small window 
is opened on the right side of the screen and a complete list of the functions is printed as a 
reminder to the user. 

The changes become permanent if the userB updates the file on disk. Otherwise, the researcher 
may Kill/Quit the session at any time without making any additional changes since the last 
time the disk file was updated. 

When the option to changes the values in a single "line” is used (i.e. pressing a "c"), all values 
are defaulted to the existing values before entering this mode. In this manner, the investigator 
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only has to strike the j return^, key repeatedly until he reaches the value(s) to be changed. This 
feature reduces the probability that typographical errors will be introduced while attempting 
to make minor changes. 


5: Analyze Column Densities vs Velocity 

In this portion of MSLAP, the investigator cycles through a series of MENUs over the plot 
window. The MENUs are ordered in the following hierarchical sequence to determine: 1) what 
to do next, 2) the species, 3) which profile (identified by wavelength), and 4) plotting symbol 
or style. After each profile has been plotted the researcher can then get additional profiles or 
make adjustments to the one just plotted. Figure 11 shows an example of the Fe II column 
density as a function of velocity for several profiles. 


Input Fila: loutst.TAU Fa II May 20, 1990 



— 

1081.680 

1.196 

— 

1125.440 

1.093 

• — - 

1144.946 

2.080 


Fig. 11 showing the output from option 5: 


Figure 11 shows the column-density results of two weaker profiles of Fe II, which are in good 
agreement with each other, plotted as continuous lines. The other over plot (solid line with 
dots) is for an intrinsically stronger absorption line. The disagreement indicates the presence 
of narrow, unresolved, saturated structure. Figure 11, also showing two of the different plot 
options, has error dots above and below for the strong line. A key of the plot symbols is 
provided at the bottom of Figure 11, showing the wavelength and log(fA) values. 
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6: Compare Data to Curve of Growth 


The Curve of Growth routine is very similar to the one to analyze column densities as a 
function of velocity. The investigator configures the plotting through a series of MENUs that 
appear over the plot window. The first species to be plotted is selected, its plotting symbol, 
and the number of theoretical curves. All of these can be changed dynamically at any time to 
assemble a figure as shown in Figure 12. The software then presents the data and the inves- 
tigator moves these data points horizontally by using the cursor to note starting and ending 
positions. Once satisfied with the fit, the researcher can then get additional data points to 
fit, or has a host of other options by using the right mouse button. 



Obs. doto from file; test.DTL 
Theo. COG: 0 I 1302.169 


Moy 20, 1990 

r*0.CK86 CAMMAw5.72E08 0.0 5 3 1 601 


too NfCN R 
too N(CN R 
to$ N(Fe ll 



10.9610 

10.5549 

13.1610 


Fig. 12 showing the output of the Curve of Growth 


IV. Comments on Correct Statistical Techniques 

The subject of Correct Statistical Techniques is one that requires numerous pages to be com- 
plete. Unfortunately, time pressures on the author do not permit an adequate treatment 
in the current volume. It is hoped that future updates of the MSLAP Documentation will 
contain substantially more topics and details on the various aspects of using good statistical 
techniques as they apply to the MSLAP analysis package. There are several issues, however, 
that must be stated, even if only briefly. 

1) Getting a good fit of the continuum is particularly important. The uncertainties in all of 
the measurements are determined from the residuals of the real data minus the fitted con- 
tinuum for each point in fitting region. If there is a systematic error produced over some 
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portion of the spectra, this may manifest itself in the form of unrealistically large values of 
the random uncertainties. To avoid this difficulty, one might be tempted always to fit with a 
polynomial of the highest order available on the assumption that such a fit should be equal to, 
or better than that from a lower order polynomial. However, such action is not without risk. 
If a polynomial with too high of an order number is selected, the continuum fit occasionally 
may be ill behaved, especially over the regions where it is not constrained. In other words, 
the fitted continuum may not represent the real continuum at precisely those wavelengths 
where it is needed mo8t, over the absorption profile. The author did not have to search very 
far through his own data to see several examples of this problem. Figure 13 shows a first 
(solid line) and fifth (dashed line) order fit to the same Fe II profile. The first order fit is 
the appropriate choice based on the F Distribution Tests for selecting polynomials of various 
order. The first, second, third, and fourth order fits are almost indistinpishable from each 
other. If one were blindly to use the highest (5th) order available, the equivalent width would 
be over estimated by more than 1 sigma of the best estimate. Other profiles exhibited far less 
dramatic continuum shapes, but just as severe differences in the measured values. 



Fig. 13 showing different polynomial fits. 

2) In MSLAP, the user specifies which regions of the spectra that he believes to be free of any 
spectral feature for the purposes of performing a continuum fit. The subjective responses of 
the researcher can make it very easy for him to avoid preferentially any 2-, or 3-sigma, de- 
viations in choosing these regions. This action, of course, skews the uncertainty calculations 
towards smaller errors. While 3-sigma events are uncommon, these still have an impact on the 
total uncertainty since the contribution to the sum from each data point is quadratic in nature. 

3) In the standard integrating program, the user specifies the starting and stopping points. 
As with issue 2, there is a human temptation capable of introducing systematic errors. In the 
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ra p* of absorption features, many researchers choose the end points of the integration to be 
the places where the spectrum crosses the fitted continuum closest to the center of the profile. 
This action leads to a predictable, systematic over estimation of the equivalent widths and 
second moments (Joseph 1989, PASP, 101, 623). 


4) There are a n um ber of common oversights. For example, the background uncertainty in 
many applications is the dominant source of error, but it is often completely ignored by the 
investigator. Likewise, the researcher needs to pay attention to the coherence length of the 
noise. The easiest way to visualize this coherence length is to consider first a Bpectra produced 
by an instrument in which each pixel is totally linearly independent of all the others. These 
data would then have a coherence length of unity. If, instead, the value in each pixel in the 
instrument was equally dependent on the values of its adjacent pixels on either side, or if the 
previous data have been smoothed by a 3-point running boxcar, then the coherence length 
would be 3. For various reasons of data handling and instrumental affects, most data to be 
analyzed usually have a coherence factor greater than unity. 


V. Customizing MSLAP 

MSLAP makes use of a IDL feature where subsequent calls to compile a routine with an iden- 
tical name overlays the previous. The routines have been organized into several files in part 
according to the frequency in which each is expected to be modified. This organization is an 
attempt to minimize the amount of code the typical researcher has to wade through in order 
to make his desired changes. For example, 90% of all alterations are expected to occur in the 
file mslap.pro, containing less than 350 lines of code compared to the more than 4,000 lines 
for the entire MSLAP package. All source code is listed in Appendix A, including a table of 
contents with brief description of each routine. 

As already stated, most modifications to MSLAP are expected to be implemented inside 
mslap.pro. Before running MSLAP, each user copies the mslap.pro file to his own directory. 
This file contains a number of dummy programs that have been commented out, but with 
instructions showing how to install altered forms of these routines. Thus, many investigators 
ran have their own customized version of MSLAP without having to keep a complete copy 
of all of the source code. In fact, occasionally a researcher will have two or more customized 
forms of MSLAP in his account. 

The rest of this section is divided into several specific applications. First, a detailed descrip- 
tion of the mslap.pro file is given. Next, the data-getting routines and the file dgets.pro is 
presented. Finally, the structure and use of the Look Up Tables is discussed. 


a. The mslap.pro File 

The file mslap.pro contains many useful dummy programs that have been broken out of the 
main software package. Some of these are real dummy programs, allowing the user to cus- 
tomize MSLAP merely by inserting a few lines of code, while others are actually comment 
fields for the purpose of being able to overlay different individualize routines. In addition, 
there is a routine at the beginning of mslap.pro called: usermenu, which provides the char- 
acter strings used to print the primary MENU of options found above the graphic plot when 
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the right mouse button is pressed. This MENU is the principle means of branching once the 
main portion of MSLAP is running. Occasionally, a researcher will develop a mental block 
regarding the meaning of a question or menu option, which can be especially problematic if 
the investigator uses the software infrequently. In these cases, the user is invited to modify 
these character strings to make them more meaningful. The functionality, however, is defined 
elsewhere and remains unchanged. 


There are 3 dummy routines, userprogl, userprog2, and userprog5, reserved for the user to 
create specialize programs. These routines are accessed during the execution of MSLAP by 
first pressing the right mouse button and then selecting the appropriate menu option. User- 
progl and 2, either of which can be executed from Options 1: or 2: in the FIRST MENU, are 
particularly useful in cases where the set of calculations are not always performed. If a set of 
specialized calculations iB to be performed each time a profile is measured then auto_int_sav 
may be the more appropriate dummy routine to use. The latter is called every time the various 
moments of the profile are calculated. In a similar fashion, userprogS is used for inserting cus- 
tomized calculations during the portion of MSLAP that analyzes the Optical Depths/Column 
Densities (Option 5: of the FIRST MENU). 

Furthermore, there are 2 routines, plotlabl and plotlab2, that are automatically called during 
the plotting of spectra and expanded spectra, respectively. These routines can be used, for 
example, to place additional labels on the plots. 


b. The dgets.pro File 

There are 5 slots available in MSLAP for the data-getting routines used to read the files con- 
taining the input spectra. These routines, called dgetl, dget2, dget5, are formally outside 
of MSLAP, and therefore, are the sole responsibility of the user to write or obtain from some 
other source. DGET routines for various data sets, however, are being written all the time 
and these contributions will be continually added to the library. While MSLAP only supports 
up to 5 dget routines at any one time, each investigator chooses the 5 that best meets his needs. 

These files are "linked” into MSLAP by issuing indirect compile statements in the mslap.pro 
file. For instance, the syntax to compile the dget2.pro file is: ”@dget2.pro”. The user need 
only include as many (5 max.) dget routines as he needs. For convenience, a file called 
dgets.pro, containing 5 sample dget routines in a single file, has been included in the library. 
An indirect compile statement of this file avoids including 5 similar statements in mslap.pro 
and minimizes the number of file names in the researcher’s directory. 

To install a dget routine, the user must place the indirect compile statement in mslap.pro as 
already mentioned and must include the following statements at the beginning of the program: 

if 0 th#n b#gin 

ap . dg«t ■ 1 

return 

•nd 

These statements form the mechanism that MSLAP uses to sense that a valid dget routine 
is present and Bupercedes the dummy one that is otherwised supplied. As an added feature, 
the string variable called "ID” can be set to a brief, descriptive text, indicating the nature of 
the routine. This variable will also appear as part of the print-out of available data-getting 
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routines. For example, adding the line: 


U 5 » ***** This i* mj ftTorit# routin# 

to dgetS.pro just after the procedure definition statement will produce the following entry 
during the set up phase of MSLAP: 


DATAGET Option*: 

1) IUE Standard CO Filaa for Hi^h-fU*. Sp*ctra (DISIGET) 

2) IDT B#iag Uaad 

3) ASCII Formit of ¥*T#l*ngth-Sp#ctra-Quality-BG 

4) 1034-El on ant Stand-Alona Data 

5) **•* Thi* i* ny favorit* routin# **** 


Vhich ona would you liAa ? 


Finally, the user should place in the parameter mp.FNAM, the file name being opened by 
the dget routine. MSLAP displays this parameter as the input data source above some plots. 
Other features include: mp.order, mp.CAM, mp.bg, and mp.bgerr, which indicate echelle or- 
der number, the camera number, the background level, and the error in the background eve , 
respectively. Note: all parameters in the three main data structures are global in nature and 
thus, can be accessed from any routine. See Appendix B for a complete listing along with 
descriptions of functionality. 


c. Look Up Tables 

MSLAP creates a file called ulut.tab in the working directory of the investigator when he first 
requests to store an species that he has identified. If this file exists, MSLAP searches it for 
close matches every time an identification is requested. MSLAP defines a close match as being 
within 10 pixels or 0.1 A which ever is smallest. The value of the search agreement can be 
selected by setting: mcntrl.wtol. Note: removing the file ulut.tab destroys all user-ID entries. 

Software exists to convert any look up table that may exist in ASCII format into one that 
can be used by MSLAP. For the time being, interested parties should contact the first author 
regarding this utility. 


VI. Installing MSLAP (for the system manager) 

Obviously, the location of the software and tabled information can be organized in a number 
of ways. The present discussion, however, will only describe the simplest method of installing 
MSLAP. Note that MSLAP is written in the IDL language and must run on a workstation 
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supporting IDL. 

First select or create a directory to hold the MSLAP programs and change to that directory. 
The entire package requires about 500 kilobytes of space. Next create a subdirectory called 
*tabdata" to hold the tabled information. Standard MSLAP expects the look-up tables to be 
located in a subdirectory of the directory holding the programs. All files containing tabled 
data have a ” .tab” extension. The rest of the files should be placed in the parent directory. 
Make sure that all of the files have read-only permission. IDL programs are always compiled 
in real time from the source code. Thus, it is unnecessary for individual users to have either 
execute or write permission. 


Edit the files: main.mslap, master. aux, and mslap.pro so that all occurrences of the string 
sequence: /u/clj/mslapdir are replaced with a string sequence appropriate for the path to the 
parent directory holding the MSLAP software. The new string should be sufficiently complete 
that users ran access this directory from any directory where they store their programs and 
data. Note: be sure to keep the symbol since this is a linking or call to compile command. 


Turn to the section on Getting Ready (§11 of this Manual) and write in the directory path 
where you have just installed the MSLAP software. Users will want to copy the mlsap.pro file 
from this directory. 

All plot co mm ands that axe hardware specific are located in a single file called plotconfig.pro. 
In standard MSLAP, this file is designed to handle automatically either sunview or X windows 
and to open one of several types of graphical or text windows, each with a specific location 
and size. Some minor adjustments to the values specifying size and location, therefore, may 
be necessary from one machine to another. The plotconfig.pro program also configures the 
hardware so that IDL plot commands produce PostScript plot files, which serve as hard 
copies of graphical output appearing at the console. The user of MSLAP can request that 
these "plots” be sent to the laser printer from inside MSLAP. If the default laser printer is 
not a PostScript printer, the line of code with the command: 

Bpawn,’lpr temp.ps’ 

must be changed so that ”lpr” is changed to ”lpr -P[name]”, where name is the device support- 
ing PostScript. If your computer supports another device such as Hewlett-Packard Graphics 
Language (HP-GL) which IDL is equipped to handle, then you will need to consult the IDL 
manual for the appropriate commands to substitute in plotconfig.pro. Only the code inside 
plotconfig.pro needs to be altered, however. If you do not have a hard copy unit that IDL rec- 
ognizes then you should comment out those portions of plotconfig.pro. Note: the appearance 
of a semicolon causes IDL to consider the remainder of that line to be a comment field. 


The file dgets.pro contains a set of 5 data-getting routines, all of which are formally outside 
of MSLAP. Actually, it is the responsibility of the researcher using MSLAP to write or substi- 
tute as necessary various dget routines so that his data can be read into MSLAP. However, if 
you represent a guest-user facility, you may wish to substitute several dget routines with ones 
that are suitable for your specific data formats. The new combination can either replace the 
dgets.pro file or an alternative file can be generated. Then a guest user can simply copy one 
of these files to his own directory and run MSLAP with minimal amount of effort to get started. 

Finally a word of caution. MSLAP is designed so that customized routines can be inserted 
with minimal difficulty. The way that this is accomplished is to make use of a feature of IDL 
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where a subsequent compilation of & routine with the same name overlays the previous. A 
number of dummy routines exist as place holders in the event that no user- written programs 
•re supplied. One may be tempted to consolidate the UBe of indirect compilation calls (e.g. 
the use of ©program .name), but this may eliminate the modularity of MSLAP. The sequence 
of compilations is im portant and you should be careful not to change this order inadvertently. 


VII, Trouble Shooting 

The most common problem faced by the user of MSLAP appears to be either the failure to 
register a cursor reading or the registering of multiple readings from a single stroke of the 
mouse. For this reason, MSLAP normally attempts to reward the investigator with some type 
of output every time that he inputs a cursor location. The continuum fitting routine can be 
the most problematic since many cursor locations are read to set the intervals to be used for 
the polynomial fit. During this activity, the routine marks each selection with a colon and 
a number beside it. The regions to be fit, thus, are from 1 to 2, 3 to 4, 5 to 6, and so on. The 
researcher is encouraged to monitor these numbers while selecting the ranges. If two numbers 
are superimposed or if one is missed, the continuum fitting routine should be started again. 
If improper registering of cursor readings becomes severe, then IDL "WAIT” statements will 
have to be included. 

The data-getting routines (dgetl, dget2, ..., dget5) are formally outside of MSLAP. Some un- 
expected behavior in MSLAP can occur if these routines do not make use of all of the required 
variables in the MSLAP-Parameter (mp) structure. For example, when the user requests new 
data in MSLAP, the software goes to one of the dgets provided by the researcher. That dget 
routine could solicit a new file name and read new spectra, but not place the updated file 
name in the mp.FNAM parameter. Then subsequent plots will not reflect the true input file 
name. Note that in some applications such as echelle data, it is desirable to configure the 
dget routine so that subsequent calls to dget return spectra from the next order and not open 
a new file. For other types of data, it is more natural to have a single spectrum for each file. 
Hence, the responsibility to return the necessary parameters to MSLAP has been left to the 
user. (See the section on the dgets. pro file.) 


APPENDIX A - The Source Code 


Routine Name 


Description 


Page 


details .pro 

flagset .pro 
graphs .pro 

set up. pro 
contim.pro 
store .pro 

f store. pro 
expand. pro 

the main prog. 


In the file: main.mslap 


The MAIN Program for Options #1 i #2 from 

the First Menu 30 

A set-up routine used by details. pro 33 

Does the plotting/labeling for details. pro, 

but does not plot expanded portions 34 

A set-up routine used by details. pro 36 

The continuum fitting routine 38 


Stores the result from details. pro in the dtl 
data structure and on disk (see fstore) . . . 

Stores optical depths to disk (see store) . . . 

Extracts a suitable portion of the spectra 

for an expanded plot for contim.pro ..... 

Provides the FIRST MENU and calls, 

details. pro, posto.pro, edatdtl.pro 
mantau.pro, or cog. pro 


dget (1-5) .pro 
curf it .pro 
eqsol .pro 
ucursor .pro 

dataget .pro 


In the file: master. aux 

The Dummy Versions as place holders 

Calculates a single polynomial fit 

Special routine called only by curfit.pro 

Provides cursors with special features. Also, 

provides Continuum Menu Options 

An interface routine between MSLAP and dget. pro ... 


52 

53 

54 

55 
58 


compare .pro 
userids .pro 


In the file: compare. pro 

- To find Look Up Table (LUT) entries, main driver .... 60 

- To find User LUT entries or make new ones 62 


posto .pro 
srt .pro 
publ.pro 
mkcusttab .pro 
appnd .pro 
next st ring . pro 


In the file: posto.pro 

Organizes the tabling of the data base 

Sorts the data for posto.pro 

Lists (Publication Form) the data for posto.pro 

Makes Customized Tables in ASCII format 

Called by mkcusttab.ro to append entries .... 
Called by APPND & MKCUSTTAB to get the 
size of an entry and to convert to string . . . 


64 

65 

66 
67 
71 

73 


In the file: edatdtl.pro 

edatdtl.pro - Main program for editing the "dtl" data structure ... 74 

getsdata.pro - Get single set of measurements '6 


In the file: mantau.pro 


79 

86 


t^ay .pro 


intgrt .pro 
equivw.pro 


- Manipulate (tau) optical depths 

- A set up routine (GetReady) . . . 


In the file: intgrt .pro 

- The primary integrating routine 

(does significant plotting) 

- Calculates the moments plus optical depths . . . . 

(originally just an Equivalent Width routine) 


88 

92 


In the file: plotconf ig .pro 

plotconf ig . pro - Sets up the hardware for plotting functions 
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cog. pro 
fireup.pro 
lines .pro 
grab. pro 
plot em. pro 
whichcog.pro 


In the file: cog. pro 

- Curve-Of-Growth primary routine 

- A set up routine • 

- Gets and plots the theoretical curves 

- Grabs the observed data points to be fit 

- Plots the data from GRAB. PRO with various symb 

- Produces MENUS to input # of theo. COGS 
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100 

102 

103 

104 

105 


Modular Spectral Line Analysis Program 
.S.L.A.P. version 1.0 


copyright (c) 1991 by Charles L. Joseph and Edward B. Jsnkins 
All rights reserved. A license may be obtained from the 
first author or from an authorized distribution center. 


software has been distributed through the GHRS Data Center. 


THE FOLLOWING OPTIONS ARE AVAILABLE: 


1: Measure Profile Moments Only 

(Equivalent Width, Profile Centroid, etc.) 
2: Measure Moments and Optical Depths 

3: Manipulate Tabled Data 

4: Edit Tabled Data 

5: Analyze Column Densities vs Velocity 

on Option t2 Data 

6: Compare Data to Curve of Growth 

10: TO EXIT FROM THIS PROGRAM 


COPYRIGHT (c) 1984, by Charles L. Joseph 

COPYRIGHT (c) 1991, by Charles L. Joseph k Edward B. Jenkins 


Released: Ol-February-1991 

Alteration History: None Bince release date. 


PRO DETAILS, sff ,timedate,libr, up 

; ******************************* DETAILS. PRO ****************************** 

P 

; To control the measurement of equivalent widths or integrated fluxes 
! BY CHARLES L. JOSEPH 09-SEPT-80 11 -MAY-90 

; sff 0 *> Do not do Optical Depths 1 *> Do Optical Depths (TAU) 

; timedate Variable containing juBt the date. 

• up User Parameter - 30 element floating vector t which is global. 


up is rsssrvsd for the exclusive UBe by the user. 


f 


; Set up dete structures — "mp" ■> MSLAP Parameter "dtl" »> for .DTL files 

■p ■ ^ MS LAP , dget: 0, cntrl: -1, order: 0, CAM: 0, SMO: 0, cohfac: 0., bg: 0.,$ 
bgerr: 0.0, window: 0.0, DTY : 0, STAR: » ». FI AM: » poly: 0, SIR: 0. > 
■cntrl ■ { MSLAPc, I: 0, date: ”, intopt: 1, IHE: 0, ESAV: 0.0, eO: 0., el: 0.,$ 
wtol: 0., e2: 0., libr: ”, »tot: 5, fin: 0., f2n: 0., fdist: fltarr(3,18) > 


fdist « [0..1., 2., 3., 4.. 5., 6.. 7.. 8.. 9., 10.. 15.. 25., 50., 100., 500., 1000., 
■cntrl. f diet (0,*) ■ fdist 

fdist • [ 161 . , 161 ., 18. 51, 10. 13, 7. 71, 6. 61, 5. 99, 5. 59. 5. 32, 5. 12, 4. 96, 4. 54, 
4.24,4.03,3.94,3.86,3.85,3.84] 

■cntrl. f diet (1,*) * fdist 

fdist « [200. ,200. ,19.00,9.55,6.94,5.79,5.14,4.74,4.46,4.26,4.10,3.68, 
3.38,3.18,3.09,3.02.3.00,2.99] 

■cntrl .fdist (2,*) * fdist 


10000 .] 

$ 

$ 


■cntrl. date ■ timedate 
■cntrl. libr * libr 


nulld ■ { noda, el: bytarr(15), iaf : long(O), el: 0., f: 0., oel: 0., eqw: 0., 
■e: 0., fm: 0., fme: 0., sm: 0., sme: 0., com: bytarr(lO), up: fltarr(30) > 
dtl ■ replicate ({ noda }, 200) 
sdata * nulld 

■data * replicate({ noda }, 4) 

COMS « ' . . . ’ 


Make array for storage 
Make working copy 
and for multiple meas. 
Make dummy comment 


SETUP , ML ,RS , FF , dtl .nulld .nside ,mp ,mcntrl .up 


; Set everything up A 


SEIT : dat aget , WA , SPECT , J J , IHDR , FF , mp .mcntr 1 , up 
IF mp. cntrl eq -99 THEN goto, FINISH 
WW * 0.0 

!x. range * [ 0. , 0. ] 

!y. range * [ 0. , 0. ] 


; Get next spectrum. 

; How to get out. 

; Set for AUTO Scaling. 


IT: 


GRAPHS, SPECT, WA, IHDR, JJ.mp.mcntrl, up ; Plot routine 

np . cntrl * 1 * Reset CNTRL flag 

■data « replicate({ noda >, 4) I Reset multiple meas. 

print,’ ’ 

print, 'Left Mouse to Locate Feature to be Measured — Need HOT be centered’ 


print , ’Middle Mouse to GO ON, get new data’ 

print, ’Right Mouse to bring up UserMenu of other options.’ 


MORE: TVCRS, 0.5, 0.5, /NORMAL 

ucursor , IX , IY , opt , iopt , 0 ,mp 
; control options. 


Set action to graph. 
Get graphical input A 


; »»>>>>>>>»»> Master Branching Section <<<<<<<<<<<<<<<<<<<<<<<<<<< 


IF opt «q 9 TEES begin * STOP ft goto, MORE ft end 
IF (opt «q 6) OR (!ERR EQ 2) THEN GOTO, IEXT 
IF opt «q 3 THEN BEGIN 


Stop for keybrd input, 
lext spectrum to plot. 
Delete last measurement 


I ■ mcntrl. I 
dtl(I) » nulld 

PRIST,’ Lina ’,1,’ Las baen DELETED! 

mcntrl . I*mcntrl .1-1 

IF mcntrl. I LT 0 THEB mcntrl. 1*0 


’ .STRING (7B) ; Print notice. 

Set I to over write. 
Can't so negative. 


ESDIF 

IF opt aq 4 then begin 
•data * dtl (ncntrl . I- 1 ) 
read.COMS 

•data. con • byte (string(COMS, format*’ (alO) ')) 
dt 1 (ncntrl. I- 1) ■ »data 
•data ■ nulld 


; Comnents? 

Get info on last neas . 
Read comment. 

Put it into byte array 
Update master data set 
Set existing to null. 


goto, Storit 
•ndif 

IF opt eq 2 THEN GOTO, FINISH 
if opt le 4 then goto, MORE 
if opt eq 7 then begin 
print , ’ * 

usrprogl ,WA .SPECT, JJ ,mp ,up 
•ndif 

if opt eq 8 then begin 
print,’ ’ 

usrprog2 , WA , SPECT , J J ,mp , up 


; End session option 
; Rem. , Del . .Nothing 
; Call User Prog. #1. 
; Just to flush cursor. 

I 

; Call User Prog. #2. 
; Just to flush cursor. 


andif 

if (opt eq 7) OR (opt eq 8) then goto, IT ; UserProgram called, 
if opt eq 5 then begin mp. cntrl * -10 ft goto, IT ft endif ; Make hardcopy, 
if (opt eq 11) OR (opt eq 12) then goto, IT ; Adjustments only. 


; >>>>>>>>>>>>>>>>>> End Master Branching Section <<<<<<<<<<<<<<<<<<<<<<<<< 


Begin Normal Flow of Routines 
"The Guts of this Program" 


MEAS : WW* IX 

J*0 

J « TOTAL (WA LT WW) 

COMPARE , WW ,RS , sdat a ,mp .mcntrl 
IF mp. cntrl eq 0 THEN GDTD.IT 
!x. range ■ [0. ,0.] 

!y. range * [0. ,0.] 


Save wavelength. 

Set up ft get closest 
index for EXP AND. PRO 
Get LUT entry. 

Sot a good ID ? 

Reset to auto scaling 
for the plots . 


EXPAND , SPECT , ESP , WA . EWA , JT , JE , J , ns ide 


Data for expanded plot . 


Then, Fit a continuum. 

C0NTIM .ESP ,EWA , CNT , WA , J J , XF , XI , JE , XLIMIT .mp .mcntrl .up , sdata .mdata . $ 
nulld, RS, ns ide 
IF mp. cntrl eq 0 THEN begin 


; Flag set to abort . 


; Reset to auto scaling 
; for the plots. 


!z. range « [0. ,0.] 

!y. range ■ [0. ,0.] 
goto, IT 
•ndif 

IF mp.cntrl gt 2 THEN GOTO. NEXT ; Flag for next spectrum. 

TAU ■ fltarr(6,8)+100. » dummy just in case 

IiTGRT,ESP,EWA,CNT,TAU,JE,XF,II,ILIHIT,mdata.mp,mcntrl,up ; Calls EQUIV.PRO A 
plots results. 


!x. range ■ [0. ,0.] 

!y. range ■ [0. ,0.] 

> 

Storit: IF mp.cntrl eq i THEN STORE, WV, ML, dtl.mdata.mp.m 
if (mp.cntrl eq 1) and (sff eq 1) then begin 
f store .TAU.sdata ,mp .mcntrl ,up 
end if 

IF Copt ne 2) AND (mp.cntrl ne -99) THEN GOTO, IT 

f 

FINISH: wdelete.O 
wdelete, 1 

if sff eq 1 then close, 5 

close, 1 » Close 

openu , 1 ,mp . STAR* * . DTL ’ 
writeu, 1 ,dtl 
close, 1 


Reset to auto scaling 
for the plots. 

icntrl.up ; Store .DTL 
Optical depths too? 

If so, plot A store. 

Do some more or go 
back to MSLAP 
Remove plot windows. 
Remove plot windows. 
Close .TAU file. 

DTL file. 

Open again to catch 
any last minute c h a n ges 
Close .DTL file again. 


RETURN 

END ; DETAILS 


PRO FLAGSET.FF.mp.up 

;************************* FLAGSET.PRO ************************************* 

I 

; by Charles L. Joseph June 29, 1990 

» 

; This procedure determines which look up table and data-getting 
; rountine are to be used. 


FF 

mp 

up 


flag for which DGET routine to use 
MSLAPparmeter - a structure 
UserParameter - 30 pt floating array 


******************* ********************************************************** 
DTY«1 

DETR: for k « 0,3 do print,* * ; print lookup tables 

print , 

print , ’Available lookup tables are:’ 
print , * ’ 

print,’ 1) Interstellar Lines (Morton and Smith 1973 plus updates) 

print,’ 2) Molecular Hydrogen, HD, and CO* 

print,* 3) Options 1 A 2 Combined’ 
print , ’ 4) Hot-Star Lines ’ 
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print,’ 

for k«l,2 do print,’ ’ 

READ, 'Which on* would you like ?’,DTY 
IFCDTY LT 1) OR ( DTY CT 4) THEM GOTO.DETR 


i 


; find out which one 
; no such table 


■p .DTY « DTY 

RFL: up » fltarr(30) ; define User-Param. 

iok « intarr(6) 
for k ■ 0,3 do print,’ ’ 

print, ’ 

print , ’DATAGET Options : ’ 
print , ’ ’ 

■p.dget ■ 0 

dget 1 , FI AM , w , spect , q , ihdr , ID ,mp , up 
if ap. dget *q 1 then iok(l) ■ 1 
print,’ 1) *,ID 
■p.dget ■ 0 

dget2 , FMAM , w , spect , q , ihdr , ID ,mp , up 
if up. dget eq 1 then iok (2) * 1 
print,’ 2) ’.ID 

■p.dget * 0 

dget3 , FMAM , w , spect , q , ihdr , ID ,mp , up 
if Bp. dget eq 1 then iok (3) ■ 1 
print,’ 3) ’,ID 

ap. dget * 0 

dget4 , FMAM , w , spect , q , ihdr , ID ,mp , up 
if mp.dget eq 1 then iok (4) ■ 1 
print,’ 4) ’,ID 

mp.dget * 0 

dget5 , FMAM , w , spect , q , ihdr , ID ,mp ,up 
if mp.dget eq 1 then iok (5) ■ 1 
print,’ 6) ’,ID 

print, ’ 

FF « total(iok) ; How many valid progs? 

if FF eq 0 then begin ; If none, complain, 

print .string (7B) , ’There are no valid data-getting routines loaded.’ 
print, 'At least 1 valid routine must be put in the file: dgets.pro* 


; Sense which data- 
; get routines are 
; present and show 
; options. 


FF « ’ ’ 

read, 'Press <Retum> to exit’.FF 
retail 

end if 

for k*l,2 do print,’ ’ 

READ. ’Which one would you like ?’,FF 
if (FF It 1) OR (FF gt 5) then goto, RFL 
if iok(FF) ne 1 then goto, RFL 
if iok(FF) eq 1 then mp.dget * 1 
RETURN 
END 


; Pause so user realizes 
; the trouble t abort . 


Find out which one. 
Valid Choice? 

Do this one exist? 
Insure flag is set. 


PRO GRAPHS , SPECT , WA , IHDR , J J ,mp .mcntrl ,up 
;******************************** GRAPHS. PRO 


******************************** 
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TO PLOT THE SPECTRUM. By Charles L. Joseph Sept. 1980 lov, 1990 


SPECT 

VI 

IHDER 
I HD El 
JJ 

FMAM 

■cntrl 

■P 


up 


flux vector 

wavelength vector 

header vector 

plot or replot status 

data quality vector 

data file name 

structure (see DETAILS. PRO) 

KSLAP structured variable: mp . cntrl 

mp . order 
■p.cam 


control flag 
order number 
camera number 


UserParmeter - 30 pt floating array, exclusive for user. 


What-to-do HardCopy Flag. 

Config for Hardopy. 
Input: ’+FNAM 


****************************************************************************** 

iy.margin(O) ■ 4. 

FHAM - mp.FUAM 
star ■ mp.STAR 
kwtd ■ 0 

if mp. cntrl eq -10 then begin 

kwtd * smenuC ['Copy?’ , ’Yes’ , ’Ho*] ,title=0, init*2) 
if kwtd eq 1 then plotconf ig, i , ’ ’,’ ’, -l.kdev,” ; 
endif 

tlO ■ ’Output: ’+star+’ M.S.L.A.P. 

tlO * tl0+’ ’+mcntrl.date 

xxxs ■ fltarr(2) 
yyys ■ fltarr(2) 
sz * size(WA) 
jcnt * 1 

if sz(l) gt 16 then jcnt * f ix(sz(l)/16) 
if mp. cntrl eq -I then plotconf ig.O.tlO, ’ ’.l.kdev,” 
ymn * min (SPECT, J5, MAI*ymx) 
if ymn gt 0 then ytst * ymn ELSE ytst = ymx 
if ytst It 0.75*ymx then ymn * 0. 
twit * max(WA,jmx, MIN=nnn) 
vmn * mp.bg - mp.bgerr - 0. l*(ymx-ymn) 
if ymn eq 0. and vmn gt 0 then yt s t * 1 ELSE ytst=0 
if ymn le vmn then vmn = ymn - 0 . 1* (ymx-ymn) 
if ytst eq 1 then vmn * 0. 
ymx * i.l*ymx 

xxxs(O) * xmn k xxxs(l) * xmx 

yyys(0) * vmn A yyys(l) * ymx 

if !x.range(0) OR t x. r ang e (1) OR !y.range(0) OR !y.range(l) then begin 

if !x. range (0) OR !x.range(l) then begin ; Use these values instead? 

T-win e !x. range (0) ; If on® is nonzero *> user 

xmx * !x.range(l) ; wants to over ride, 

endif 

if !y .range(O) OR !y.range(l) then begin ! Use these values instead? 

vmn ■ !y.range(0) 
ymx * ly.range(l) 


; Set up graphics device. 
Find min and max flux. 

For bottom of plot use 0 
unless featureless. 
Min k max plus max (index) 
Get min. needed for BG 
Prevent double jeopardy? 
Less than min. spect? 
Correct double jeopardy. 
Adjust for more space. 

Put limits into small 
over plot array k 
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•ndif 

•ndif 

plot , xxxs , yyy s , /NODATA , rtitle- ’Wavelength (1)', ytitl*»’R«l. Flux', $ 
xrange«[xmn,xmx] ,yTange*[vmn,ymx] ,XBtyle*l,ystyle*l 
if mp.SHO than begin 

OPLOT,WA,SMOOTH(SPECT,3) , $ ; Plot Mioothed flux. 

Ititle* ’Wavelength Angstroms’ ,Ytitle*’Rel. Flux’ 

•ndif alas begin 

OPLOT, WA.SPECT, $ ; Plot unsmoothed flux. 

It itle* ’Wavelength AngBtroms’ ,Ytitle*’Rel. Flux’ 

•ndelse 

OPLOT.WA.JJ ; Plot data quality, 

if (J5 gt jcnt) than xxxs(O) * WA(J5-jcnt) ELSE xxxs(0)**WA(0) ; BG range, 

if (J5 It jmx-jcnt) than xxxs(l)*WA(J5+jcnt) ELSE xxxs(0)*WA(jmx) 
yyyB(0:l) * mp.bg 


Over plot BG. 

Over plot BG-BGerror. 


oplot ,xxxs ,yyys 
yyyB (0 : 1) * mp.bg - mp . bgerr 
oplot ,xxxs ,yyys ,linestyle*2 
yyys(0: 1) * mp.bg + mp .bgerT 
oplot , xxxs ,yyys ,linestyle=2 
plotlabl.WA. SPECT , J J ,mp .mcntrl ,up 
if mp. order gt 0 then PRINT,’ Order’ 
if mp. order gt 0 then PRINT, mp. order 
if kvtd then begin 

xyouts , 0 . 5 , 1 . 1 , FNAM , /NORMAL 
plotconf ig,-l , ’ ’,’ ’,-l,kdev,’’ 
endif 

IF mp . cam NE 0 THEN BEGIN 

IF mp.cam LE 2 THEN PRINT,’ long’ ELSE PRINT,’ short’ 

END ; 

mp.cntrl B l » Set to norm. 

RETURN ! Go back 

END ; 


Over plot BG+BGerror. 
Customized labels . 
Print order # . 


Put file name on plot. 
Send plot A -> term. 

Valid IUE CAM i? 


PRO SETUP, ML, RS.FF.dtl.nulld.nside.mp, mcntrl, up 

-***************************** SETUP. PRO ************************************ 


Gets or creates a new .DTL file as well as other set up information 
for DETAILS. by Charles L. Joseph June 1981 May 1990. 


ML 

RS 

dtl 

nulld 


nside 

mp 

mcntrl 

up 


«i»ti ttiiiib number of measurements that can be made. 

Red Shift (Doppler velocity) used in LUT . 

Structure that holds contents of .DTL file 
A nulled or zero version of the variable dtl 
* of points on each side of the expanded plots - set to 60 
MSLAPparameter - a structure, used to pass many arguments, 
structure (see DETAILS. PRO) 

UserParmeter - 30 pt floating array, exclusive for user. 
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I 

j ****** ******************************* ***************************************** 

— « j 


KL*200 
1-0 

ns id# * 150 
START: STAR-’ 1 


• of lines 
■at autir index 
• of expansion points, 
name loop 


print , ’ ’ 

READ, ’What is the output file name ? ’ .STAR 
IF STAR EQ ’H’ THEM BEGIN 
SPAWN,’ Is *.DTL’ 

GOTO, START 

END 

PRINT,’ ’ 
close, 1 

openr , i , STAR+ ’ . DTL ’ , ERROR « errtst 
close, 1 

epp - * » 

if errtst eq 0 then begin 
CLOSE, 1 

0PENU,1,STAR+’ .DTL’ 

ML*200 

readu,l,dtl 

sz ■ where (dtl.iaf ne 0.) 

■z * size(sz) 

if (sz(0) eq 0) then 1*0 else I * sz(l) 
close, 1 

print ,string(7B) , ’File Already Exists’ 
print, I,’ measurements have already been made 
print , ’ ’ 

print , ’Append new data to the old?’ 
read, ’(Note: a NO will erase old data.)’,app 
endif 


; need help ? 

; list old files 
; start again 


; Does file exists? 

; Set append flag to "NO" 

; File Found 

; Get the old file. 

; MAX meas. that can be made 
; Read the file. 

; Get t of measurements. 

; If some made, set counter. 
; Close file for now. 

; Flash Screen and print. 

’; Indicate # already made. 

; Append or erase A start? 


if (app ne ’Y’) AND (app ne ’y’) then begin 
close, 1 

0PENW,1,STAR+’ .DTL’ 

dtl * replicate({ noda }, 200) 

writeu, 1 ,dtl 

CLOSE, 1 


If APPEND flag NOT set in 
"file exists" section. 
Hake new file. 

Empty array for storage . 
Write it to file. 

Close until latter. 


1*0 
endif 
print , ’ ’ 

print, ’The radial velocity is UBed only to help identify species.’ 


print, ’Enter 0 if you are uncertain of the real value.’ 

READ, ’Enter Delta-Lambda/Lambda or Radial Velocity for the source: ’ ,RS 


IF ABS(RS) LT 1.0E-02 THEN BEGIN 

PRINT, ’The Delta-Lambda/Lambda is: ’ ,RS 
END ELSE BEGIN 

PRINT, ’The Radial Velocity is: ’,RS,’ km/s’ 
RS=RS/299792.5 


assume D-l/1 
tell assumption 
else rad-vel 
say so 

convert to D-l/1 


END 
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FLAGSET.FF.mp.up 


CLOSE. 1 

®p. STAR « STAR 
mcntrl. I * I 
EXIT: RETURN 


END 


get data type 


; go back to DETAILS 


PRO CONTIM . EISP , EXVA , CNT , WA . EPS , IF , XI , JE , ILIMIT ,mp .mcntrl .up . adat a . mdat a , $ 
xmlld.RS .aside 

; **************************** CONTIM.PRO *********************************** 

9 

; TO DETERMINE A CONTINUUM AND END POINTS FOR THE INTERGRATION 

; By Charles L. Joseph June, 1982 

; By Charles L. Joseph and Edward B. Jenkins May 12, 1990 


EXSP 

EXWA 

CNT 

VA 

EPS 

XF 

XI 

JE 

ILIMIT 

mp 

mcntrl 


up 

s data 
mdat a 
nulld 
RS 

aside 


expanded flux vector 

expanded wavelength vector 

continuum vector 

wavelength vector 

data quality vector 

right stop X value for integration 

left stop X value for integration 

index of the line center 

array containing start/stop to fit regions. 

MSLAPparameter - a structure, used to pass many arguments, 
structure (Bee DETAILS. PRO) 

UserParmeter - 30 pt floating array, exclusive for user. 

structure containing a single dtl(k) element 

same as sdata, but for 4 elements. 

same as sdata, but all elements are Bet to 0 

red shift (Doppler velocity) for LUT 

the number of data points on each side of the expansion. 


**************************************************************************** 


mdata(O) = sdata 

m * 0 

smx ■ mu (EXWA, jwmx, MIN=wmn) 
START: imes ■ 0 
csav * 0 

ernes * fltarr(100,2) 
kget * 0 
iopt ■ 2 

AC « dblarr(9,10) 

BC * dblarr(16) 
besav * dblarr(10, 16) 

XWRK * fltarr(602) 

YVRK * XWRK 
VSAV ■ XWRK 


Set Null Polynom. Fit 
Find max A min lambdas, 
i of descrete pts to fit. 
Do Not Overplot last CNT. 
Holds cursor readings . 
Index and # for XLIMIT. 
Cntrl option, UCURS0R.PR0 
For CURFIT.PRO 
For CURFIT.PRO, fit para. 
To save various BC arrays 
Working wavelength vector 
Working flux vector. 
Temporary holding vector. 
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• 1 else lcs * sz(l) 


CTT ■ 0.0*EXWA 
XLIMIT « fltarr(30) 

KM ■ intarr(9) 

COHFAC » 1.0 
sz ■ size( JE) 
if sz(0) eq 0 than lcs 
»tt ■ fltarr(lcs,2) 

777 • fltarr (2) 

777(1) ■ !y.crange(l) 

777(0) ■ !y.crange(0) 

7half » 0.5* (777(1) -777(0)) + 777(0) 
fdsav ■ fltarr(2,mcntrl.mtot+l) 

kdev ■ get env( ’TERM’) 

if kdev ne ’eun’ then kdev * 'xterm’ 

— Plot the spectrum, either smoothed or not 


; Make empty continuum. 

; Fitting end points. 

; For CURFIT.PRO. 

; Find the t of line cntrs . 

; Array for overplot cntrs. 
; Array for overplot cntrs. 
; Set values for top to 
; bottom over plots . 

; For cursor positioning. 


; I Window or Sunview? 


; Find min and mar flux. 

; Min k max plus max (index) 
; Adjust for plotting. 


IF mp.SMO ne 0 then PL0T,EIVA,SM00TH(EXSP,7) ,/YNOZ 
if mp.SMO eq 0 then begin 

ymn ■ min(EXSP,J5, MAX=7mx) 
xmx « max(EXWA, jmx, MIN=xmn) 
ymn ■ ymn - 0.1* (ymx-ymn) 
ymx * 1 . l*ymx 

if !x. range (0) OR !x.range(l) OR !y.range(0) OR !y.range(l) then begin 

if !x. range (0) OR ! x. range (1) then begin ; Use these values instead? 

xmn ■ !x. range (0) ; If one is nonzero «> user 

twit x !x. range (1) ■ wants to over ride, 

endif 

if ?y r ang e (0) OR !y.range(l) then begin ; Use these values instead? 

ymn * !y.range(0) 
ymx * !y.range(l) 
endif 
endif 

plot, exwa,exsp,/YROZ,xrange= [xmn, xmx] ,yrange=[ymn,ymx] ,xstyle=l ,ystyle=l 
asy * f indgen(16)*(!PI*2/16.) 
usersym,0.5*cos(asy) ,0.5*sin(asy) ,/FILL 
oplot , exwa , exsp ,psym=8 
endif 

OPLOT, WA,EPS<0 
for k*0, lcs-1 do begin 

if lcs eq 1 then xxx(k,0:l) * EXWA(JE) 
if lcs gt 1 then xxx(k,0:l) * EIWA(JE(k)) 
oplot, xxx(k, 0:1) ,yyy,linestyle=l 
endfor 

for k*0,3 do print,’ 


Show individual points. 

Plot data quality. 

If mult, profiles, then 
flag them on the plot . 

Use dotted line. 

Let user know in C0NTIM. 


print,’ »»»> START of Continuum Fitting Routine «««« 


for kk=0 , 1000 do begin 


Do a large # of attempts . 


if (ki It 1) then tvcrB, 0.5, 0.5, /NORMAL 
if (iopt eq 2) then begin 


; 1st time? Place cursors. 
; 1st time? Instructions. 



40 



print , ’ * 

print, ’Left Mouse to locate END POINTS of featureless continuum segments’ 
print,’ - Only horizonatal position will be used.' 

print , ’ - Up to 15 segments are allowed . ’ 

print,* Right Mouse button provides additional MENU options.' 
print, 'Click Center Mouse when finished.’ 


•ndif 

CHANGE: ucursor.IX,IY.opt,iopt,l,mp ; Graphical input A contrl. 

if II gt wmx then II « EIWA(jwmx-l) ; Don’t except out of range 

iy XI it no then II * wain ; Don’t except out of range 


Branching in Continuum Routine 

if ! ERR eq 2 then goto, calc ; Middle Mouse used in 

; UCURSOR.PRO. 

if !ERR eq 1 then begin ; Left Mouse used. 

XLIMIT(kget) ■ II l Store Fitting Limits, 

k g at * kg fiti + 1 | Upd&to counter ■ 

yout « 0 . 07* ( ! y . crange C 1 ) - ! y . crange (0) ) + !y.crange(0) 
xyouts.ix.yout, ’ : ’+strtrim(kget,2) ; Show the fitting limit, 

endif 


if iopt eq 2 then begin * Define discrete 

; continuum points . 
imes * 0 

plotconf ig,0, ’ ’,’ ’ ,3,kdev, ’ ’ ; Open small instr. window. 

xout3 ■ !x. crange (0) 

yout * 0.5*(!y.crange(l)-!y.crange(0)) + !y.crange(0) 
xyouts,xout3, yout, ’Use Left Mouse - Continue until Middle Mouse’ 
wset,0 ; Control back to plot. 


•ERR « 0 


while (!ERR ne 2) AND (imes It 100) do begin 
ucursor , ix , iy , opt ,kchoice , 0 ,mp 
; if kdev eq ’xterm’ then wait.l 

if !ERR ne 2 then xyouts.ix.iy, ’x’ ; Show discrete points. 

if !ERR ne 2 then ernes (imes, 0) * ix ; Update discrete param. 

if !ERR ne 2 then ernes (imes, 1) * iy 

if !ERR ne 2 then imes « imes + 1 ; Update total number, 

endvhile 


wdelete,3 

wset.O 

print , ’ ’ , string (7B) 

print, ’»» Resume finding regions 


; Remove instruct . window . 
; Make sure back to graph. 
; Bell and give instruct, 
of featureless continuum <«<<’ 


endif 


Redefine line center(s). 


if iopt eq 3 then begin ; 

print , ’ ’ 

print, string(7B),’»» Use LEFT Mouse to locate New Line Center ««’ 
print,’ NOTE: If the RIGHT Mouse iB UBed first - then up to 3’ 

print,’ additional line centers can be defined.’ 

yout ■ 0.5*(!y.crange(l)-!y.crange(0)) + !y.crange(0) 
if les gt 1 then JE * JE(0) 
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TVCRS ,EXWA ( JE) .yout , /DATA 
cursor, ix.iy 

if kdev eq ’xterm’ than wait,l 
xtst * !err 
!err * 0 

▼tap * where (exwa ge ix) 

JE » Ttmp(O) 

kk ■ 0 

if xtat eq 4 than begin 
xxx (0:1) - ix 
oplot , xxx , yyy , linaatyl a* 1 
Ttmp ■ JE 
JE * intarr(4) 

JE(0) ■ Ttmp 
for k«l,3 do bagin 

rado: w-tmp ■ strtrim(string(4-k)) 


; Tamp, atora kay stroke. 

; Gat indax of exwa cor- 
; responding to input. 

; If RT. Mouse, do more. 

• 

; Dotted line of nee canter 
; Save JE to make JE(k). 

; Make JE into an array. 

; Put 1st line center in. 

; String of remaining #. 


print , ’ ’ 

print.’Use LEFT Mouse to find up to ’.vtmp,’ more line centers, 
print, ’Use MIDDLE Mouse if finished.’ 
cursor, ix.iy 

if kdev eq ’xterm’ then wait,l 
if !err eq 2 then goto.lcfin 
vtmp * where (exwa ge ix) 

JE(k) « vtmp(0) 

COMPARE , ix ,RS , sdata ,mp .mcntrl 
if mp.cntrl eq 0 then goto, redo 
xxx(0:l) * ix 
oplot , xxx , yyy , 1 ines tyl e* 1 
mdata(k) « sdata 
kk « kk + 1 


Simple graphical input. 

; X window? Slow down. 
Finished? Get out . 

Get index of exwa cor- 
responding to input. 

; Findout which LUT entry. 
If Not a good ID. 

Good ID, get ready to 
plot it, dotted line. 

Update Multiple Meas. 

Increment counter. 


endfor 

endif ; 

lcfin: if (kk gt 0) then JE * JE(0:kk) 

goto, START 
endif 


END Redefine Line Cntrs. 

; How many line cntrs? 
Now go fit continuum. 


if iopt eq 4 then begin ; Adjust f of points in the expanded plots? 

for k=0,2 do print,’ ’ 

print , ’MOVE CURSOR into the COMMAND WINDOW to input the number of points’ 
print, ’in the Expanded Plots. MSLAP starts with 150 points on either’ 
print,’ side for 301 total. NOTE: PROGRAM RETURNS TO CALLING ROUTINE’ 
n.side ■ FlX(nside) 

print ,’ Current setting is:’,nside,’ for a total of ’ ,2*nside+l, ’ points’ 
print , ’ ’ 

read, ’Enter the number of pointB on each side ’ ,ntmp 
if ntmp gt 0 then nside * ntmp 

Aside * FIX (nside) ; Make sure an integer, 

goto, DONE 
endif 

if iopt eq 5 then goto. START ; Restart CONTIM.PRO. 

IF opt eq 9 THEN begin k STOP ft goto, CHANGE ft end ; Stop for kbd input. 

IF opt eq 6 THEN GOTO, DONE ; Get out of CONTIM.PRO. 
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IF opt eq 3 THE* BEGIN 
dtl(mcntrl.I) ■ nulld 


PRINT,' Lino ’.mcntrl.I,’ has boon DELETED! 


■cntrl . I«mcntrl .1-1 
IF mcntrl.I LT 0 THEN mcntrl.I»0 
END IF 

IF opt oq 4 than bogin 
road.COMS 

■data.com « byte(string(COHS,format-’ (alO) ’)) 
•ndif 

IF opt «q 2 THEN GOTO, DONE 
if opt lo 4 than goto, CHANGE 
if opt aq 7 than usrprogl.WA.SPECT, JJ.mp.up 
if opt aq 8 than usrprog2,WA,SPECT, JJ.mp.up 
if (opt aq 7) OR (opt oq 8) then goto, CHANGE 
if opt aq 5 then begin 

print,’ ’ ,otring(7B) ; 

print,’>>> HARD COPY ia HOT available HERE «<’ 
print,’ One can be made after each measurement ’ 
print , ’ ’ 

•ndif 

if opt aq 11 then goto, ST ART 
andfor 


; Delete last .DTL measure. 

’ .STRING (7B) ; Print notice 
; Set I to over write. 

; Can’t go negative. 

m 

9 

• Comment on .DTL measure.? 


Return to Main Prog 
Ram. .Del . .Nothing 


; Usrprog called. 

; Hardcopy requested, 
but is not allowed 
; Give Warning 


New plot parameters. 


End Branching A Selecting Continuum Constraints 


calc: print,’ ’ 

if imes ge 1 then IWRK(O) « ernes (0:imes-l,0) 
if imes ge 1 then YWRK(O) * ernes (0:imes-l,l) 

if (kget It 2) AND (imes It 2) then begin ; »> HOT ENOUGH END POINTS? «< 
k choice ■ wmenu([’Exit Continuum Fitting?’ , ’YES' , ’NO ’3 ,title*0,init*l) 
if kchoice ne 1 then goto, CHANGE ; Go get more option OR tell 

; calling routine ABORTED t 

; get out . 

Calculating various polynomial fits and the statistics to perform 
F TESTS for the validity of using higher order polynomials.’ 

— Please Wait — ’ 

Strip out those portions that are featureless and put compactly 
into vector XWRK. Start at XWRK(jjj) each time through. 


mp . cntrl * 0 
goto, OUT 
•ndif 
print , ’ 
print , ’ 
print , ’ 


k « 0 

if imes gt 0 THEN jjj=imes-l ELSE jjj*0 
for ii*0,kget/2-l do begin 
xtst - XLIMIT(2*ii) ; ing 

while EXWA(k) It xtst do k=k+l 


; Start index after discrete pts 
; Search for beginning and end- 
indices for featureless 

; portions and stuff into XWRK. 


jji * k 

xtst ■ ILIMIT(2*ii+l) 

while EXWA(k) It xtst do k=k+l 


332 - k 
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XWRK(jjj) - EXWA(jjl:jj2) 
YWRK(jjj) » EXSP(jjl:jj2) 

jjj • jjj ♦ JJ2 ’ jj* 

andfor 


; Put into working array. 

; Put into working array. 

; t of wlwBwnts in working array 


if (imes It 3) AID (kget eq 0) than bagin ; Straight lina fit? 

VSAV » whara (EXVA ga ernes (0,0) AID EXWA la cmes(l,0)) 

DYDI » (c»es(l,l)-cmeB(0,l))/(cBes(l,0)-CBe8(0,0)) 

XWRK » EXWA (VSAV) 

TWRK » cmes(0,l) ♦ DYDX*(XWRK-cme8(0,0)) ; Linear fit. 

•z * size (XWRK) 
jjj ■ sz ( 1) - 1 
andif 


xml ■ max(IWRK(0: jjj) ,min«xmn) 

VSAV « where (EXWA ge xmn and EXWA la xmx) 

J1 * VSAV(O) 
sz - size (VSAV) 

J2 - VSAV(sz(l)-l) 

XO * (zmx + xmn) / 2.0 
XWRK » IWRK(0 : j j j) - XO 
YWRK « YWRK(0:jjj) 

VSAV - YWRK 

jjj - iii ; - i 

if (jjj It 6) then begin 
print ,string(7B) 

print , ’Warning: Poorly Defined Continuum’ 
print , ’Returning to Previous Level’ 
return 
andif 

Find the best chi-squared fit 

EE ■ dblarr(16) 5 To hold residuals of var. fits 

ntot * mcntrl .mtot+2 
if mtot gt 9 then mtot = 9 

for m * 2,mtot+l do begin ; order polynomials. 

MP1 * m + 1 
M2 * 2*m - 2 
BC(*) ■ 0.0 

CURFIT,XWRK,YWRK,AC,BC,MM,m,MPl,M2, 1, jj j .2, 0.0,0 
bcsav (m-1,*) ■ BC ! Save the fit parameters. 

EE(m-l) « total(YVRK(i : j j j-3)“2) ; Sum of square of residuals. 

YWRK * VSAV ; Replace destroyed spectra. 

endfor 

ENE * float (jjj-3)/C0HF AC 5 * of independent samples. 

EE « EE/COHFAC 

mcntrl. HE « ENE \ Save for posterity. 

c.choice ■ string(bytarr(31 .mcntrl .mtot+6)) 
fdist « mcntrl. fdist 

m * 0 


Find range of abscissa. 

VSAV is temporary array. 
Determine starting and ending 
points for continuum vector 
and store in J1 and J2. 

Adjust range to center on 0 
AIbo adjust size of vectors 
to remove trailing zeros. 

Make extra copy of YWRK 

; Mot enough points to fit? 
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BStp ■ 0 

flw ■ where (ene ge fdist(0, 1 : 17) ,f ent) 
if «ne eq f diet (0, font) than begin 
fin ■ fdist(l,fcnt) 
f2n - f diet (2, f ent) 

•ndif ELSE begin 

fin ■ (ene-f diet (0 ,f ent) )/(f dist(0,f cnt+l) _ f diet (0,f ent)) 
f2n ■ fln*(fdist( 2 ,fcnt+l)-fdist( 2 ,fcnt))+fdi 8 t( 2 ,fcnt) 
fin ■ fln*(f diet (l,fcnt+l)-fdist(l,fcnt))+fdist(l t f ent) 


•ndelse 

flp ■ string(format*’(f6.2)’ ,fln) 

f2p ■ string (format*' (f6.2) * ,f2n) 

for k*l,mcntrl.mtot do begin 
fahd ■ (EME-k-3) * (EE (k) -EE (k+1) ) /EE (K+l ) 
fahd2 « (ENE-k-4) * (EE(k) -EE (k+2) ) / (2 . *EE (k+2) ) 
fdsav(O.k) * fahd 
fdsav(l,k) * fahd2 
if mstp eq 0 then m * k 

if (fahd le fin) AND (fahd2 le f2n) then mstp * 1 
fahd * string (format*’ (f6.2) ’ .fahd) 
fahd2 * string (format*’ (f6.2) ’ ,fahd2) 

c_choice(k) * ’ ’■*-strtrim(string(k) ,2) + ’ *+fahd+f lp+f ahd2+f2p 


endfor 

■istp * string (format*’ (i2) ’ ,m) 


final :BC ■ bcsav(m,0: 15) 

MP1 * m + 1 


; Again with the best poly . 


M2 * 2*m - 2 


XVRK * EXWA(J1:J2) - 10 ; 

YWRK * EXSP(J1:J2) J 

jjj ■ J2 - J1 ; 

CURF I T . X WRK , YWRK , A C , B C , MM , m+ 1 , MP 1 , M2 , 1 , j j j . 3 , 

CNT * float (YVRK(l:jjj-3)) I 

plot ,EXWA ,EXSP , /YNOZ » 

oplot , exwa , exsp ,psym=8 I 

sz * size (CNT) I 

xc - EXWA(J1 : Jl+sz(l)) \ 

oplot, xc, CNT » 

if esav eq 1 then oplot ,xc ,cntsv,linestyle*2 
vtmp * where (XLIMIT gt 0,nxs) 
for k*l,nxs do begin l 

ix - XLIMIT(k-l) 


Take only valid part of 
wavelength and Bpectra. 

Total f of points in the fit. 
0.0,0 ; Fit A make continuum. 

YWRK replaced with continuum. 
Plot original spectrum. 
Overplot with large dots . 

How many elements are there? 
Size wavelength accordingly . 
Over plot fit continuum. 

; Overplot prev. continuum. 

Show range of continuum fit. 


xyouts.ix.yout, ’ : ’+strtrim(k,2) 
endfor 


• _ — Ge t the choice of Polynomical and Type of Integration 


getch: print,’ ’ A print,’ ’ 

mstr * string (format*’ (i2) ’ ,m) 
if mstr eq mstp then begin 

print, ’The recommended Order iB’,mstp,’ (highlighted), but you may select another.’ 
endif ELSE begin 


oplot,xxx(k.O:l),yyy,lineBtyle«l ; Flag the plot too. 

TVCRS, xxx (k, 0) ,yhalf, /DATA ; Repeat Cursor command, insure, 

print , ’ * 

print, ’Cross-Hairs INPUT the left edge of the integration range.’ 
print,’ - Only horizonatal position will be used.’ 

CURSOR, IX, IY ; R«ad left (beginning) point, 

if kdev eq ’xterm’ then wait,l ; X windows too fast, slow it. 

PRINT, STRING (7B) ; Ring the bell, 

if lcs gt 1 then XI(k)«IX else XI«IX ; Store start pt as array? 
print , ’ ’ 

print. ’Cross-Hairs IHPUT the right edge of the integration range.’ 
print,’ - Only horizonatal position will be UBed.’ 

CURSOR, JX,JY ; Read right (final) point, 

if kdev eq ’xterm’ then wait.l ; X windows too fast, slow it. 

PRIHT, STRING (7B) ; Ring the bell, 

if lcs gt 1 then XF(k)*JX else XF*JX ; Store stop pt as array? 


if (xxx(k,0) It IX) OR (xxx(k,0) gt JX) then begin 
print, string (7B) 

print , ’Warning: Integration Window does not contain’ 
print,’ Center Wavelength defined before.’ 

print,’ ’ 
goto, setmg 
end if 
endfor 

andif i End kchoice < 6. 


if kchoice eq 6 then begin ; kchoice=> Use Fixed Window 

if mp . window eq 0 then begin ; Has one already been selected? 

vtmp ■ 0.0 ; If not, read window l store, 

read, ’Enter Full Width of window in km/s ’.vtmp 
mp. window * vtmp 
end if 

print, ’Will use standard window of ’ ,mp. window 

vtmp ■ mp. window/2. 9979e5 ; Convert to unitless value. 

XF*EXWA(JE)+EXWA(JE)*vtmp/2. ; Wavelength integration 

XI“EXWA(JE)-EXWA(JE)*vtmp/2. ; window around line center, 

endif 


if kchoice eq 7 then mcntrl . intopt=0 else mcntrl. intopt=l ; Use ALT.INT. 


JE « JE - J1 


; Adjust line center index. 


sz « size(CNT) 

EXWA « EIWA(Jl:Jl+sz(l)) ; 
EXSP ■ EXSP ( J 1 : J i+sz ( 1 ) ) ; 
mcntrl. ESAV * sqrt(EE(m)/(ENE - float (m))) ; 


Stip out only necessary parts 
Stip out only necessary parts 
Save error of the fit. 


DONE: PRINT,” 

mp.cntrl * 1 

IF (opt eq 2) THEN mp.cntrl ■ 0 
IF (iopt eq 4) THEN mp.cntrl * 0 


Let’s get out of here. 
Set flag to OK 
or EXIT, No Measurement 
or EXIT, No Measurement 
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print, 'REMINDER: the recommended Ordar is’.mstp 
•ndelse 

print, ’To assist in other choices: X la the observational diffaranca in the’ 
print, ’raducad chi squares divided by tha raducad chi square , shich if larger' 
print, 'than tha theoretical F(l,n) indicate* that going to tha next higher order’ 
print,' polynomial is justified statistically. Y is similar to X mxcapt it is’ 
print, 'for comparison to F(2,n), an order that is 2 higher. The F Distributions’ 
print, 'are at the 6X confidence level for: ’.fix(ene),* points.’ 

print , ’ ' 

print, 'Select tha Ordar of tha Polynomial’ 

print,' - Other Polynomials may be examined before deciding.’ 

print,’ - Selecting order ’ ,mstr , ’ implies use that polynomial.’ 

c_choice(0) ■ 'Order X F(l,n) Y F(2,n)’ 

c_choice(6) ■ 'Order: '+mBtr+’, Predetermined window’ 

c_choice(7) * 'Order: *+mstr+’, Ubs ALT_INT Routine’ 

if csav eq 0 then c_choice(8) ■ ’Toggle ON Overplot Last Cont.’ 

if csav eq 1 then c_choice(8) * 'Toggle OFF Overplot Last Cont.’ 

c_choice(9) ■ 'Restart Fitting’ 

c_choice(10)« ’Abort - RETURN to Prev. Level* 

kchoice « wmenu(c_choice,title«0, init*m) 

if kchoice It i then goto, getch ; Not valid - do it again. 


. - — Implement the choice 


if kchoice eq 10 then opt * 2 
if kchoice eq 10 then goto, DONE 
if csav eq 1 then cntsv * CNT 
if kchoice eq 8 then begin 
vtmp * csav 

if vtmp eq 1 then csav * 0 
if vtmp eq 0 then csav * 1 

if csav eq 1 then cntsv * CNT 
goto, getch 
andif 


; Abort - pretend set by USRMENU 
; then go to the bottom. 

; Toggle Overplot last continuum 


TVCRS, 0.5, 0.5, /NORMAL 
XI « min(XLIHIT) 

XF « max(XLIMIT) 

if kchoice eq 9 then goto, START 
if (kchoice ne m) and (kchoice It 6) then begin ; Different order? 


Put Cross-hairs on screen. 
Set integration range to the 
maximum allowed by spectrum, 
kchoice ■> restart 


m ■ kchoice 
goto, final 
sndif 

if kchoice It 6 then begin 
TVCRS , 0 . 5 , 0 . 5 , /NORMAL 
if lcs gt 1 then XI*fltarr(lcs) 
if lcs gt 1 then XF=fltarr(lcs) 
for k*0,lcs-l do begin 
setmg: if lcs eq 1 then xxx(k,0:l) 

if lcs gt 1 then xxx(k,0:l) 
TVCRS, 0.5, 0.5, /NORMAL 
TVCRS, xxx (k,0) ,yhalf ,/DATA 


Set order to user choice. 
Go display different order. 


else 

else 


XI' 

IF 1 


kchoice *> Settled on order I 
Put Cross hairs on screen. 

*0. ; Multiple features? If 

=0. ; so, get ready. 


EXWA(JE) ; Get feature center. 
EXWA(JE(k)) ; Do Bame if multiple ones 
Put Cursors on the graph. 

Put Cursors on the feature. 


IF opt eq 6 THEN mp.cntrl ■ 2 


; or EXIT to NEXT spect . 


OUT: print,' ' 

up .poly ■ m ; Store in MS LAP structure. 

Bcntrl.fln ■ fd»av(0,m) 

Bcntrl.f2n ■ fd*av(l,m) 

RETURN 
END ; CONTIK 


PRO STORE ,WW , ML , dtl .Bdata ,mp .mcntrl ,up 
;******************************** STORE .PRO 

9 

; TO WRITE OUT THE FILES FOR DETAILS 

; by Charles L. Joseph June 1982 May 1990 

9 

; ML maximum number of measurements that can be made. 

9 

******** ************************ ****************** **************************** 


; go back to DETAILS 


******************************** 


iafs * Bdata. iaf 
iafs * where(iafs ne 0,sz) 
if sz eq 0 then goto, RET 
IF (mcntrl. I+SZ) GE ML THEN BEGIN 
ERROR * 66 

PRINT, 'Too many lines measured!’ 
YN= ’ ’ 

READ, 'use <RETURN> to exit’.YN 
GOTO, RET 

END 

IF mcntrl. 1+10 GE ML THEN BEGIN 
PRINT, STRING (7B) 

PRINT, 'Within 10 measurements of 
WAIT, 4 


; Get the Ion Codes. 

; How much incoming data? 

; forced exit if too 
; many lines, set flag 
STRING (7B) ; tell them 


; exit 

; go back to DETAILS 
the maximum allowed’ 


END 

IF mcntrl. I LT 0 THEN mcntrl. I * 0 ; can’t have I < 0 

openu , 1 ,mp . STAR+ ’ . DTL ’ 
readu.l.dtl 
close, 1 

for k*0,sz-l do begin ; Add all new data to old 

•data * mdata(k) 
dtl (mcntrl . I ) * Bdata 

print, FORMAT®’ ("Storing measurement:" ,i3,3x,al0) ’ , $ 
mcntrl. I, string (sdata. el) 

mcntrl . I*mcntrl . 1+1 ; increment masterindex 

IF mcntrl. I GE ML THEN mcntrl . I®ML-1 ; catch over run 

endf or \ 


openu , 1 ,mp . STAR+ ’ . DTL ’ 
writeu,l,dtl 


;or write file A exit 
; store data on disk 
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close file 
or write file t exit 

go back to DETAILS 
to exit normally 


pro f store, tau,sdata,mp,mcntrl,up 

;******************************* FSTORE . PRO ******************************** 

V 

; To store the flux, wavelength, continuum, and optical-depth vectors. 

; by Charles L. Joseph May 1990 


; ************************************************************* **************** 

; Plot Optical Depths before storing them 

aax * fltarr(5) 

aay * aax 

aax* [0. ,0. ,-0.25,0. ,0.25] 
aay * [0., 3. 0,2. 6, 3. 0,2. 6] 

!y. margin (0) ■ 4. 

v * 2.9979e5*(tau(0,*)- sdata.wD/sdata.wl 
asy * findgen(16)*(!PI*2/16.) 
user8ym,0.5*cos(asy) ,0.5*sin(asy) ,/FILL 
plot,v,tau(2,*) ,xtitle*’velocity’ ,ytitle=’log tau’,psym=8, $ 
yrange=[-2.5,1.5] ,ystyle*l 
usersym,0.S*cos(asy) ,0.5*sin(asy) 
oplot,v,tau(l,*) ,psym=8 
oplot,v,tau(3,*) ,psym=8 
usersym , aax , aay 

ttmp * where(tau(2,*) gt 1.0, tent) » Which are limits? 

if tent gt 1 then oplot ,v(ttmp) ,tau(2, ttmp) ,psym=8 ; Over plot lower limits 

if tent eq 1 then oplot, [v(ttmp) ,v(ttmp)] , [tau(2, ttmp) ,tau(2, ttmp)] ,psym=8 
ttmp * where (tau (3,*) gt 1.0, tent) ; Which axe limits? 

if tent gt 1 then oplot ,v(ttmp) ,tau(3 , ttmp) ,psym*8 ; Error’s Lower Limits, 

if tent eq 1 then oplot, [v(ttmp) ,v(ttmp)] , [tau(3, ttmp) ,tau(3, ttmp)] ,psym=8 
fval * string (format* ’ (f6.4) ’ .sdata.f ) 
fval * ’ f * ’+fval 
print , ’ ’ 

print, 'Log Tau Plot Rest Wavelength: ’ , sdata.wl , ’ ' ,f val 

kwtd * wmenu([’Copy?’ , ’YeB’ , ’No’] ,title=0, init=2) 
if kwtd then begin 

plotconf ig, 1, ’ ’,-l.kdev,” 1 Config. for Hardcopy, 

asy * findgen(16)*(!PI*2/16.) 
usersym,0.8*cos(asy) ,0.8*sin(asy) ,/FILL 

plot,v,tau(2,*) , xtitle* ’velocity * ,ytitle=*log tau’ ,psym=8, $ 
yrange= [-2 .5,1.5] ,ystyle*l 


flue userBymuui vi 
arrow for L.L. 
Values for an arrow 
(Lower Limit) . 
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user sym, 0.5* cos (asy) ,0.5*sin(asy) 
oplot, v,tau(l,*).psym«8 
oplot ,▼ ,tau(3 ,*) ,psym*8 
usersym , iu , aay 

ttmp - where(tau(2,*) gt 1.0,tcnt) ; Which are limits? 

if tent gt 1 then oplot ,v(ttmp) ,tau(2, ttmp) ,psym*8 ; Over plot lower limits 
if tent oq 1 then oplot, [v(ttmp) ,v(ttmp)] , [tau(2,ttmp) ,tau(2, ttmp)] ,psym*8 
ttmp » vhere(tau(3,*) gt 1.0, tent) • Which axe limits? 

if tent gt 1 then oplot ,v(ttmp) ,tau(3, ttmp), psym*8 ; Error’s Lower Limits, 

if tent eq 1 then oplot, [v(ttmp) ,v(ttmp)] , [tau(3, ttmp) ,tau(3,ttmp)) ,psym*B 
xxo»0.0*(!x.crange(l)-!x.crange(0)) + !x. crange (0) 
yyo«1.06*(!y.crange(l)-!y.crange(0)) ♦ !y .crange(O) 


xy out s , xxo , yyo ,mp . FMAM 

xxo«0.4*(!x.crange(l)-!x.crange(0)) ♦ !x. crange (0) 
xyout8,xxo,yyo,8tringCsdata.el)+strtrim(string(sdata.wl) ,2) 

xxo«0 . 7* ( ! x . crange ( 1 ) - ! x . crange (0) ) + 'x.crange(O) 
xyouts ,xxo ,yyo .mcntrl . date 

plotconf ig,-l , ’ ’ ,-l,kdev, ’ ’ i Send plot k -> term, 

endif 


First Retrieve k Update Header Information 


as ■ assoc(5,fltarr(200)) 
tas ■ aB8oc(5,fltarr(2,600)) 
close, 5 

openr,5,mp.STAR+’ .T1U’ , ERROR = errtst 
close, 5 

if errtst ne 0 then begin 
openw , 5 ,mp . STAR+ * . TAU ’ 
for k=0,5 do as(k) * fltarr(200) 
close, 5 
endif 

openu,5,mp.STAR+* .TAU’ 

iafs * as(0) 

rwav ■ as(l) 

bgs * as (2) 

bgers ■ as (3) 

dates * as (4) 

lwfs * as (5) 


; Association variables for 
; random disk access. 

; Test if file exists. 

; Ho file so create one. 

; Open a new file. 

; Hake 6 dummy records . 

; How close it. 

; Open file for Update k get 
; 1st record - species codes 
; Get prev. rest wavelengths 
; Get previous BG levels. 

; Get previous BG errors . 

; Get previous date codes . 

; Get prev. log f -lambdas. 


itau = TOTAL (iafs ne 0) 


; Total number before? 


iafs (itau) * long(sdata.iaf) 

print , iafs (itau) 

if iafs (itau) It 0 then stop 

rwav (itau) * sdata.wl 

bgs (itau) * mp.bg 

bgers (itau) ■ mp.bgeir 

dates (itau) * 910112 

lwfs(itau) * aloglO(sdata.wl*sdata.f ) 


; Add new species code. 


Add current lab wavelength 
Add current BG level. 

Add current BG error . 

Add current date. 

Add new log f -lambda. 


Hext Store, Header, Optical Depths + Spectra 
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tautmp ■ fltarr(2,600) 
sz ■ size(tau) 
if sz(O) 1* 0 then begin 
print, string (7B) , ’ ’ 

print,’ WARNING: No Optical Depth Data" 


wait, 2 
close, 5 
return 


andif 
•z ■ sz(2) 


; Make temporary array. 

; Find size of incoming data 
; If no incoming , get out . 

; Ring the alarm BELL. 

; 2 seconds to see it. 

; Close the file. 

; Vo saving it, return. 

; Change info, to a scalar. 


as(0) ■ iafs 
as(l) ■ rwav 
as (2) ■ bgs 
as (3) ■ bgers 
as (4) ■ dates 
as (5) * lwfs 

tautmp(0:l,0:sz-l) * tau(0:l,*) 
tas(3*itau+l) * tautmp 

tautmp(0: l,0:sz-l) * tau(2:3,*) 
tas(3*itau+2) * tautmp 

tautmpCO: 1 ,0:sz-l) * tau(4:5,*) 
tas(3*itau+3) « tautmp 


Store updates of record- 
keeping and header info, 
back to the disk file. 


Stuff Wavelength It TAUs(-) 
and write to disk. 

Stuff TAU’s and TAUs(+) 
and write to disk. 

Stuff Spectrum t Continuum 
and write to disk. 


close ,5 > Close file and return. 

return 

end 


PRO EXPAND, JO, JS,WA,WS,JT,JE, J, ns ide 

;************************* EXPAND. PRO *********************************** 

9 

; TO EXPAND THE FLUX AND WAVELENGTH VECTORS. 

; by Charles L. Joseph Sept 1980. 

f 

; JO flux vector 

; JS expanded flux vector 

; WA wavelength vector 

; WS expanded wavelength vector 

• jt number of points in expanded vectors 

; JE central index of vectors ( expanded ) 

; J index of cursor position 


************************************************************************** 


nsave * ns ide 
NPTS=T0TAL(WA GT WA(0)) 
JT ■ 2*nside + 1 


; Save the request. 

; # of points in vector 
; # of expansion points 
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if (HPTS It JT) than aside « FIX(irPTS/2) - 1 ; 

if J It aside than jfact * aside - J else jfact * 0 ; 
if (NPTS-J) It aside than jfact » IPTS - J - aside ; 
JS ■ JO(J*jfact-aside: J+jf act+nside) 

VS * WA(J+jf act-aside: J+jfact+aside) 

JE » aside ♦ 1 - jfact ; 


aside ■ asave 


Adequate # of points? 
Too dose to 1ft edge? 
Too close to rt edge? 


Index of line center. 
Restore the request. 


RETURN 

END 


fMtft##fttttttt#t##f#t#tttt«t##tf#ttftttftttt#ttf #i«tff*f fititffttt* 

MAIN PROGRAM 

by Charles L. Joseph May 12, 1990 


date*’Fri Apr 13 00:00:01 1990’ 
date * systime(O) 

date * strmid(date,4,7)+strmid(date,20,4) 
up * fltarr(30) 

ST: for k«l,13 do print,’ * 

print , date 
HUM - 0 

libr ■ ’/u/clj/mslapdir ’ 
spawn.’head -26 ’+libr+’/main.mslap’ 
for k=l,3 do print,’ * 

READ, ’Which option would you like ? ’ ,NUM 
IF NUM EQ 10 THEN GOTD.OUT 


Get time and date 
Strip out date 


reset NUM 

Defines the MSLAP lib. 
print 1st 31 lines. 

what to do 
to get out 


DETAILS OR POSTO 


IF NUM EQ 1 THEN DETAILS, 0, date, libr, up 
IF NUM EQ 2 THEN DETAILS , 1 .date ,libr, up 
IF (NUM ge 3) AND (NUM le 6) THEN BEGIN 


run DETAILS no tau’s 
run DETAILS with tau’s 
POSTO, EDATDTL , COGS 


I 

GSTAR: STAR® ’ ’ ; set up loop 

READ, ’What is the file name, created in Option #1 or f2 ?’,STAR 

IF STAR EQ ’H’ THEN BEGIN ; Want help? 

SPAWN, ’Is *.DTL’ ; Give a listing and 

GOTO, GSTAR : go back try again. 

END ; 

if NUM ne 5 then openr , 1 , STAR+ ’ . DTL ’ , ERROR * errtst ; See if file exist, 
if NUM eq 5 then openr , 1 , STAR+ ’. TAU ’ , ERROR * errtst ; See if file exist, 
close, 1 ; Close again, 

if errtst ne 0 then begin ; Oops no file 

YN0=” ; 
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if IUM eq 4 than read, 'Create Dummy File?’,YNO ; EDATDTL? Use dummy? 
if (YNO eq 'Y') OR (YNO eq 'y’) then goto, CHOICES 
ygO*’' J Try another fils? 

READ, 'Would you like to try another FILE ?’,YN0 
if (YIO eq 'Y') OR (YHO eq ’y’) then goto,GSTAR 


GOTO, ST 

END 

END 

CHOICES : 

IF HUM EQ 3 THEM POSTO.STAR 
IF IUM EQ 4 THEM edetdtl.STAR 

• 

IF IUM EQ 5 THEN montau, STAR, date 

f 

IF IUM EQ 6 THEN cog, STAR, date ,libr 


else goto start 


run POSTO 
run EDATDTL 


GOTO, ST 

OUT: ; 

END 


go back to start 
; tty mode 


j *************************************** ************************ ************* 

; This file contains many procedures or indirect callB to compile procedures 
; used by M.S.L.A.P. These contain some of the standard as well as the 

; dummy routines that may be over layed by editing the mslap.pro routine 

; in the user’s local directory. 

; Charles L. Joseph December 24, 1990 

j **************************************************************************** 

•/u/cl j /mslapdir/mantau 

•/u/ cl j /mslapdir/ intgrt . pro 

C/u/ cl j /mslapdir/ edatdtl 

•/u/cl j /mslapdir/posto 

•/u/cl j /mslapdir/cog 

•/u/cl j /mslapdir /plotconfig 

•/u/cl j /mslapdir/compare .pro 

pro dgetl,FNAME,w,spect ,q,ihdr,ID,mp,up 

mp.dget * 0 ! These are dummy routines 

ID * ’NOT Being Used’ ; which reserve slots for 

return i read ones that will be 

and : plugged in. 

pro dget2,FNAME,w,spect ,q,ihdr,ID,mp,up 
mp.dget * 0 
ID * ’NOT Being Used’ 
return 
end 

pro dget3,FNAME,w,spect,q,ihdr,ID,mp.up 
mp.dget * 0 
ID * ’NOT Being Used’ 
return 
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end 

pro dget4,FNAME,w,»pect ,q,ihdr,ID,mp,up 
■p.dget ■ 0 
ID ■ ’NOT Being Ueed’ 
return 
•nd 

pro dg®t 5 , FNAME . w , epect , q . ihdr ,ID,mp,up 
mp . dg«t * 0 
ID « *I0T Being Used’ 
return 
end 


***************************************************************************** 
PROGRAM CURFIT -calculates ploynomial fits for CONTIM. 

By Edward B. Jenkins 

Translated to IDL and Edited by Charles L. Joseph May 1990 

ORDER OF POLYNOM. » M ; MP1 * M + 1 AND M2 - 2*M - 2 
ARRAY WILL BE PROCESSED FROM INDEX "NS" TO INDEX "N" 

IF SOME POSITIONS IN THE ARRAY BETWEEN INDICES NS AND N ARE TO BE IGNORED, 
JUST SET THE APPROPRIATE Y VALUES GREATER THAN 1.D30 
MODE * 0 JUST EVALUATES ANSWER VECTOR B 
MODE * 1 EVALUATES B AND REPLACES Y BY BEST FIT CURVE 
MODE *= 2 EVALUATES B AND REPLACES Y BY RESIDUALS 
MODE = 3 COMPUTES BEST FIT CURVE FROM EXISTING X AND B ARRAYS 
IMUL * O.DO TELLS SUBROUTINE TO REFERENCE X ARRAY; OTHERWISE X IS IGNORED 
LOGIND « 0 DOES ORDINARY FIT; IF LOGIND .GT. 0, COEFF’S ARE IN TERMS OF LOG Y 
IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION YCN) ,A(M,MP1) ,B(M2) ,MM(M) ,X(N) 

********** ************************************************************ ******* 


PRO CURFIT, X.Y, A ,B, MM, M.MP1, M2, NS, N, MODE, XMUL, LOGIND 
IF MODE EQ 3 then GOTO, COMP 
A * 0.0*A 


; Go compute polynom. 
; Zero A and B arrays 


B * 0.0*B 

A(0 ,0) * N - NS + 1 ; A(0 ,0)= # of points 

FOR I « NS,N DO BEGIN 

IF Y(I-l) GE l.e70 THEN A(0,0) « A(0,0) - 1. ELSE BEGIN 
IF IMUL EQ 0 THEN XP * X(I-l) ELSE XP = double (I )*XMUL 


II * IP 


B(0) ■ B(0) ♦ XP 

IF LOGIND GT 0 THEN Y(I-l) * double(alog(Y(I-l))) 
A(0,HP1-1) * A(0,MP1-1) + Y(I-l) 

FOR J * 2, M DO BEGIN 

A(J-l.MPl-l) = A(J-l.MPl-l) + Y(I-1)*XP 
IP * IP* XX 
B(J-l) « B(J-l) + XP 
ENDFOR 

IF LOGIND GT 0 THEN Y(I-l) * double (EXP (Y (I- 1) ) ) 
FOR J= MP1.M2 DO BEGIN 


IP - IP*IX 

B(J-l) « B(J-l) + IP 
ENDFOR 
ENDELSE 

ENDFOR I ® CONTINUE 

MS - 2 

FOR I-l.M DO BEGIN 

FOR J - MS, N DO BEGIN 
XPJ ■ X * J - 2 
1(1-1, J-l) - B(IPJ-l) 

ENDFOR 
MS » 1 
ENDFOR 

EQSOL ,1,M,MP1 ,B,MM ; CiLL EQSOL 

IF MODE EQ 0 THEN RETURN 
COMP: FOR I - IS.N-1 DO BEGIN 

IF IMUL EQ 0 THEN IP « 1(1-1) ELSE IP * double (I-1)*XMUL 
II « IP 
YCOMP ■ B(O) 

FOR J « 2,M DO BEGIN 

YCOMP « YCOMP + B(J-1)*XP 

IP - IP*XX 

ENDFOR 

IF LOGIND GT 0 THEN YCOMP = double (EXP (YCOMP)) 

IF (MODE EQ 2) AND (Y(I-l) LT 1.D70) then Y(I-l) * Y(I-l) - YCOMP 
IF MODE NE 2 THEN Y(I-l) « YCOMP 
ENDFOR 
RETURN 
END 


************************* ****************************************** ******* 

PROGRAM EQSOL - To solve Differential Equation for CURFIT. 

By Eduard B. Jenkins 

Translated to IDL and Edited by Charles L. Joseph May 1990 
************************************************************************** 


pro eqsol,a,n,nn,x,in 
; a « fltarr(n.nn) 

; x ■ fltarr(n) 

; ■ ■ intarr(n) 
for i * 1 ,n do begin 
m(i-l) - 1 
am ax « a(i-l ,0) 
for j * 2,n do begin 

if ((abs(a(i-l, j-l))-abs(amax)) gt 0) then begin 
amax * a(i-l , j-l) 
n(i-l) * j 
end if 
andfor 

for j » l,nn do begin 
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if(amax eq 0) than print, 'about to divide by amax*0 in EQSOL’ 
a(i-l, j-1) » a(i-l,j-l)/amax 
andf or 

for ip » l,n do begin 
if (ip - i) na 0 than bagin 
■mm ■ mCi-1) 
moult ■ aCip-1 

for j"l ,nn do if (j-«mm) aq 0 than a(ip-l,j-l) ■ 0 $ 

alaa a(ip-l.j-l) « a(ip-l , j-l)-zmult*a(i-l, j-1) 
and if 
andf or 
andf or 

for i « l,n do begin 
■mm » m(i-l) 
i (mnm-1) * a(i-l,nn-l) 
andf or 


return 

and 

; PROGRAM UCURSOR to provide continuous readout of cutbot position unitl 
; the mouse is clicked. The left and middle mouse buttons return to calling 
; routine, while the right button brings up a menu of options. If the 
; graphics device is not a SUN, control information is can be passed by the 
; key used to input the cursor location. 

I 

; cflg * 0 use Main UserMENU. Selection is returned in opt parameter. 

; cflg = 1 use Continuum MENU first. Selection is returned in iopt param. 

; cflg < 0 BO Menu, reserve right mouse button for other purposes. 

; Coordinates are returned in IX and IY . 

; See main.mslap for definitions of the mp structure. 

; By Charles L. Joseph 6/4/90 Latest Revision: 1/7/91 


pro ucuxsor,IX,IY,opt,iopt,cflg,mp 
▼fig * 0 
START: ; 


opt * 9999 
iopt * 9999 
kdev * getenv( ’TERM ’ ) 
if kdev eq ’vt300’ then kdev = ’xterm’ 
if (kdev eq ’xterm’) OR (kdev eq ’sun’) then 
xoutl * 0.0*(!x.crange(l)-!x.crange(0)) + 
xout2 * 0.5*(!x.crange(l)-!x.crange(0)) + 
youta * 0.5*(!y.crange(l)-!y.crange(0)) + 
tst ■ 0 

while (tst ne 1) do begin 
cursor, IX, IY, /CHANGE 


begin 

•x.crange(O) 

!x.crange(0) 

!y .crange(O) 

; Continuous readout 
; coordintes. 


of cursor 


56 


if !ERR eq 
if !ERR eq 


then tst 
than tst 


if !ERK eq 4 than tat * 1 

vast , 1 


•rut 


if true - left mouse 

if true - middle mouse 

if true - right mouse 

Set to cursor printing window 

Clear and print new coords. 

IX 


if wflg then xcrd*2.9979e5*(IX-wavrest)/wavrest ELSE xerd 
xyouta ,xout 1 ,youta , string (format* ’ (f 9 . 3) ' ,xcrd) 
if (youta It 600) then xyouts,xout2, youta, stringCfoxmat*’ (f9.3) * ,IY) 
if (youta ge 600) then xyouts ,xout2, youta, string (format*’ (f9.0) ’ ,IY) 
vaet,0 • Return to plotting window. 

•ndwhile 

if kdev eq ’xterm’ then wait,l 

if !ERR eq 2 then return ; Exit this routine if middle 

if ! ERR eq 1 then return • or left mouse used, otherwise 

a MENU is wanted. 


OLD Tektronics type Cursors 


•ndif else begin 

CURSOR, IX. IY : Hhat to do? 

PRINT, STRING (7B) 5 4025 

ondelse 


Get Various UserHenu Options 

if cflg It 0 then return i N° Menu Options Wanted 

; under any circumstances. 

If cursor reading is out side of plot area, or the right mouse button used, 
then set menuflg=l to get a menu. If calling routine is the Continuum 
fitting routine (cflg-1), then call special menu defined below. Otherwise 
(cflg =0) , call Main UserHenu. 


menuflg * 0 

if((II It !x.crange(0)) or (IX gt !x .crange(l))) then menuflg ■ 1 
if ( (IY It !y.crange(0)) or (IY gt !y .crange(l))) then menuflg * 1 
if (!ERR eq 4) and (cflg eq 0) then usrmenu.opt ; Not in Continuum Routine, 

if (!ERR eq 4) and (cflg eq 1) then menuflg * 1 ; In Continuum Routine. 


if menuflg then begin 

c_opts * string(bytarr(22,7)) 
c_opts(0) * ’Continuum Options’ 
c_opts ( 1) * ’Take No Action’ 
c_opts(2) * 'Descrete Cont. Points’ 
c_opts(3) * ’Define Profile Cntr(s)’ 
c_opts(4) * ’Change EXPANSION size’ 
c_opts(5) * ’Restart Continuum Fit’ 
c_opts(6) * ’Call MAIN userMENU’ 
iopt * wmenu(c_opts , title=0, init=l) 
if iopt eq 6 then usrmenu.opt 
end if 


; First MENU for Continuum 
; routine. Option 4 calls 
; Main UserMenu. 


Get iopt option. 


57 


if (opt eq 10) OR (opt eq 12) then begin 
print, string (7B),’ * 
if opt eq 10 than begin 

print, 'Find Wavelength of Interest and Click any Mouse ’ 
print, 'I-Coordinte will than read in velocity’ 


Set up to display velocity 


endif 

if opt eq 12 then begin 

print , 'Locate Backround Level and Click LEFT Mouse ■then’ 
print, 'offset vertically by 1 standard deviation and Click again.’ 
print, 'If you wish to enter the numbers manually. Click RIGHT Mouse.’ 


endif 
tst ■ 0 

while (tst ne 1) do begin 
cursor , IX , IY , /CHARGE 
if !ERR eq 1 then tst * 1 

if !ERR eq 2 then tst ■ 1 

if !ERR eq 4 then tst * 1 

wset , 1 
erase 

xyouts.xoutl ,youta, string (format® ’ (f 9. 3) ’ ,IX) 
xyouts ,xout2,yout a, string (format*’ (f9.3) ’ ,IY) 


Continuously read wave- 
length position of cursor. 

Cursor reading provides 
reference wavelength. 


WBet.O 

endwhile 

wavrest « IX • Start this routine over, 

. i ; using velocity instead of 

if opt eq 10 then goto, START > wavelength, 

if opt eq 12 then begin 
if !ERR eq 4 then begin 

print,’ ’,string(7B) 

Print, 'Move cursor into COMMAND window for input.’ 

Read , ’ Ent er background : ’ , bg 
mp.bg ■ bg 

read, ’Enter uncertainty : ’ ,bg 
mp.bgerr * bg 
endif ELSE begin 
print, ’ ’ 

print, ’Now locate a 1 sigma deviation.’ 
mp.bg * IY 
cursor , IX, IY 

if IY gt mp.bg THEN mp .bgerr=IY-mp .bg ELSE mp.bgerr=mp.bg-IY 
endelse 

endif 

endif 

if opt eq 11 then begin 
print,’ ’,string(7B) 

Print, 'Move cursor into COMMAND window for input of New Plot Parameters.' 
print,’ ’ 

read, 'Enter min and max values for the X axis: ’,xi,x2 
read, 'Enter min and max values for the Y axis: ’,yl,y2 
!x. range ■ [ xl, x2 ] 

!y. range * [ yl, y2 ] 
endif 


RETURN 

END ; UCURSOR 


PRO dataget ,WA,SPECT, JJ,IHDR,FF,mp,mcntrl,up 

f 

;***************************** dataget .pro a***************************** 

9 

; To gat the Spectrum Data. — formerly inside GRAPHS. PRO 

; This routine interfaces the data-getting routines of DGETS to DETAILS 

; By Charles L. Joseph May 1990 


SPECT 

flux vector 

WA 

wavelength vector 

IHDER 

header vector 

FF 

data type flag 

JJ 

data quality vector 


************************************************************************** 
COMMON BLOCK l.DTY.SMO ; get common 


; 1st time through 
; create blank name 


FNAME « mp.FNAM 
IF mp.cntrl EQ -1 THEN begin 
START: FNAME * ’ » 

for k*0,l do print,’ ’ 

print , ’Enter complete INPUT Data Filename’ ; 
read, ’Including the path if necessary ’.FNAME; read file name 
if (FNAME eq ’end’) OR (FNAME eq ’END’) then begin 
mp.cntrl * -99 
return 
andif 

if FNAME eq ’H’ then begin 
for k=0,2 do print,’ ’ 
spawn, ’Is’ 

for k=0,2 do print,’ ’ 

print, ’»>>> Note: entering the word "end" exits <<«’ 
goto, ST ART 
endif 
close, 1 

openr, 1 ,FNAME, ERROR = errtst 
close, 1 

if errtst ne 0 then begin 
print,’ ’, string (7B) 
print , ’Warning: Unable to read file’ 

yno ■ ’ ’ 

read, ’Try to input again? ’ ,yno 
if (yno eq *1’) OR (yno eq ’y’) then goto, ST ART 
mp.cntrl * -99 
return 

endif 


; close logical unit 
; check that it exists 

; a problem if nonzero 
send beep 
print warning 


; return to MSLAP 
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COHFIRD: cohfac ■ 0 

print, string (7B) , ’ * 
print , * ’ 

print, ’IMPORT AIT: Tl»e coherence length is used to calculate ALL uncertainties . ’ 
print, ’It ia the number of pixels influencing the value in a given pixel. 1 
print, ’Thus, the coherence length must be 1 or greater. For example, data 
print ,’ smoothed by a 3-point Running Box Car has a coherence length of 3.’ 
print .string (7B) , * ’ 

read, ’What is the coherence length of the present data?’ .cohfac 

■p. cohfac * cohfac 

if mp. cohfac It 1 then goto, COHFIND 

up. bg « 0.0 ; Assume background t 

ap.bgerr « 0.0 ; error initially, 

andif 


I 

mp.FIAM * FRAME 

IF FF EQ 1 THEN dgetl , FRAME, VA .SPECT, JJ .IHDR.ID.mp.up 
IF FF EQ 2 THEN dget2, FNAME.WA, SPECT, JJ, IHDR.ID.mp.up 
IF FF EQ 3 THEN dget3, FNAME.WA, SPECT, JJ, IHDR.ID.mp.up 
IF FF EQ 4 THEN dget4, FNAME.WA, SPECT, JJ, IHDR.ID.mp.up 
IF FF EQ 5 THEN dget5, FNAME.WA .SPECT, JJ, IHDR.ID.mp.up 
if (FF It 1) OR (FF gt 5) then print , ’Invalid File Option’ 
print , ’ ’ 

print, 'Rev Data Have Been Read — Data Getting Option: ’ ,fix(FF) 


if mp.cntrl eq -1 then begin 
sz * size(WA) 

if sz(0) eq 0 then mp.cntrl * -99 
if mp.cntrl eq -99 then begin 

print .string (7B) , ’WARNING: Data-Getting 
print,’ that has ZERO elements! 

vait ,3 
end if 

if sz(0) eq 0 then return 
wavetol * 10*(WA(sz(l)-l)-WA(0))/sz(l) 
if wavetol It 0.5 then wavetol * 0.5 
mcntrl.wtol « wavetol 
endif 


Determine tollerance for 
finding entry in ULUT. 

Just in case problem with DGET, 


Routine returned a vector’ 
— Returning to MENU’ 


indicate an abort A exit . 
Tollerance * +/- 10 pixels, but 
never less than 0.5 Angstrom. 
Save for COMPARE. 


print,’ ', string (7B) 

print , ’ 

| The background uncertainty frequently has a major impact on the | 

| uncertainty of the various measurements. I 

| Current background is: ’.mp.bg,’ with an error: ’ .mp.bgerr, ’ I’ 

| Use the RIGHT Mouse Button if these are unsatisfactory. I 

| Rote: all graphical input is performed by placing the cross-hairs | 
| at the point of interest and pressing a mouse button. I 


print 

print 

print 

print 

print 

print 


print, 


return 

end 

; •*#**#*«##****#*****#** END of MASTER. AUX Package ########################## 
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PRO COMPARE. VV.RS.edata,mp,mcntrl 
*************************** COMP ARE. PRO 


********************************* 


TO PHD INTERSTELLAR LIVES IV VARIOUS LOOK UP TABLES 

VW wavelength of feature to be identified 

RS red shift (Doppler Velocity Features) 

sdata data structure for single measurement 

■p MSLAPparameter structure 

Charles L. Joseph 3/28/79 Latest Rev.: 4/5/90 

*************************************************************************** 

WF»WW/(RS+1.) : •PP 1 ? r#d *hift 

nullid ■ { notdat, el: bytarr(15), iaf: long(O), wl: 0.0, f: 0.0 > 


uidl * nullid 

aids ■ replicate ({ notdat }, 2000) 
uids * replicate ({ notdat }, 200) 
wrk * replicate ({ notdat }, 10) 
vrk2 * replicate ({ notdat )•, 10) 
wrk2 . el « -32000. 


Make a "single" entry. 
Make MSLAP LUT array. 

Do the Bame for user’s. 
Make small working 
structures . 

Unused ones to odd val . 


close, 3 

if (mp.DTY eq 1) OR (mp.DTY eq 3) then $ 

openr,3,mcntrl .libr+’/tabdata/ilut .tab’ ; Open LUT file, 
if mp.DTY eq 2 then openr,3,mcntrl . libr+’/tabdata/mlut.tab’ 
if mp.DTY eq 4 then openr.S.mcntrl.libr+’/tabdata/qlut.tab’ 


Open LUT file. 
Open LUT file. 


readu,3,mids 
CLOSE, 3 

OPENR , 3 , ’ ulut . t ab ’ , ERROR = errtst 
if errtst eq 0 then readu,3,uids 
vtst * where (mids .wl le WF, N) 
if V le 2 then V « 2 
vrk2(0:4) ■ mids(H-2:H+2) 


Read LUT structure. 

Get User LookUpTable 
if one exists. 

Which is best entry? 

Into working structure. 


if mp.DTY eq 3 then begin 
close, 3 

openr ,3 ,mcntrl . libr+ ’ /tabdata/mlut . tab ’ 
readu,3,mids 

vtst * where (mids .wl le WF, N) 
wrk2(5:9) ■ mids(N-2:N+2) 
endif 


; If two tables are to be 
; used, replace strange 
; values with real ones . 

; Which is best entry? 

; Stuff into working str. 


W - wrk2 . wl 

for k*0,9 do begin ; Order combined small 

vtst ■ min(W.ibest) ; structure according to 

wrk(k) * wrk2(ibest) ; wavelength. 

W(ibest) « 32000. 
endfor 


W = wrk.wl 
W « abs(W-WF) 


; Find absolute diffs. 
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rtst « min(W,ibest) 
if ibest It 1 than ibest*l 
OK: print,’ ’ 


; ibest is array index. 


Set up to shoe choices found in the tables 


nuids ■ 0 

USERIDS , WW ,HC ,WF .nuids ,n_names ,ns ,uidl ,uids ,ncntrl ; Any User ID’s? 

ajMM t(O) • ’Which Entry?’ 

n_names(l) ■ ’lone - Return to Spectra’ 

; found in user Lookup. 

for k « 0,2 do begin 

uidi * wrk( ibest- 1+k) ! Get next entry. 

DD « (WW-uidl . wl) /uidi . wl*299792 . 5 ; Convert to velocity. 

AI0K - string (uidi. el) + ’ ’ 

AI0H * strmid(AI0N,O,15) 

sk • string(format*’ (f 9.3) ’, uidi .wl) 

fk ■ string (format*’ (f 7.4) ’, uidi .f) 

n.names (k+2) * wk+’ ’+string(f ormat* ’ (f6. 1) ’ ,DD)+* ’+AI0N+’ f*’+fk 

endf or 

n.names (nuids+5) * ’Input Identity’ ! nuids * f of lines 

print, ’Select Look-Up Table Entry — Mote: Hone is an option.’ 

print,’ — 2nd number is the Doppler velocity (km/s).’ 

WHICH ■ wmenu (n.names , title=0,init=3) 

NNN * WHICH - 3 
mp.cntrl * 1 

IF mp.cntrl EQ 0 THEN NNN=0 ; 

if (WHICH le 1) then mp.cntrl * 0 

Three Bources of input 

if (WHICH EQ (nuids+5)) then begin ; ================ ; Input new ULUT entry. 

HNN * 999 

print, ’>>>> Requesting ID <<<<’ 

USERIDS, WW ,HC ,WF , HNN, n.names ,ns ,uidl ,uids .mcntrl 

sdata.el * uidi. el ; Put codes into struct- 

sdata.iaf - uidl.iaf I for MSLAP measurements, 

sdata.wl * uidl.wl 
sdata.f * uidl.f 
WL*uidl .wl 

AION - string (uidi .el (0: 1)) 

MUL * 0 
HNN * 999 
endif 

IF ((HNN GE 2) AND (HNN le 100)) THEN BEGIN ; =======; Get ID from user table 

H - ns (WHICH) 

print, ’Using ULUT entry #:*,N 
uidi * uids(N) 

sdata.el * uidi. el 

sdata.iaf * uidl.iaf 

sdata.f * uidl.f 
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■data. si ■ uidl. si 
1I0N » atring(uidl.el(0:l)) 
KUL-0 


CLOSE, 3 

EVD 

if (WV It 2) than begin ; ■« 
ibaat ■ ibaat ♦ VBICB - 3 
if ibaat It 0 than ibaat ■ 0 
nidi ■ wrk(ibeat) 

•data. el » uidl.al 
sdata.iaf * uidl.iaf 
•data. el » uidl.al 
adata.f * uidl.f 


EVD 

if adata.iaf It 0 then print , ’adata. iaf : 
if sdata.iaf It 0 then atop 
CLOSE, 2 


1 .sdata.iaf 


RETURN 

END 


; Gat from stndrd LUT. 

; Gat requeated index. 

; Make sura range ia OK. 

; Gat ionization atate. 

« 

• Close look up table. 

; Go back to DETAILS. 


* ********************************************************************** 
i 

; PROCEDURE USERIDS to allow the user to identify lines 

I 

; by Charles L Joseph 13-Aug-84 

9 

• ********************************************************************* 
f 

PRO USERIDS ,W ,NC ,VF .nuids ,n_names ,ns ,uidl ,uids .mcntrl 

ncnt * whereCuids .wl ne 0., NC) i get # of lines 

W ■ uids.wl 

if (nuids gt 100) then goto, IDIT ; Input identity? 

; This part searches for existing entries 

i.nanes * string (bytarr (26 ,26)) • Big enough to hold max 

pa ■ intarr(26) » V of Menu selections, 

for K*0,99 do begin 

IF ABS(WOO-WF) LT mcntrl. wtol THEN BEGIN ; Close to that observed? 

uidl * uids(K) 

DB® (WW-W (K) ) /W (K) *299792 . 5 ; Convert to velocity. 

AI0N « string (uidl. el) + * * ; Start building string 

AION * atrmid(AION,0, 15) ; for the menu selection, 

ek * string^ormat®’ (f9.3) * ,W(K)) 

i_names(nuids+5) * sk+* ’+string (format® ’ (f6. 1) ' ,DB) 
fk * atring(format= > (f7 .4) * .uidl .f ) 
i_names (nuids+5) * i_names(nuids+5)+’ ’ +AI0N+’ f=’+fk 

ns(nuids+5) * K > Store ULUT addresB. 

nuids ® nuids + l ; Increment the counter . 
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if (nuids gt 20) than goto, OUT 

END 


endfor 

OUT: n_namee ■ i_names(0:nuids+5) 

close ,3 

RETURN 


; Can only hold 20 entries 

; Reduce Menu Selections 
; to the correct size. 


lev entries to the User LookUpTable (LUT) are made here 


ID IT: print,’ ’ 

print, ’Enter Complete Species Name - include ionization and fine structure.’ 
print, ’Do not leave any preceding blanks. Enter for example: "Fe II*" for' 
print, ’iron singly ionized and excited Fine Structure. Tliis entry will only’ 
print, 'be used as an element identifier and print label. The ionization will’ 
print, ’have to be entered later.’ 
vdum ■ ’ ’ 
read , vdum 
elm«vdum 

uidl.el * 0*uidl.el 

elm ■ strtrim(elm,l) + ’ ’ 

uidl.el ■ byte(strmid(elm,0,15)) 
vdum * ’ ’ 

READ, ’What is the IONIZATION state [ 2 = II 4 * IV ] ?’,vdum 

IA*fix(vdum) 

vdum ■ ” 

if IA le 0 then begin I Encode as molecule or 

print , ’ ’ 

print , ’Species is assumed to be a molecule’ 
read, ’What is the vibrational level?’, vdum 
iaf * 10*fix(vdum) 
vdum * 9 9 

read, ’What 1 b the rotation level? (9 Max) ’, vdum 
iaf * iaf + fix (vdum) 
vdum = st raid (elm, 0,2) 

if (vdum eq ’H2’) OR (vdum eq ’HD’) then begin 
read, ’Enter 1000 if this is part of the Werner system ’,ly 
if ly eq 1000 then iaf « iaf + 1000 
endif 

endif ELSE begin • else encode as ion. 

READ, ’What is the FINE STRUCTURE state [ 3 « *** ] ?’,vdum 

FS*fix(vdum) 

iaf * FIX(IA*10 + FS) ; Ionization t Fine Str . 

sndelse 

READ, ’What is the Rest WAVELENGTH of the line ?’,vdum 

uidl.wl - float (vdum) J Rest Wavelength. 

READ, ’What is the OSCILLATOR strength of the line ?’,vdum 

uidl.f * float (vdum) I Oscillator Strength. 

tst * uidl.el (0:1) ; Encode as Periodic Table A Molecules 

els « ’H HeLiBeB C N 0 F NeNaMgAlSiP S ClArK CaScTiV CrMnFeCoNiCuZn’ 

els * els+ ’ H2H2HDHDC0CHCNC20HU U U’ ; Elements + Molecules 

els * byte(els) 
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3 - o 

•tat * 0 

while atat aq 0 do bagin 

•tat • ((tat(O) aq ala(j)) AND (tat(l) aq ala(j+l))) 

3 - 3*2 

if 3 g* 78 than goto.axitat 
•ndsh.il a 

•xitat: uidl.iaf - long(iaf) ♦ long (1000) *(3/2) 


; Include ionization. 


CLOSE, 2 


; close look up table 


READ.’Entar alanent permanently in table?*, ID 
IF (ID EQ *y *) OR (ID eq *Y’) THEN begin 
uida(NC) « uidl 
close, 3 

openr , 3 , * ulut . tab * , ERROR « errtat ; See if file exists and 

close, 3 I make one if necessary, 

if errtst eq 0 then openu,3, ’ulut. tab’ ELSE opens, 3, ’ulut .tab’ 
sriteu,3,uids i Update ULUT. 

•ndif 
CLOSE, 3 
RETURN 
END 


PRO POSTO .STAR 

; ****************************** POSTO. PRO ********************************* 

V 

; To list the measured lines from DETAILS and create ASCII files of 

; this data. 

; by Charles L. Joseph December 19, 1980 

; Latest Revision: January 7, 1991 

; STAR file name 

» 

****************************************************** ********************* 
LPF«0 I flag t o pause for LP 

STRT: close, 1 

OPENU , 1 , ST1R+ ’ .DTL' ; Open the .DTL file, 

nulld * { noda, el: bytarr(15), iaf: longCO), wl: 0., f: 0., owl: 0., eqw: 0., $ 
me: 0., fm: 0., fme: 0., am: 0., ame: 0., com: bytarr(lO), up: fltarr(30) > 


dtl * replicate ({ noda }, 200) 
edata * nulld 
readu,l,dtl 
close , 1 

sz * where (dtl . iaf ne 0.) 
sz * size(sz) 

if (sz(0) eq 0) then NTL = 0 else NTL = sz(l) 

IF NTL EQ 0 THEN GOTO, EXIT 
LOOP : VTD*0 

Bp awn, ’clear’ 

for k-0,3 do print,' 9 

print,' *■*«**■■«*■»*■* Program POSTO. PRO =*=«=* 


Make array for storage 
Make working copy 
Read the data. 

Get f of measurements 

t store in NTL. 
return if none 
Bet up loop 


I I 


print, 

print,' It is best to resize this window to the full screen while in POSTO.’ 
print, 'This can be done now end a reminder to unzoom the window will appear’ 
print, 'at the appropriate time.’ 
for k*0,4 do print,’ ’ 

PRINT, ’ Options to MANIPULATE TABLED DATA:’ 


PRINT,’ 

1: 

Reorder .DTL file by Laboratory Wavelength’ 

PRINT,’ 

2: 

Reorder .DTL file by Observed Wavelength’ 

PRINT,’ 

3: 

Reorder .DTL file by Ion (Use Opt 1 or 2 first) ’ 

PRINT,’ 

4: 

List existing .DTL file’ 

PRINT,’ 

5: 

Create Customized Table in ASCII format’ 

PRINT,’ 


Note: all .DTL files should be reordered’ 

PRINT,’ 


in the same fashion first.’ 

PRINT,’ 

10: 

EXIT - Back to MSLAP ’ 

print , ’ 

1 


READ,’ 

Which do you want ?’,WTD ; which one 

IF WTD EQ 10 

THEN GOTO, EXIT 

LPF-0 


; lots of time for LP 


IF WTD EQ 3 THEN PRINT, ’NOTE: It may be necesary to use Option 1 or 2, first’ 
IF (WTD GE 1) AND (WTD LE 3) THEN SRT , dtl , WTD , STAR 
IF WTD EQ 4 THEN BEGIN 

PUBL, sdata, dtl .STAR, NTL 

YN0=’ » ; want a listing? 

READ, ’Press <ENTER> when ready to proceed’ ,YN0 

END 

if WTD eq 5 then begin 

mkcusttab, STAR, dtl ; Make Custom Table, 

goto, STRT ; Go re-read orig. data, 

endif 

if WTD It 6 then goto, LOOP 

EXIT: print, string (7B) ; Flash screen. 

print , ’Unzoom the this window IF it has been expanded to full screen.’ 
yno * ’ ’ 

read, ’Press the <RETURN> key when ready to continue. ’, yno 
RETURN 

END ; 


S 

PRO SRT, dtl, WTD, STAR 

;**************************** SORT. PRO ********************************* 

f 

; TO SORT .DTL DATA FILES by Charles L. Joseph 22-Jan-1980 

j dtl the primary structure holding the measurements of DETAIL . 

; WTD What-To-Do Flag 1 *> sort by Wavelength 

; 3 => sort by Periodic table then molecules. 

; STAR name of the file 

9 

;***** ************************************************************ ******** 

if WTD eq 1 then begin ; Sort by Lab. Wavelength 

nslot ■ dtl.wl 
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nslot ■ 10000. - nslot 
•ndif 

if WTD #q 3 than nslot ■ dtl.iaf 

n*lot ■ (nslot > 0) ♦ 99999.* (nslot eq 0) 

tdtl ■ dtl 


kk * 0 

for k » 0,199 do begin 
▼y ■ nin(nslot.ky) 
if tj It 99999. than begin 
tdtl(kk) « dtl(ky) 
nslot (ky) • 99999. 

kk * kk + 1 


andif alee goto, out 
andfor 

out: dtl » tdtl 

openu , 1 , STAR* ’ .DTL’ 
sriteu.i.dtl 
close, 1 
RETURN 
END 


; Ordering is in NSLOT 

; Sort by Ion «> NSLOT 
; Flag to ignore null data. 
; Make copy of data. 

; Perform sort from min to 
; max. 


; Make changes permanent by 
; writing them to file. 

; back to POSTO 


I 

PRO PUBL.sdata, dtl, STAR, NTL 

; ******************************* PUBL .PRO ******************************* 

i 

; To print the results for POSTO. 

; By Charles L. Joseph November 1980 

; sdata a structure holding the measurements of a single profile. 

. dtl a structure bolding the primary data >> many sdata’ s. 

; STAR name of the file 

* NTL number of measured lines 

• ******* ********************************************************* ********** 


PHDR® ’ Ion Lab f 

N*-l 

PG=0 

NT=Q 

for k*0,4 do print,’ ’ 

PAGE: if Id.name eq ’TEK’ then erase 
PG*PG+1 

starp * STAR*’ ’ 

starp * strmid(starp,0,15) 
pgp « string (format®’ (i2) ’ ,PG) 
PRINT, starp, ’ 

PRINT,’ ’ 

PRINT, PHDR 

PRINT, FORMAT® ’ (14x.a4,5x, "Value 
PRINT,’ ’ 

NT-NT+44 

IF NT GT NTL-1 THEN NT«NTL-1 


D-lam EQW error comments’ 
; counter for pages 
; page # 

; total counter 

; start loop 4010 
; increment page # 


Page : ’ ,pgp 

; page header 
mA".5x," Eqv Width")’ 

; and 2nd line 
; # on page 

; how many lines ? 
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LOOP: 1«H+1 

dt » dtl(N) 

AIOH » string (dt. si) 

COM ■ string Cdt.com) 

IF I CT ITL-2 THE* *T-45*PG 


; printout data 


; for pagss 


PLMN 

PRINT, FORMAT-’ (A16.1x,F7.2,lx.F7.3,lx,I5,ix,f8.3,lx.f8.3,2x.A10) ’ , $ 
AI OH , dt . si , dt . f , dt . osl , dt . sqw , dt .me , COM 


IF IT LT 45*PG THE* GOTO, LOOP 
FOR 1-0,11 DO PRINT,’ ’ 

IF I LT ITL-1 THEM GOTO, PAGE 
CLOSE, 3 
RETURN 
END 


don’t page yet 
move out of the 
set up nev page 
close .LST up 
back to POSTO 


say 


**************************************************************************** 
MKCUSTTAB - to make a customized table (ASCII format) . 

by Charles L. Joseph November 29. 1990 

******* ***************************************************** **************** 


pro mkcusttab.STAR.dtl 
frstim * 1 
sdots *’ ’ 

for k*l,6 do Bdots-sdots*’ : I* 

sdots * ’ I ’+sdots 

dash - ’ 

uline * ' ’ 

for k=l ,5 do uline * uline+’ ’ 

mstr * string(bytarr(120,308)) 

; each el: 120 char. 

miaf ■ long(intarr(301))+long(99999) 
vtmp * share (dtl . iaf gt O.tmax) 
nstrt * 0 

if tmax ge 6 then nstp * 5 ELSE nstp * tmax - 1 
ITOT * 36 
L - 0 

copt * intarr(20) 

uopt * intaxr(20) 

upc * string CbytarrC 10, 20)) 

prvopt * ’ ’ 

. * 

spe * ’ 

mstr(306) * spe+STAR 

mstr(304) * ’Species Wavelength f 

mstr(303) * 

for k*0,tmax-l do begin 


; First time flag. 

; Make an axis marker. 

; Index starts at 0. 

; Make long underline. 

; 308 el. string array 

; To hold the ID codes . 

; Table size -> tmax. 

; Starting line-display. 
; Stopping line-display. 
; t needed - Species ID. 
; Index for copt . 

; Column options, a .DTL 
; For options of UP. 

; Holds format strings . 

; Holds previous opts . 

» 

; Last few lines of mstr 
’ ; hold header info. 


sdata * dtl(k) 

spe * string(sdata.el) ; Build a string of 

spe * spe+8tring(f ormat* ’ (f 9.2) ’, sdata. wl) ; ID’s and assoc, info 

spe * spe+’ ’+string(f ormat= * (f 9 .4) * , sdata. f )+’ ’ 

mstr(k) * spe • But into mstr. 

miaf (k) - sdata. iaf ! Save ID codes, 

endfor 
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; Then print instructs. 


wrk • •»»»»»»»»»»»» CUSTOMIZED TABLE WORK SPACE * 

jjJ; * |jk+ 1 ’ 

tbl • '»»»»»»»»»»»» Partial Display of the TABLE * 
tbl » tbl+'«««««« <<<<<<<<<<< ’ 
cnt » ’ * 

for k«l,6 do cat ■ cnt+’ 12345678901234567890’ 
cat • ’O’+cnt 

MORE: ; » L ♦ 1 » * ncr * op** 0118 iadex 

if (tnAi-1 gt nstp-nstrt) AKD (nstp-nstrt It 4) then begin 
nstp ■ nstrt ♦ (tmax - 1) 
if tmnx gt 5 then natp ■ nstrt + 5 
andif 

spawn, ’dear’ • E 1188 * 

for k«0,2 do print,’ ’ » 

print, tbl : print table haader * 

print,* » ; print part of table. 

for k*0,3 do print ,mstr(306-k) ; 

for k«nstrt,nstp do print ,mstr(k) I 

print ,s dots 
print, cnt 

for k-0,3 do print,’ ’ ’• P rint instructs 

print, wrk # ^ 

print, ’First three col umn s include species name, lab. wavelength, and ’, $ 

'oscillator strength.’ 

■TOT ■ FIX(NTOT) 

print, ’Total number of characters used is:’,NT0T,’ out of 120 max.’ 
print, ’Total number of lines in the table: ’ ,tmax 
print , ’ * 

print, ’Options for Adding/Subtracting to/from table are:’ 
print , ’ ’ 

print,’ 1) Get new .DTL file’ 

print,’ 2) Delete character columns 

print,’ 3) Add Observed Wavelengths of current .DTL file’ 

print,’ 4) Add Equivalent widths of current .DTL file’ 

print,’ 5) Add Equivalent width ERRORS of current .DTL file’ 
print,’ 6) Add lBt moments (DOPPLER VELOCITIES) of current .DTL file 

print,’ 7) Add 1st moment ERRORS of current .DTL file’ 

print,’ 8) Add 2nd moments of current .DTL file’ 

print,’ 9) Add 2nd moment ERRORS of current .DTL file’ 

print,’ 10) Add comment field (Requires 10 characters)’ 

print,’ 11) Add One of the UPs (UserParameters) ’ 

print,’ 12) Finish - Create file’ 


print,’ ’ 

READ, ’Enter Option Humber’, iopt 
if iopt gt 2 then begin 
prvopt ■ prvopt+’ ’+string(f ix(iopt)) 
copt(L) * iopt 
wup * 0 

if iopt eq 11 then begin 

print, ’Which element of the UP is to be added?’ 
read, ’(Enter -1 for none)’, wup 
if wup It 0 then goto, MORE 


; Input the option #. 
; If iopt > 2 then: 
keep a record A 
store options. 

To use a UP option. 

Oops, none wanted. 


If not the same, 
find out which 
ones . 


, f . \ ; Save choice. 

uopt(L) » wup • 

_ . . , . ; Input format, 

format in: print,* ’ 

print, ’Supported Formate are: f, e, A i having a max field of 10. 
print , ’Examples : (f8.3), (el0.3), (flO.O), (i6) ’ 

wptc ■ ’ ’ 

read, ’Enter format: ’ ,wptc 
print , * 9 

tat - atrpoa(wptc,’>>) • Parenthesis used? 

If tat eq ”1 then print, 'Parenthesis must be included’ .string (7B) 
if tst eq -1 then goto.formatin 
upc(L) ■ wptc 
endif 

if iopt ne 11 then L ■ L + 1 * Incr. later if UP. 

tfaopt EQ 1) than b.gin i co1 * for 

print ,,triag(7B) , * Did f il« bad options: ■ .prvopt : Shot previous ones. 

yno “ ’ ’ 

read, ’Add all of these columns in new .DTL data?’, yno 
if (yno ne ’Y’) AND (yno ne ’y’) then begin 
tst * 0 
kst ■ 0 

icnt ■ 0 __ 

print, ’Sequentially enter each option for next .DTL file’ 

print, ’Enter -1 when finished’ 

while tst ge 0 do begin J Get options desired 

read, ’Enter next option’, tst 
if tst gt 2 then copt(icnt) = tst 
if tst gt 2 then icnt * icnt + 1 
wup « 0 

if iopt eq 11 then begin '• To use a UP option, 
read, ’Which element of the UP is to be added?’, wup 

uopt(icnt) * wup • Save choice - 

print, ’Formats are: f, e, A i having a max field of 10.’ 
print , ’Examples : (f8.3), (el0.3), (flO.O), (i6) ’ 

wptc * ’ ’ 

read, 'Enter format : ’ .wptc 
upc(icnt) * wptc 
endif 

if icnt gt 0 then L * icnt 
endwhile 

if icnt gt 0 then begin 
coptCicnt : 19) * 0 
uopt(icnt : 19) * 0 
upc(icnt:19) * ” 
prvopt * ’ ’ 

for k*0,icnt-l do prvopt=prvopt+string(copt(k)) 

endif 

if icnt eq 0 then goto. MORE • choices - abort, 

ondif 

vtmp * mstr(306) 
bsz * strlenCvtmp) 

if (ssz gt HTOT+4) then vtmp = Btrmid(vtmp ,0,NTOT+3)+’__’ 


Hew master count. 


Catch any stray 
left overs from 
from before. 
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mstr(306) ■ rtmp 

appnd.dtl.mstr.miaf ,copt,uopt,upc,HTOT,tmax 
endif 


if iopt eq 2 than begin 

print , ’Which character columns ara to be deleted?’ 
raad, ’Enter starting and ending columns’ , cl, c2 


; Gat new .DTL info. 

; Delate some columns 
; of characters. 


for k*0 ,307 do begin 

tst ■ mstr(k) 

tst ■ strmid(t8t,0,cl-l)+strmid(tst,c2,120-c2) 


mstr(k) » tst 
endfor 

■TOT » FIXCirrOT - (c2-cl) - 1) ; Adjust counter. 


endif 

if (iopt GE 3) AID (iopt LE 11) then begin 
if iopt eq 11 then begin 
wup * uopt(L) 
upf • upc(L) 

L « L + 1 


; Add an entry. 
UP option. 


Still Seeds incr. 


endif 

for k«0,199 do begin 
sdata * dtl(k) 

next string , iopt+2 ,nstrg , st itl , strsz , sdata , sup , upf 
mstr(k) * mstr(k)+’ ’+nstrg 
endfor 


ssz * 0 

if frstim then begin 
ssz * strlen(STAR) 

if (ssz gt 8trsz+2) then ssz=strsz+2 
frstim * 0 
endif 

mstr(306) ■ mstr (306) +strmid(uline,0,strsz+2-ssz) 
mstr(304) * mstr(304)+stitl 

mstr(303) ■ mstr (303) +8trmid (dash, O.strsz+2) 

■TOT * STOT + strsz + 2 
endif 

if (iopt It 12) then goto, MORE 


If first time, then 
adj . character size 
accordingly . 

No longer first. 

Adj. to file name 
* new column names 
+ more dashed line . 
Adjust counter. 


Outfile: Ytmp ■ ’ ’ 

read, ’Enter file name of the table: (type: "none" for no file) ’ ,vtmp 
if (vtmp ne ’none’) AND (vtmp ne ’NONE’) then begin 
close, 1 

openr.l.vtmp, ERROR * errtst 


close, 1 

if errtst eq 0 then begin 


yno * ’ ’ 

read, 'File already exists. Try another? ’ ,yno 
if (yno eq ’ Y ’ ) OR (yno eq ’y’) then goto, Outfile 
endif 

opens, 1, vtmp 
wrk * mstr (306) 

mstr(306) ■ strmid(srk,0,NT0T) 

for k«0,3 do printf , 1 ,mstr(306-k) ; Top of table, 

for k*0,tmax-l do printf , 1 ,mstr(k) \ Table entries. 
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printf , 1 ,mstr (303) 
close, 1 
•ndif 


; 2nd dash line. 


RETURN 

END ; mkcusttab 


¥*** 


APPXD is used by MKCUSTTAB to append additional data files to the 
customized table. 

by Charles L. Joseph Dec. 1, 1990 


; dtl structure holding the contents of the . DTL file 

; mstr an array of character strings holding the master data 

; miaf master iaf used to slot the appended file 

; jcnt the counter pointing to entry in the dtl structure. 

; for example, 3 *> 

. ************************************************************************** 
pro appnd , dtl ,mstr ,miaf , copt ,uopt ,upc ,NT0T , tmax 


STRT : ndtl * ” 
nnew * 0 
nold * 0 
read, ’Enter new 
close, 1 


.DTL file name (leave off 


openr, l.ndtl*’ .DTL* , ERROR * erxtst 


; Next .DTL filename. 
; No-New Counter. 

; No-Old Counter. 

.DTL extension) : ’ ,ndtl 

; Does file exist? 


dose.l 

if errtst ne 0 then begin 
yno * ’ ’ 

read, ’Error opening file -- try another?’ , yno 
if (yno eq ’Y’) OR (yno eq ’y’) then goto, STRT 
return 
•ndif 

openu, 1 ,ndtl+ ’ .DTL' 
readu.l.dtl 
close, 1 

niaf * dtl. iaf 
vtmp » where (niaf ne O.ktst) 
niaf * 99999* (niaf eq 0) + niaf 
dtl. iaf ■ niaf 

•dots * ’ 

•dots ■ sdots+ ’ 

•dots ■ 

dash * 

vtmp * where (copt gt 0,L) 
k « 0 
kk - 0 


If file NOT pres. 

Try again? 

Yes, get another or 
else return. 

Open the .DTL file. 
Read the data. 

Close the .DTL file. 
Get species codes. 
How many — > ktst. 
Set all 0’s -> 99999 

Make a dummy filler 
string. 

To continue dash. 
Find L*# of options. 
Index for new data. 
Index for old data. 


Hext : if k gt 199 then return 
•data * dtl(k) 
iafnxt * sdata.iaf 


; Get next line of new 
; Next species code. 
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Ttanp ■ 0 
tmpstr * * ’ 
tmptit ■ ’ ’ 

for lcnt*0 ,L-1 do begin 
jcnt ■ copt(lcnt)+2 
vup • uopt (lent) 
upf * upc(lcnt) 

next string , j ent ,nstrg . st itl , etrsz , sdat a . wup ,upf 
tjipatr » tmpstr+’ ’+nstrg 
taiptit » tmptit+stitl 
Ttnp » Ttttp + etrez + 2 
endfor 

netrg * tmpstr 
etitl « tmptit 
etrez * rtnp 
Loop: iaftet • miaf(kk) 

if iafnxt eq iaftet then begin 
mstr(kk) » mstr(kk)+nstrg 
kk « kk + 1 
k ■ k + 1 

if (k It ktst) OR (kk It tmax) then goto, Bert 
endif 

if iafnxt gt iaftet then begin 
nnew * nnev + 1 

nstr(kk) ■ mstr(kk)+’ '+etrmid(edots ,0,8trez-2) 
kk • kk + 1 
goto. Loop 
endif 

if (iafnxt ne 0) AND (iafnxt It iaftet) then begin 
nold ■ nold + 1 

for n*kk,299 do niaf (300+kk-n) * miaf (299+kk-n) 

for n*kk,299 do mstr (300+kk-n) = mstr (299+kk-n) 

niaf(kk) ■ 8data.iaf 

netr(kk) « ’ ’ 

epe ■ etring(edata.el) 

epe ■ spe+string (format* ’ (f9.2) ’ ,sdata.wl) 
epe * epe+* , +string(format*’ (f9.4) ’ ,sdata.f)+’ 
epe * epe+8trmid(edot8 , O.NTOT-36) 
netr(kk) * epe 
tmax * tmax + 1 
goto, Loop 
endif 


Get all entriee for 
next epecieB. 

In case UP get info- 
nation t format. 
Entry — > string . 
Build up all entries 


Put accumulation 
variables into 
permanent onee . 
Get next line of old 
If species agree: 

- add on 

- increment each 
counter A 

- Any more? 

; Missing line of new. 

; Fill entry with dots 
; Set index to get 
; next line of old. 

; New entry not in old 

; Shift to make room. 

; Shift to make room. 

; Add species code. 

; In case new << old. 

; Build a string of 
; ID’s and assoc. 

; information pluB 
; fill in dots. 

; Add new entry. 

; Table has 1 more. 

; Should * next time. 


Fin: uline * ’ ’ 

for k*l,5 do uline=uline+ ’ ’ 

ndtl * ndtl+uline 

mstr (306) * mstr(306)+strmid(ndtl,0,strsz) 

mstr (304) ■ mstr (304) +st itl 

mstr(303) * mstr (303) +strmid (dash, O.strsz) 

NTOT * ITOT + strsz 

print,’# of blank lines in OLD: ’.nold.’ in NEW: ’.nnev 


• 2s added to strsz 
; before. 

• Incr . column count . 
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m*«****M***«***«**M***M*****************M*»»****«M»MMff««*M 

■EXTSTRIIG it uitd to got the size ud the conversion to string of the 
next set of entries to be put in the table. Called by APPWD or MKCUSTTAB. 
by Charles L. Joseph lovember 1990 


; jcnt the index of the requested portion of the sdata structure. 

• (see the Manu al for the sdata or dtl structures.) 

; nstrg the returned Hext String 

; stitl the title of the entry (e.g. EQW, ERROR, 1st MOM., etc.) 

strsz string size, the number of columns required for the next entry, 
sdata a structure holding the measurements for a single profile. 

; wup Which User Parameter requested < — only if jcnt * 13 

; upf Holds the format if a UserParameter is requested. 

pro next st ring , j cnt , nstrg, stitl , strsz , Bdata.wup ,upf 


nstrg * * ’ 
stitl ■ * ’ 
strsz * 0 

if jcnt eq 5 then nstrg ■ string(format=’ (f8.2) ’ .sdata. owl) 

if jcnt eq 5 then stitl * ’Obs. Wave.’ 

if jcnt eq 6 then nstrg * string (format®’ (f8.4) ’ .sdata. equ) 

if jcnt eq 6 then stitl * ’ EQW (A) ’ 

if jcnt eq 7 then nstrg ■ string (format** (f8.4) ’ .Bdata.me ) 

if jcnt eq 7 then stitl * ’ ME (A) ’ 

if jcnt eq 8 then nstrg * string(format=’ (f8.2) ’ , sdata. fm ) 

if jcnt eq 8 then stitl * ’ 1st MOM.' 

if jcnt eq 9 then nstrg * string(format=* (f 8.2) ’ .sdata. fme) 

if jcnt eq 9 then stitl * ’ ERROR ’ 

if jcnt eq 10 then nstrg * string(format=’ (f8.3) ’ .sdata. sm ) 

if jcnt eq 10 then stitl * ’ 2nd MOM.’ 

if jcnt eq 11 then nstrg ■ string (format*’ (f8.3) ’ .sdata. sme) 

if jcnt eq 11 then stitl * ’ ERROR ’ 

if jcnt eq 12 then nstrg * string(sdata. com)+’ 

if jcnt eq 12 then nstrg * strmid (nstrg, 0,10) 

if jcnt eq 12 then stitl * ’ Comments ’ 

if jcnt eq 13 then begin 

val « strpos(upf ,’f ’) 
if val eq -1 then val = strpos(upf , ’F’) 

if val eq -1 then val * strpos (upf , * i ’ ) 

if val eq -1 then val * strpos (upf ,’ I ’ ) 

if val eq -1 then val * strpos (upf ,’ e ’ ) 

if val eq -1 then val * strpos (upf , ’E’) 

if val eq -1 then return 
strsz ■ strpos (upf ,’.* ) 

if strsz eq -1 then strsz = Btrpos(upf , ’) ’) 
strsz * f ix(strmid(upf ,val+l ,strsz-l)) 
val * sdata. up 


rom ur o-l a aj • 
’loating format? 
Toating format? 
integer format? 
integer format? 
Ixponential format? 
ilxponential format? 
’an't find format, 
-ind end of format, 
[f no dot. 

String size. 
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t» 1 ■ val(wup) 

nstrg * string (format *upf ,val)+ ’ • 

if strsz It 6 than nstrg - » ’♦nstrg 
nstrg ■ strmid (nstrg ,0,10) 

stitl • ’ < UP’+string(format«’(i2)’,wup)+’ > 

sndif 

if jcnt ga 12 than strsz « 10 ELSE straz » 8 


; Get array element . 
; Use format in upf ■ 
; Pad with 4 blanks? 
; Make 10 characters 
; Make column title. 

; f of characters. 


RETURN 

END 

; ••«»**********•«*«***« END of P0ST0 PACKAGE «#######*#######»#«############# 


; **************************** EDATDTL .PRO *********************************** 

9 

; To edit the .DTL files created by the Modular Spectral Line Analysis Program 

9 

; by Charles L. Joseph 6/12/84 

; Latest Revision: 8/10/90 

;****************************************************************************** 
pro edatdtl.star 


nulld * { noda, el: bytarr(15), iaf : long(0) , si: 0., f: 0., osl: 0., e<^w. 0., $ 
me: 0., fm: 0., fme: 0., sm: 0., sme: 0., com: bytarr(lO), up: fltarr(30) > 
dtl * replicate({ noda >, 200) J Make array for storage 

•data - nulld > Make working copy, 

close ,1 

OPENR, 1, STAR-*-’ .DTL’ , ERROR * errtst ; Check if file exists, 

dose, 1 

if errtst ne 0 then begin I I* no fil®. mak® one. 

close, 1 

opens , 1 , STAR+ ’ . DTL ’ 
writeu.l.dtl 
close, 1 
endif 

openu , 1 , STAR* ’ . DTL ’ i Open fil® for updating, 

readu.l.dtl * ^®®d data file. 


plotconf ig,0, ’ ' , ’ ’ ,-2,kdev, * * ; Set Hardware defaults. 

; kdev ■ getenv( , TERM’) 

if kdev eq ’xterm’ then set.plot.’X’ 

window, 0,color*2, title*’ ’ ,xpos=760,ypos=300.xsize=380,ysize=300 
!y .margin(O) * 2 J 11118 window for print - 

plot ,f indgen(lO) » infi Reminders, 

if kdev ne ’xterm’ then xO = _ 0.25*(!x.crange(l)-!x.crange(0)) + Ix.crange(O) 
if kdev eq ’xterm’ then xO = 0.0*(!x.crange(l)-!x.crange(0)) + fx.crange(O) 
yo * 0.5*(!y.crange(l)-!y.crange(0)) 
yO * fy.crange(O) 
erase 

for k*0,2 do print,’ ’ 

print,’ »» EDITING the .DTL File ««’ 
for k*0,2 do print,’ ’ 


; Send Reminders to 
; extra window . 
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xyouts ,xO ,2 
xyouts, xO,l 
xyouts ,xO,l 
xyouts ,x0 , 1 
xyouts ,x0, 1 
xyouts ,x0, 1 
xyouts, xO.O 
xyouts.xO.O 
xyouts.xO.O 
xyouts.xO.O 
xyouts.xO.O 


.0*yo+y0,’ 

,8*yo*y0, ’<CR> 
.6*yo+y0, ’ 
,4*yo+y0, ’ 
.2*yo+y0, ’ 
.0*yo+y0, ’ 

. 8*yo+y0 , * 
.6*yo+y0, * 
.4*yo+y0, ’ 
.2*yo+y0, ’ 
.0*yo+y0, ’ 


START: ML* 200 

A ■ share (dtl.wl ns 0.0, NTL) 
if (FTL le 0) than ITL » 1 

■1 • » IAF Code:’ A »2 • 

•4 * » Obs:’ A *5 « 

•7 » » 1st Mom:’ A s8 « 


Editing Reminders — ’ 

or n - goto NEXT antry ’ 

u - go UP to prev. antry’ 
c - change currant entry’ 
i - insert new data’ 
d - delate currant antry* 
t - goto TOP of data’ 
b - goto BOTTOM of data’ 
or q - lill/Quit* 

s - Write (Update) File* 
a - Exit and write changes’ 


t 

b 

k or q 


’ Rest Wave: 
•EQW: ’ 

’ 2nd Mom : ’ 


; Maximum f of measures. 


; lull data set ? 

A s3 * ’ f:’ 

t s6 * ’ Error : * 


TOP: I • -1 . 

LOOP: f*N+l • Eact time through, 

dt « dtl(N) ; P rint new line of data ' 

AION * string(dt.el) 

print .FORMAT* ’ (al5,all ,i6,al2,f9.3,a4,f6.4,a6,f9.3) ’ , AION.bI, $ 
dt.iaf ,s2,dt.wl,s3,dt.f ,s4,dt.ovl 

print, FORMAT*’ (15x,3(al2,f9. 4)) ’ , s5,dt.eqw,s7 ,dt.fm,s8,dt.sm 
print, FORMAT*’ (15x,3(al2,f9. 4)) ’ , s6,dt.me,s6,dt.fme,B6,dt.Bme 
print , ’ Comment : ’ , string (dt . com) 

print , ’ ’ 

print , ’UeerParameters ’ 
print, dt. up 
A*’ ’ 

; Read keyboard option 

A*BYTE(A) 

A*A(0) 

if (A gt 90) then A « A - 32 ! To handle lower case. 


IF ((A LE 13) OR (A EQ 78)) THEN GOTO. LOOP 
IF (A EQ 67) THEN BEGIN 
sdata * dtl(N) » 

getsdata , sdata ,N , A 
dtl(N) * sdata 

1 * 1-1 • 

END 

IF (A EQ 69) OR (A EQ 87) THEN GOTO, DONE 
IF A EQ 83 THEN STOP 
IF (A EQ 72) OR (A EQ 83) THEN N=N-1 
IF A EQ 85 THEN BEGIN 
M=N-2 

IF N LT -1 THEN BEGIN 
N*-l 

PRINT, ’TOP OF FILE’ .STRING (7B) 
END 


<CR> N Next line 
C - correction 
Set defaults 

Get Updates 
Insert in dtl 
Adj . Pointer 

I E W Exit/Write 
I S Stop 

U Up one line, 
Adj . Pointer 
Top of File? 
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END 

IF A EQ 84 THEN GOTO, TOP 
IF A EQ 66 THEN i-NTL-2 
IF (A EQ 76) OR (A EQ 81) THEN GOTO, KILL 
IF (A EQ 68) than begin 
dtl(N) - dtl(N+l :ML-1) 
dtl(ML-l) » nulld 

*TL - NTL - 1 ; 

•ndif 

i 1 (A EQ 73) THEN BEGIN i 

dtl(N+l) • dtl(N:HL-2) 
getsdata , sdata , N , A 

dtl(N) ■ adata : 

l » I - l : 

itl ■ ftl + i : 


i 


T Goto Top of file 

B Goto Bottom 

K Q Kill/Quit 
D Delete Line 


Idj . lumber tot . 

Insert Line. 

Make room in dtl 
Get new line. 
Insert in dtl . 

Adj. Pointer 
Adj . lumber tot . 


END 

GOTO, LOOP 
DONE: close, 1 

0PENU.1.STAR+’ .DTL* 
writeu, 1 ,dtl 

IF A EQ 87 THEN GOTO, TOP 
KILL: ; 

EXIT: CLOSE, 1 

•delete, 0 


; Write changes to .DTL 

; Just a write "W" go on. 

; Abort - Make no changes. 
; Close .DTL file. 

; Remove extra window. 


END 


pro getsdata.sdata.N.wtd 

; *************************** getsdata. pro ********************************** 

\ Gets new single-line of data (sdata) for the .DTL file from the keyboard. 

; If "wtd" eq 67 <=> make a correction. Then previous data is taken as the 

; default values. 

; By Charles L. Joseph 8/10/90 

; Latest Revision: 8/10/90 

.**♦*********************************************************♦***************** 


if wtd eq 67 then begin I C - corrections 

print , string (7B) 

print, ’Existing data will be taken aB the Default Values in the’ 
print, ’following questions.’ 

•ndif 
print , ’ ’ 

print, ’Enter Complete Species Name - include ionization and fine structure.’ 
print, ’15 characters maximum and DO NOT leave any proceeding blanks.’ 
print, ’Enter for example:’ 

print,’ Fe II* for iron singly ionized and excited Fine Structure’ 

print,’ H2 L 7,0 P(5) — Lyman By stem of molecular Hydrogen’ 

print,’ ’ 


if rtd eq 67 then print, 'Default: ’ ,»tring(Bdata.«l) 
vdum » * ’ 
read , vdum 

if (vdum «q ") AID (rtd eq 67) then •lm®8tring(sdata.el) «ls« •lm=vdum 
•data. el ■ 0*adata.al 

•1b ■ strtrim(elm, 1) + ’ ' 

•data. el ■ byte(strmid(elm,0,15)) 

iafsv ■ long (sdata. iaf ) 

i if ■ longCedata.iaf ) - long(1000)*FIX (sdata. iaf/1000) 

iaf ■ FH(iaf) 

vdum ■ * ’ 

if rtd *q 67 then print, ’Default: ’ ,iaf/10 
print, ’What is the IONIZATION state [ 2 « II 4 - IV 3 ? 
read, ’If a molecule, then what is the VIBRATION number? ’.vdum 
if (vdum «q ”) AND (rtd eq 67) then IA=iaf/10 else IA*fix(vdum) 
rdum * ’ ’ 

if rtd *q 67 then print, ’Default: ’ ,iaf- 10* (iaf /10) 

print, ’What is the FINE STRUCTURE state [ 3 * ***]?> 

read, ’If a molecule, then rhat is the ROTATION number? ’.vdum 

if (vdum eq ”) AND (rtd eq 67) then FS=iaf-10*(iaf /10) else FS=fix(vdum) 

iaf-FIX(IA*10 + FS) 

vdum * strmid(elm,0,2) 

if (vdum eq ’H2’) OR (vdum eq ’HD’) then begin 
ly « ” 

read, ’Is this part of the Werner Bystem’ ,ly 
if (ly eq ’y’) OR (ly eq ’Y’) then iaf = iaf + 1000 
•ndif 


if rtd eq 67 then print, f ormat=’ (a8,f 10.3) ’, ’Default: ’ .sdata.wl 
READ, ’What is the Rest WAVELENGTH of the line (Angstroms)? ’ .vdum 
if (vdum ne ’’) OR (rtd ne 67) then sdata.rl * float(vdum) 

if rtd eq 67 then print, ’Default: ’ .sdata.f 

READ, ’What is the OSCILLATOR stength of the line ?’,vdum 

if (vdum ne ’’) OR (rtd ne 67) then sdata.f * float(vdum) 

if rtd eq 67 then print, format*’ (a8,f 10.3) ’, ’Default: ’ .sdata. owl 
READ, ’What is the OBSERVED WAVELENGTH of the line (Angstroms) ?’ ,vdum 
if (vdum ne ”) OR (rtd ne 67) then sdata. owl * float (vdum) 

if rtd eq 67 then print, ’Default: ’ .sdata. eqw 

READ, ’What is the EQUIVALENT WIDTH (Angtroms) of the line ?’,vdum 
if (vdum ne ’’) OR (rtd ne 67) then sdata. eqw * float(vdum) 

if rtd eq 67 then print , ’Default :’ .sdata. me 

READ, ’What is the ERROR of the measurement ?’,vdum 

if (vdum ne ’’) OR (wtd ne 67) then sdata. me * float(vdum) 

if rtd eq 67 then print , ’Default :’ .sdata.f m 

read, ’What is the First Moment (Velocity km/s) ?’,vdum 

if (vdum ne ”) OR (rtd ne 67) then sdata. fm = float(vdum) 


if rtd ®q 67 then print,’ Default: ’ .sdata.fme 

r®ad,’What is the ERROR in th® First Moment T’.Tdum 

if (vdum n® ”) OR (rtd n® 67) then sdata.fme » float(vdum) 

if rtd ®q 67 th®n print, ’Default: ' .sdata.sm 

read, ’What is th® S®cond Moment (km/s/s)?’ ,rdum 

if (vdum n® ’’) OR (rtd n® 67) then sdata.sm * float(vdum) 

if rtd ®q 67 then print, ’Default: ’ .sdata.sm® 

r*ad,’tfhat is th® ERROR in th® Second Moment T’.rdum 

if (vdum n® ”) OR (rtd ne 67) then sdata.sme » float(vdum) 

if rtd ®q 67 then print, ’Default: ’ .string (sdata.com) 

read, ’C om ments on Measurement (Up to 10 characters) ?’,vdum 

if (vdum ®q ”) AMD (rtd eq 67) then com=string(sdata.com) else com«vdum 

•data.com * 0* 8data.com 

com * 8trtrim(com, 1) + * ) * 

sdata.com ■ byte(strmid(com,0,10)) 

vdum * ’ ’ 
print , * ’ 

if rtd ne 67 then sdata.up * 0.*sdata.up 

read, ’Stop to change UserParameters? (Y/N Default is N) ’ ,vdum 
print , ’ * 

if (vdum eq ’y’) OR (vdum eq ’Y’) then begin 
up * sdata.up 

print , 1 The UserPaxameter is a 30 point floating array. 
print.'The program has Btopped and values may be entered interactively, 
print,’ For example: the statement "up (2) ■ 31.", puts 31 in the 
print,’ third element of the array.’ 

print, ’Mote: the addresses range from 0 to 29.’ 

print, ’The command: ".con <Retum>" must be entered when finished." 
stop 

sdata.up * up 
endif 


Encode as Periodic Table/Molecules 


tst * byte(strmid(elm,0,2)) 

•Is - *H HeLiBeB C N 0 F NeNaMgAlSiP S ClArK CaScTiV CrMnFeCoNiCuZn’ 

•Is « ®ls+ ’ H2H2HDHDC0CHCN C2 OHUiUiUi * ! Elements + Molecules 

els * byte (els) 


j - o 

atst * 0 

while atst eq 0 do begin • 

atst - ((tst(O) ®q els(j)) AND (tst(l) eq els(j+l))) 

j ' j + 2 

if j gt 78 then goto.exitst 
•ndwhile 

exitst: sdata.iaf * long(iaf) + long(1000)*(j/2) ; 


Search for its order. 


Include ionization. 


RETURN 


EH) 


pro mantau, STAR, date 

;******************************* HAHTAU.PRO ******************************** 
; By Chariot L. Joooph Juno, 1990 

0 

1 

; To manipulate the flux, wavelength, continuum, and optical-dapth vectors. 


HC » 0 
MKA » 0 

kcnt « 0 

kdev ■ get env( ’TERM’) 

tas * a88oc(5,fltarr(2,600)) 

ebtst * 0 

cnct ■ 1 

vshft * 0. 


HardCopy Flag ■ HO. 
Hake ASCII Flag ■ HO. 
File access counter. 
Get terminal type. 

For random disk access 
Ho error bars plotted. 
Connect the dots plot . 


for k*0 ,4 do print,’ ’ 
print , ’ 
print , ’ 
print , ’ 
print , ’ 
print , ’ 
print , ’ 
print , ’ 
print , ’ 
print , ’ 

getrdy .STAR.niaf .rwav.lwf s ,bgs ,bge ,n_opt ,c_opts ,B_opt .ctd.opt , iafb4 .knone 


The basic philosophy of this routine is as follows. The researcher 
selects either new species or new profiles to be examined. For each 
selection, there are options as to how the data will be displayed. 
Plotting options affect the current and all subsequent data. Thus, 
changes in the display style from one profile to the next should be 
imposed after the profile has been selected. 


aax * [0. ,0. ,-0.25,0. ,0.25] 
aay * [ 0 . , 3 . 0 , 2 . 6 , 3 . 0 , 2 . 6 ] 

STRT: wtd ■ wmenu(c_opts ,title=0,init 

if wtd eq 4 then stop 
if wtd eq 4 then goto, STRT 
if wtd eq 7 then begin 
vtmp * ebtst 

if vtmp eq 0 then ebtst • 1 
if vtmp eq 1 then ebtst * 0 
if ebtst eq 1 then c_opts(7) 
if ebtst eq 0 then c_opts(7) 
goto, STRT 
endif 

if wtd eq 8 then begin 

tenet * wmenu(ctd_opt , title* 
if (tenet le 0) OR (tenet gt 
cnct * tenet- 1 
if cnct eq 2 then ebtst * 0 


Values for an arrow 
(Lower Limit) . 

= 1 ) ; Input What -To-Do Flag 

Stop - Temporary Halt. 

- Start again. 

Toggle plot error bars. 

; Plots to have errors. 

; Stop plotting errors. 
’Toggle Error Bars OFF’ 

’Toggle Error Bars ON’ 

; Go 8 tart again. 

; Diff. Connect-Dots? 
0,init*l) ; Get choice. 

3) then goto, STRT ; Silly choice-ignore , 
; Set current status . 

; Line only -> Ho errors . 
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if ebtst eq 0 than c_opts(7) - ’Toggle Error Bars OH* 
goto, STRT '• Co ,tart «« ain ’ 
•ndif 

if vtd eq 9 than begin 

print, string (7B),” ! *“8 BELL, 

print,’ Mova cursor to thia work area’ 
read, ’and enter Telocity offset (km/s) : ’ ,wshft 
goto, STRT I Go ,tart »8 ai11 ' 

•ndif 


if wtd eq 11 than begin 

userprog5,STAR,niaf .rwav.lwfs.bgs.bge.tau.kcnt ; Call user’s program 


goto, STRT 
•ndif 

if (wtd It 1) OR (wtd ge 12) then goto, EXIT 
if wtd gt 2 then wtd « 2 

symnum ■ intarr (19) 
snum ■ intarr (19) 

iTvnm « intarr (19) 

ksav » intarr(19) 

pleb ■ intarr (19) 

cplt ■ intarr(19) 

▼s ■ fltarr(19) 

spei ■ string(80) 


Exit Options. 

Temp, void options. 

; Zero plot symbol table 
Zero species table. 

Zero profile table. 

To save record #’s. 

Plot which Error Bars? 

For Connect-dots plot . 

For velocity shifts. 

For window label. 


START: ; — — - — Starting place for 1st species 

L ■ 0 t E * 0 J Auto Plot ranges. 


B * 0 A T * 0 „ v 

TKUM « 0 I Initialize Image Number 


if (wtd eq 1) OR (wtd eq 2) then begin 

which * wmenu(n_opt ,title=0 , init=l) 
if which eq knone then goto, EXIT 
snum(INUM) * which 
iaf * iafb4(which) 


; Start with 1st data. 

; Get species choice. 
; Choice: "No More". 

; Save Species Number 


•pei « «... ’+n_opt (which) 

plotconfig,0,spei, * ’ .O.kdev, ’large’ ; Open type 0 window, 

endif 


Next Get, Header, Optical Depths + Spectra 


LOOP: ; 
lch « 1 

nent « where (niaf eq iaf) 
sz * size(ncnt) 

if sz(0) eq 0 then goto, BRNCH 
; if sz(l) eq 1 then kent ■ fix(ncnt(0)) 
if sz ( 1) ge 1 then begin 

l_choice B string(bytarr(22,sz(l)+2)) 
l_choice(0) * ’Which Line? Log fW’ 
l_choice(sz(l)+l) * ’None’ 
for k=l ,sz(l) do begin 

•wav ■ string (format®’ (f8.3) ’ ,rwav(ncnt(k-l))) 


; Just in case no choices 
; Find all matches. 

; Total # in sz(l) . 

; Problem, Ho Match. 

; Only one profile. 

; Several profiles of the 
; same species setup to 

; make a selection. 

; For each profile: 

-list rest wavelen. 
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*» (f5.3) ' ,lwfs(ncnt(k-l))) 
’♦siwf : 


-list log f-lambda, 
•put into KENU . 


slwf ■ »tring (format* 
l.choice(k) » swav+* 

•ndfor 

lch « smenu(l_choice,title*0,init*l) ! Cat KEBU profila choice 

if (lch la 0) OR (lch ga sz(l)+l) then goto, BRNCH ; lot valid aalaction. 
kcnt • ncnt(lch-l) i choice into variab. 

andif 


lnum(IBUM) - lch 
ksav(IBUM) * kcnt 
if cnct na 2 than symnum(IBUM) ■ 
if symnum (IFUM) eq 9 then begin 
cnct ■ 2 
abtst » 0 

if abtat aq 0 than c_opte(7) ■ 


; Store profile number. 

; Store for later plots. 
wmenu(s.opt,title*0,init«l) ; Symbol choice? 

; Line OBLY: 

; ■> Set to Line. 

; *> Bo errors. 

'Toggle Error Bars OB’ ; ■> Reset Menu opt. 


•ndif 

if abtst than pleb(IBUM) « 1 ELSE pleb(INUM) ■ 0 
cplt(IBUM) * cnct 
▼s(IBUM) « vshft 


Read A Plot Optical Depths 


PLT: asy ■ findgen(16)*(!PI*2/16.) 

if HC then plot config.l,’ ’,-l.kdev,” 


; Used for plot symbol. 
; Config. for Hardcopy? 


for k*0,IBUM do begin 
kcnt * ksav(k) 
close, 5 

openr,5.STAR+’ .TAU’ 
tautmp * fltarr(2,600) 
tautmp * tas (3*kcnt+l) 
sz ■ TOTAL (tautmp(0,*) ne 0) 
tau ■ fltarr(6,sz) 

tau(0: l,0:sz-l) = tautmp(0:l,0:BZ-l) 
tautmp ■ 0.* tautmp 
tautmp * tas(3*kcnt+2) 
tau(2:3,0:sz-l) * tautmp (0:l,0:sz-l) 
tautmp * 0 . * tautmp 
tautmp * tas(3*kcnt+3) 
tau(4:5,0:sz-l) * tautmp(0:1.0:sz-l) 
ttmp * where (tau (2,*) gt 1.0, tent) 
mtmp * where(tau(3,*) gt l.O.mcnt) 
close, 5 


Make temporary array. 

Get wavelength A TAUs(-). 
Find the length. 

Make suitable TAU array * 
stuff Wavelength A TAUs(-) 
Zero temporary array . 

Get TAUs and TAUs(+) A 
put into TAU array. 
Zero temporary array. 

Get Spectrum A Continuum 
and put into TAU. 
Which are limits? 

Which are limits? 


lwf .fact - 14.567 - lwf s (kcnt) Convert to Column Density. 

tau(l:3,*) * tau(l:3,*) + lwf.fact 

swav - string(format=’ (f8.3) ’ .rwav(kcnt)) ; Rest wavelength string. 


if MKA eq 1 then begin 
slwf * string (format*’ (f 6.3) ’, lwf s (kcnt)) 
print ,ewav+ 9 log f -lambda: , +Blvf+ , 
printf , 1 ,swav+ * log f -lambda: * +slwf+' 
sbg * stringCf ormat s * (f 9 . 1) * ,bgs (kcnt) ) 


; IF Make ASCII file. 
< — log f -lambda. 

1 +n„opt (snum(k) ) 

,+ n_opt (Bnum(k) ) 

< — Background. 
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•bge ■ string(foimat«’ (f9. 1) ’ ,bge(kcnt)) ; <~ BG Error, 

printf ,1, 'Background: '+sbg*’ BG Error: '+sbge 

printf , 1 , ' ’ 

printf , 1 , 'Wavelength log N(min) log N log B(maz) Spact. 
for kt»0,sz-l do printf ,1, format*' (f9. 3, 5f 11.2) ' , $ 

tau(0,kt) ,tau(l ,kt) ,tau(2,kt) ,tau(3,kt) ,tau(4,kt) ,tau(5,kt) 
printf, 1,* ’ 

•ndif 


Cont. » 


if vshft na 0. than bagin 
print , string (7B) , ' ’ 

print, 'Caution: profile shifting should seldom be done’ 
print,' and then only with extreme care.* 

print, 'The *+swav+ 'profile has been shifted by: ’.vshft,’ km/s.’ 

•ndif 

▼0 ■ 2.9979e5/rwav(kcnt) 

v ■ ▼0*(tau(0,*)- rwav(kcnt)) + vs(k) ; Calculate velocity. 


if symnum(k) eq 3 then begin 
usersym.O .7*cos (asy) ,0.5*sin(asy) ,/FILL 
tsym * 8 

end else begin 

usersym,0.7*cos(asy) ,0.5*sin(asy) 
tsym * symnum(k) 
if tsym eq 9 then tsym * 0 
endelse 

ymsav * ly.margin(O) 
yht « 0.22 - 0.04*k 


; Make special symbols? 
Convert dots to filled 
circles . 

; Option for open circles 
; Most cases, take IDL’s. 
Solid line option. 


Just in case HC ■ 1, be 
ready to make a key . 

If cplt(k)=2 -> LINE 

*1 -> LINE+DOTS. 
*0 -> DOTS. 


if (k eq 0) AND (cplt(k) It 2) then begin 
if (H ne 0) OR (B ne 0) OR (L ne 0) OR (T ne 0) then begin 

plot , v ,tau(2 , •) ,xtitle= ’ velocity ’ ,ytitle= ’ log N ' ,psym=tsym, /YNOZ , $ 
xrange= [L ,H] ,yrange= [B.T] ,xstyle=l ,ystyle=l 
•ndif ELSE begin 

plot,v.tau(2,*),xtitle='velocity’,ytitle=’log N’ ,psym=tsym,/YNOZ 


endelse 

if HC then !y .margin (0) * 4 

if HC then plots, [0.15,0.20] , [yht, yht] ,psym=tsym, /NORMAL 
if HC then xyouts, 0. 25, yht-0. 01, l.choice (lnum(k) ), /NORMAL 
!y .margin(O) * ymsav 
•ndif 

if (k ne 0) AND (cplt(k) It 2) then begin 
oplot ,v,tau(2,*) ,psym=tsym 
if HC then !y.aargin(0) ■ 4 

if HC then plots, [0. 15,0.20] , [yht, yht] ,psym=tsym, /NORMAL 
if HC then xyouts, 0.25, yht-0. 01, 1. choice (lnum(k)) ./NORMAL 
!y.maxgin(0) * ymsav 
•ndif 

if (k eq 0) AND (cplt(k) eq 2) then begin 
if (H ne 0) OR (B ne 0) OR (L ne 0) OR (T ne 0) then begin 

plot,v,tau(2,*) ,xtitle= ’velocity ’ ,ytitle='log N’,/YN0Z, $ 


,yrange» [B,T] ,xstyle=l ,ystyle*l 
•ndif ELSE b«gin 

plot , tau(2 , *) ,xt itle« ’ velocity * .ytitle- » log N * . /YIOZ 

•ndelse 

If HC than !y.margin(0) • 4 

if HC then plots, [0.15,0. 20], [yht.yht] ./NORMAL 
if HC then xyouts, 0 . 25, ybt- 0 . 01 , l_choice(lnum(k)), /NORMAL 
!y .margin(O) » ymsav 
•ndif 

if (k no 0) AND (cplt(k) eq 2) then begin 
oplot ,v,tau(2,*) 
if HC then !y.margin(0) » 4 

if HC then plots, [0. 15,0.20] , [yht.yht] ,psym=tsym. /NORMAL 
if HC then xyouts, 0. 25, yht~0. 01, l_choice(lnum(k)) , /NORMAL 
!y .naxgin(O) • ymsav 
•ndif 

if cpltCk) eq 1 then begin 

oplot, v,tau(2,*) I Connect the points opt. 

if HC then !y. margin (0) * 4 

if HC then plots, [0.15,0.20] , [yht.yht] , /NORMAL 
if HC then xyouts, 0. 25, yht-0. 01, l_choice(lnum(k)) , /NORMAL 
!y. margin (0) * ymsav 
•ndif 

if pleb(k) then oplot, v, tauCl, *) ,psym=tsym, Bymsize=0.6 
if pleb(k) then oplot, v, tau(3, *) ,psym=tsym, 8ymsize=0.6 


usersym.aax.aay ; Over plot any Lower Limits 

if tent gt 1 then oplot, v (ttmp) ,tau(2,ttmp) ,psym=8 

if tent eq 1 then oplot, [v(ttmp) ,v(ttmp)] , [tau(2, ttmp) ,tau(2, ttmp)] ,psym=8 
if pleb(k) then begin 

if ment gt 1 then oplot ,v(mtmp) ,tau(3,mtmp) ,psym=8 

if ment eq 1 then oplot, [v(mtmp) ,v(mtmp)D , [tau(3,mtmp) ,tau(3,mtmp)] , $ 

psym=8 

end if 


andfor 

if HC then begin 

spe * ’Input File: ’+STAR+’ .TAU* 
xyouts ,0. 15,0 .99 , spe , /NORMAL 
spe * n_opt(l) 

xyouts ,0.50,0. 99 ,n_opt (snum(k)+l) , /NORMAL 

xyouts ,0.80,0.99, date , /NORMAL 

plot conf ig , - 1 , spei , ” , - 1 ,kdev , ’ large ’ 

HC * 0 
•ndif 

if MKA eq 1 then begin 
close, 1 
MKA « 0 
•ndif 


; If in HardCopy Mode: 

; Put File Name on header. 


Species Name on header . 

Put Date on header. 

send to plotter k config 
no longer HardCopy . 


If Making ASCII file: 
close the file, 
no longer Making ASCII. 


» 


Kuttr Branching Section inside LOOP 


HUTCH : ; 

wtd • wmenu(c_opts,title*0,init*2) 
if wtd eq 4 than stop 
if wtd eq 4 then goto, BRTCH 
if wtd eq 7 then begin 
wtmp ■ ebtst 

if wtop eq 0 then ebtst * 1 

if wtop eq 1 then ebtst * 0 

if ebtst eq 1 then c_opts(7) • 
if ebtst eq 0 then c_opts(7) ■ 
if ebtst eq 1 then pleb(INUM) * 
goto, PLT 
endif 

if wtd eq 8 then begin 

tenet * wmenu(ctd_opt ,title*0,init*l) 
if (tenet le 0) OR (tenet gt 3) then goto, PLT 

if (cnct eq 2) AND (tenet It 3) then $ 

eymnum(IHTJM) ■ wmenu(s_opt ,title=0,init*l) 
cnct * tcnct-1 
cplt(INUM) * cnct 
if cnct eq 2 then ebtst * 0 
if ebtst eq 0 then c_opts(7) = 
if ebtst eq 1 then pleb(INUM) 
goto, PLT 1 

endif 

if wtd eq 9 then begin 
print ,string(7B) , ” 
print, ’Move cursor to this work area’ 
read,’ and enter velocity offset (km/s) : ’ ,vshft 
vs(INUM) * vshft 
goto, PLT 
endif 

if wtd eq 10 then begin 
print ,string(7B) , ’ * 

print, ’Move cursor to this work area.’ 
print, ’The current value is : ’ ,lwfs (kent) 
read, ’Enter new value of log (f-lambda) ’ .lwftmp 
lwfs(kcnt) * lwftmp 
goto, PLT 
endif 

if wtd eq 11 then begin 


call COG Menu. 

Stop - Temporary Halt. 

- Start again. 
Toggle plot error bars. 


; Plots to have errors. 

; Stop plotting errors. 
’Toggle Error Bars OFF’ 

’To gg le Error Bars OH’ 

1 ELSE pleb(INUM) « 0 

; Replot k Bring up Menu. 

Diff. Connect-Dots Status? 
Get choice. 

Silly choice - ignore. 

From line to with "dots"? 
Symbol choice? 

Set current status It keep 
record for the future. 
Line only => Ho errors . 
’Toggle Error Bars DN’ 

1 ELSE pleb(INUM) * 0 


Replot k Bring up Menu. 


Ring the BELL . 


Save for later plots. 
Go start again. 

Adjust f value? 

Ring the BELL. 


Save for later plots. 


userprog5 ,STAR,niaf ,rwav ,lwf s ,bgs ,bge ,tau,kcnt ; 
goto, BRNCH J 

endif 

if wtd le 0 then goto, BRNCH 

if (wtd It 1) OR (wtd ge 12) then goto, EXIT ; 

if (wtd eq 1) then begin l 

vshft * 0. • 

goto, START 
endif 


Call user’s program 
- Start again. 


Exit Options. 

Get NEW 1st data set. 
Velocity Bhift default. 
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_ ^r(wtd aq 2) than begin 
^ I SUM ■ ISUM ♦ 1 

if ISUM gt 18 than begin 
print .string (7B) 

print, ’WARNING: Program can only 
ISUM » 18 
•ndif 
wahft » 0. 

■nun (ISUM) - which 

iaf ■ iafb4 (which) 
goto, LOOP 
•ndif 

if wtd aq 3 then begin 
HC » 1 


; Add sore sets of data. 

; Incr. Image Bumber. 

; Too many profiles! 

; Ring warning BELL, 
handle 19 profiles.’ 

; Sat to back to mar. 

; Velocity shift default. 


; Go get new set of obs. 
; Want HardCopy? 


goto, PLT 
•ndif 

if wtd aq 5 then begin 

READ, ’What are Xmin.Imax .Ymin.Ymax? ’ ,L,H,B,T 
goto, LOOP 
•ndif 

if wtd aq 6 then begin 
ASC: afile « ” 

print, string (7B) , ” J 

print.’Move cursor into this Workspace and’ ; 
read,’ Enter OUTPUT file name: afile ; 

close, 1 

openr.l, afile, ERROR * errtst ; 

close, 1 

if errtst aq 0 then begin \ 

print, string(7B) , ’Warning there is a file with that name. 


Adjust plot limits? 


Make ASCII File of data? 

Ring the BELL. 

Give directions. 

Get File Name. 

Is there one already? 

If file exists 


yno ■ ’ ’ 

read, ’Would you like to choose another name? (no *> destroy old)’, yno 
if (yno aq ’Y’) OR (yno eq ’y’) then goto, ASC 


•ndif 

openw, 1 ,af ile 
MKA * 1 
goto, PLT 
•ndif 

if (wtd eq 5) OR (wtd eq 6) then goto, LOOP 


Open the file. 

Set flag to print. 

GO cycle through data 


END Branching Section 


EXIT: close, 5 
; !x.margin(0) * 10. 
wdelete,0 


; Close file and return. 

; Left plot margin->default. 
; Delete Existing window. 


RETURN 

END 


GETRDY is a aatup Program for MANTAU 


by Charles L. Joseph 6/1/90 


STAR 

nlaf 

rwav 

lwfs 

bgs 

bge 

n_opt 

c_opts 

•_opt 

ctd_opt 

iafb4 


knone 


file same of the .TAU file, 
am array holding the list of species codes . 
an array holding the list of rest wavelengths . 
an array holding the log f "lambdas . 
an array holding the background levels, 
an array holding the BG errors, 
string array holding the different species, 
string array holding the program Control options . 
string array holding the Symbol options . 
string array holding the Connect-The-Dots? options, 
an array holding the distinctly unique species codes, 
differs from niaf in that niaf may have multiple 
entries of the same species, 
is the comparison value of the choice indicating none. 




***** 


pro getrdy , STAR , niaf .rwav , lsf s ,bgs ,bge ,n_opt , c_opts , s.opt . ctd.opt , iaf b4 , knone 


as ■ assoc(5,fltarr(200)) 
close ,5 

openr , 5 , STAR+ ’ . TAU ’ 
niaf * as CO) 
rwav » as(l) 
bgs * as (2) 
bge * as (3) 
dates * as (4) 
lwfs * a8(5) 
close ,5 


Association variables for 

Open file for Update t get 
1st record - species codes 
Get addresses in .DTL file 
Get previous BG levels. 

Get previous BG errors. 

Get previous date codes . 
Get prev. log f -lambdas. 


els - *H HeLiBeB C N 0 F HeNaMgAlSiP S ClArK CaScTiV CrKnFeCoNiCuZn’ 
els « els+ ’H2H2HDHDC0CHCNC20HU U U’ ; Elements + Molecules 

gfg S ********** X 1 

sii - 'I II III IV V VI VII VIIIIX X ’ 


plotconfig.O, ’Optical Depths’,’ ',0, 
i.opt * string(bytarr(26,99)) 
iaf s ■ where (niaf ne 0.) 

■z * size(iafs) 

if (sz(0) eq 0) then ntl ■ 0 ELSE n 
iafb4 * fltarr(99) 
kk = 1 

for k ■ 1, ntl do begin 
nxt * FIXCiaf s(k-l)) 
iaftst * niaf (nit) 
tst * where (iaftst eq iafb4) 
sz * size(tst) 

if (sz(0) eq 0) then begin 
iafb4(kk) « iaftst 


iv, ” ; Set up for type 0 window. 

; Can handle 99 species. 

* BZ(1) 

; Search through obs. data. 
; Get address of next value. 

; Put in test variable and 
; see if there are other 
; identical entries . 

; If not, add this species 
; to the list. 
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lent ■ fix (2 . * (niaf (met) / 1000 . -1) ) 

ii ■ niaf(nxt) - 1 000. *FIX (niaf (nxt)/1000) 

fa * ii sod 10. 


ii « f ix(ii/10) 
siaf • ’ ’ 

if ient It 60 than bagin 
if ii la 0 than ii » 11 

if (fa gt 0) AID (fa la 9) then f a-atittidCafa.O.fa) ELSE fa 
aiaf ■ atrmid(aii,4*(ii-l),4)+fs 


» » 


andif 

if ient ga 60 than bagin 

aiaf ■ ’Tib: ’+string (format*’ (i2) ’ ,ii) 
aiaf • aiaf+’ rot : ’+string (format*’ (i2) ’ ,fs) 
andif 

i.opt(kk) * atrmid (els , ient ,2) + ’ ’+siaf 

^ m vk ♦ 1 ; Incr. f of valid entries. 


andif 


andfor 

i_opt(0) ■ ’Species Options’ 
i.opt(kk) * ’No more’ 
knone * kk 

n_opt * i_opt(0:kk) : Take only aub-array. 

iafb4 « iafb4(0:kk) 


Load the Control Menu Options 

c.opts ■ atring(bytarr(26,13)) 
c_opts(0) * ’Options:’ 
c.opts (1) * ’Get New Species’ 
c.opts (2) * ’Add Another Profile* 
c.opts (3) ■ ’Send Plot to Laser Printer’ 
c.opts (4) * ’Stop - Temporary Halt’ 
c.opts (5) * ’Adjust Plotting Limits’ 
c_opts(6) * ’Make ASCII file of Data’ 
c_opts(7) * ’Toggle Error Bars ON’ 
c.opts (8) * ’Change PLOT Format’ 
c.opts (9) ■ ’Shift Rest Wavelength’ 
c_opts(10) * 'Change log (f-lambda) ’ 
c.opts (11) * ’Run User Program 5’ 
c.opts (12) * ’Quit and Return’ 

-------------- Load Plotting Symbol Options 

a.opt * string(bytarr(16,10)) 
a.opt(O) ■ ’Plotting Symbols’ 
a.opt (1) * ’Plus sign’ 
a_opt(2) ■ ’Asterisk’ 
a.opt (3) ■ ’Filled Circle’ 
a.opt (4) * ’Diamond’ 
a.opt (5) » ’Triangle’ 
a.opt (6) ■ ’Square’ 


^*pt(7) - 'X* 

*»_opt(8) ■ 'Open Circle’ 
•.opt (9) - 'Lino - no pts 


I 


Load Connect the Dot* Options 


ctd_opt ■ •tring(bytarr(24,4)) 
ctd.opt(0) » 'Connect the DATA points?’ 
ctd_opt(l) * 'Plot only as DOTS’ 
ctd_opt(2) ■ 'Plot as DOTS and CONNECT’ 
ctd_opt(3) ■ 'Plot data only as a LINE’ 


RETURN 

END 

PRO INTGRT , Y , X , CNT ,TAU , JE , IF , XI . XLIMIT .mdata ,mp .mcntrl .up 

I*********************************** INTGRT *********************************** 

J By Charles L. Joseph Jan. 12. 1984 June. 1990 

; DETERMINES EQUIVALENT WIDTHS OR INTEGRATED FLUXES 

\ OF ABSORPTION OR EMISSION FEATURES 

• A TRAP0Z0IDAL METHOD IS USED TO COMPUTE THE INTEGRATION 

flux vector 
wavelength vector 
continuum vector 
epsilon vector (data quality) 
index of cutbot 

final wavelength value allowed for in fits or integration 
initial or first wavelength value ( XI“left, XF=right) 
structure containing mp.wtd *1 exit 


Y 

X 

CNT 

EPS 

JE 

IF 

XI 
np 


************** 


*************************************************************** 


HC « 0 

xxxs * fltarr(2) 
yyys ■ f ltarr (2) 
sz * size(JE) 

if sz(0) eq 0 then lcs ■ 1 else lcs * sz(l) 
TmT ■ max (XLIMIT, J3) 
tbit * max(x,J2) 
jmx * fix(J2/2) 

mp.SNR * (CNT (jmx) -mp . bg) /mcntrl . ESAV 
jmx * J2 
Tmn * min(x) 

while (I(J2) gt xmx) do J2=J2-1 
XLIMIT (J3) * KJ2-3) 
kcnt * 0 


HardCopy Flag -> NO. 

For Overplot markers . 
For Overplot markers . 
Determine # of profiles 
to be measured. 

Get longest wavelength. 

Temporary use of jmx. 
Signal-to-Noise Ratio. 


; Find A adjust XLIMIT* s 
; Exclude last 3 points. 

; * of pts in integration 
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size (I) 

jcat • 1 

if bz( 1) gt 16 then jcat ■ f ix(sz(l)/16) 


; Size for BG overplot. 


IF mcntrl . intopt EQ 0 THEN BEGIN ; ««<««««« ; Use nonstndrd ALT.INT? 

The following 6 print statements may NOT be removed legally. 

print,’ ’ ,string(7B) 

print, string (7B) , 'No Alternate Integration Routine is Provided’ 
print, ’in standard M.S.L.A.P. - version 1.0’ 
print , ’ ’ 

print, ’Origin of the Alt. Integration Routine is as follows:’ 
print , ’ ’ 

ALT.INT . X , Y , CNT , TAU , JE , XF , XI , k cnt , XLIMIT .rndat a ,mp .mcntrl 

auto_int_sav ,X , Y ,CNT ,mp ,mcntrl ,up .XLIMIT ,XI ,XF ; Default call: A dummy. 

END ELSE BEGIN ; <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< ; Use Standard EQUIVW. 


kkcnt * 0 

for k*0,lcs-l do begin 

•data - mdata(k) ; Get next mult. ID. 

print, string (sdata. el) ■ Print species name. 

if lcs gt 1 then xxf=XF(k) else xxf=XF ; Get next stopping point 

if lcs gt 1 then xxi=XI(k) else xxi*XI ; Get next starting point 

if lcs gt 1 then jje=JE(k) else jje=JE ; Get next line center. 

; mp.SNR « 1. /(mcntrl. ESAV/CNT(j je)) 

mp.SNR * (CNT ( j j e) -mp . bg) /mcntrl .ESAV ; Calc. S/N from Contim. 

EQUIVW , X , Y . CNT , TAU , j j e , xxf , xxi , kcnt , XLIMIT , sdat a ,mp .mcntrl , up 
auto_int_sav.X,Y. CNT, mp, mcntrl, up, XLIMIT. xxi.xxf ; Default: A dummy, 
if (kkcnt+kcnt gt 0) then tautmp * fltarr(6,kkcnt+kcnt) 
if k eq 0 then begin 

if kcnt gt 0 then tautmp * TAU ! Save Tau’s, if any. 

END ELSE BEGIN 

if kkcnt gt 0 then tautmp(0:5,0:kkcnt-l) • tausv 
if kkcnt gt 0 then tautmp(0:5,kkcnt:kkcnt+kcnt-l) = TAU 
•ndelse 

kkcnt * kkcnt + kcnt 

if kkcnt gt 0 then tausv * tautmp ; Save any Tau’s. 

•data. up ■ up ; fold UserPar into data 

mdata(k) « sdata ; Put into mult. data, 

•ndfor 

kcnt - kkcnt : Tid Y U P counter and 

if kcnt gt 0 then TAU « tautmp ; optical depths array. 

ENDELSE ; <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< \ Standard EQUIVW. 


Plot Results Section 
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PLTSTRT: if id.siaa eq ’PS’ than device, /CLOSE ; 

ti* m * Heasureaent Results ’ ! T 

plotconf ig,BC, ' * ,tle,2,kdev, ’ * i c 

; if HC eq 1 then device , /ESCAPSULATED , FILES AME* ’ intgrt2.ps ’ 
; TiUa »> profile areas : 

•jjsh m min(tau(4,0:kcnt-l) , J5) * ^ 

yix • aax(Y) » 

if (yam It 0.2*yax) then pbg * 1 ELSE pbg • 0 ; S 

if pbg then begin * G 

wan ■ ap.bg ~ ap.bgerr * c 

if vmn It 0 then van ■ 1 . l*vmn l * 

if yan It van then van ■ yan * ^ 

yaz » l.l*yax * J 

xxxs(O) ■ xan t xxxs(l) * xmx i I 

yyys(O) ■ van t yyys(l) * ymx ; 


Title for window. 
Configure graphics, 
i’ 

Find ain and aaz flux 
to see if plot BG. 
Set plot BG flag. 

Going to show BG+error 
Get ain. needed for BG 
Adjust for aore space. 
Less than ain. spect? 
Adjust for aore space. 
Put limits into small 
over plot array A 


; let IDL make axis. 

plot .xxxs.yyys ,/BODATA, xtitle=’ Wavelength (A)’, ytitle=’Rel. Flux’ 
if ap.SNO then oplot ,I,SM00TH(Y,9) else oplot,X,Y ; Bow really plot data. 

van - where (x It tau(0, J5) . jbg) » Get r * ad y t0 over P lot 

if (jbg gt jcnt) then xxxb(O) * x(jbg-jcnt) ELSE xxxs(0)*x(0) ; BG range, 

if (jbg It jax-jcnt) then xxxs(l)=x(jbg+jcnt) ELSE xxxB(0)*x(jmx) 
yyys(0: 1) * mp.bg 

oplot, xxxs.yyys • 0ver P lot BG ‘ 

yyyB(0:l) * ap.bg - mp.bgerr 

oplot. mi ; Ovor plot BG-BGerror. 
yyys(0:l) * ap.bg + mp.bgerr 

oplot, xxxs.yyys, linestyle=2 S ° ver P lot BG+BGerror. 

endif 

if pbg ne 1 then begin » »° T g oin 6 t0 8how BG ‘ 

if mp.SMO then plot.X,SM00TH(Y,9) ,/YNOZ, xtitle= ’Wavelength (A) ’ , $ 
ytitle*’ Flux’ else plot,X,Y,/YNOZ,xtitle=’Wavelength (A) ’ ,ytitle=’Rel. Flux’ 
endif 

if mp.SHD ne 1 then oplot ,X,Y,psym=8 

oplot. X. CUT 5 Put on continuum. 

plotlab2 ,X , Y ,CNT ,mp .mcntrl ,up i Customized labels. 

xxo*0.0*(!x.crange(l)-!x.crange(0)) + ix.crange(O) 
yyo=l .05*(!y.crange(l)-!y.crange(0)) + !y.crange(0) 

T71TAU • Put on File name. 

xyoutB^xxo.yyo^mp.FNAM » ruc 

xxo*0 . 65* ( ! x . crange ( 1 ) - ! x . crange (0) ) + fx.crange(O) 

xyouts ,xxo,yyo,systime(0) 

if kcnt gt 0 then begin 

for kk«0,kcnt-l do begin Show areas that were 

xxxs(0: 1) * tau(O.kk) • integrated. 

yyys(O) ■ tau(4,kk) 
yyys(l) ■ tau(5,kk) 
oplot ,xxxs ,yyys 
endfor 
endif 

xdum * where (XLIMIT gt O.nxs) 

yout * 0 . 07* ( ! y . crange ( 1) - ! y . crange (0) ) + !y.crange(0) 

for kk-0.nxs-1.2 do begin ; Show the areas that 


Put on continuum. 
Customized labels. 


Put on File name. 


Show areas that were 
integrated. 


Show the areas that 
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■yyys(O) ■ !y.crange(0) 
yyys(i) » ly.crange(l) 
xxxs(0: 1) - XLIMIT(kk) 
oplot ,iui ,yyys ,linsstyl*-l 
i (0:1) - XLIMIT(kk+l) 


oplot ,xxxs ,yyys ,linestyle*l 
xxxb(O) - XLIMIT(kk) 
yyye(0:l) » yout 
oplot .xxxs ,yyys 
«xx if or 


Set to plot top to 
bottom . 

ware us«d to calculate 
polynomail lit of the 
continuum. 


; Connect solid line for 
; continuum areas. 


ly.margin(O) ■ 1 


; Enable text area. 


•pe b 'Measurements in file: ’+mp.STAR+’ .DTL’ 
xyoute ,0.05,0.32, ape , /NORMAL 
■pe ■ 'Real S/H * ’+string (format*’ (f6.1) ' ,mp.SNR/(mp.cohfac“0.5)) 
xyouts , 0 . 765 , 0 . 32 . spe , /NORMAL 

■pe - string (format*’ (f6.1) ’ .mp.SNR) ; Start building Btring. 

■pe * 'Apparent S/N Ratio: ’+spe+’ based on ’ 

■pe « spe+Btring (format*’ (i3) ’ , mcntrl. NNE)+’ points with a Noise ’ 

■pe ■ spe+’ Coherence Length of ’+string(format*’ (f4. 1) ’ .mp.cohf ac) 
xyouts, 0.05, 0.30, spe, /NORMAL : Print S/N results. 

■pe * Btring (format*’ (i2) ’ ,mp.poly) ; Start building string, 

■pe ■ ’Continuum was fit with a polynomial of order: ’+spe+’ ’ 

■pe ■ ■pe+string(format*’ (f9.2) ’ ,mcntrl.f In) +8tring (format*’ (f9.2) ’ ,mcntrl.f2n) 


Create Btring. 
Print File name. 


xyouts ,0.05,0.28, spe , /NORMAL 
■pe * ’BG (background) was taken to be:’ 

■pe * spe+string (format*’ (f 12. 1) ’, mp.bg) +’ +/“' 

■pe * spe+string (format* • (fl2.1 ) 9 ,mp .bgerr) 
xyouts ,0.05, 0.26, spe, /NORMAL 
•data * mdata(0) 
spe * 'For ’+string(sdata.el) 
xyouts , 0 . 05 , 0 . 24 , spe , /NORMAL 
•pe * stringCf ormat= * (f 8 .3) ’ ,sdata. wl) 
xyouts ,0.24,0.24, spe , /NORMAL 
spe * 9 error contributions from BG:' 
spe * spe+string(f ormat- * (f 9 .4 ) 9 ,mcntrl . eO) 
spe 

xyouts ,0.34,0.24, spe , /NORMAL 
spe * 9 Errors belov axe the Addition in Quadrature of* 
spe c spe+' the Background and RMS-Noise Errors* 
xyouts ,0.05,0.22, spe , /NORMAL 
spe * ’Species 1 


Start building Btring. 


Print BG t uncertainty 


BG Error contribs . 


Start building 
a string. 


spe+string (format* ’ (f 9 . 3) ’ .mcntrl . e 1) +string(f ormat* ’ (f 9 . 3) ’ .mcntrl . e2) 


BG Error contribs 
Build string for a 
note on errors. 
BG Error contribs . 


xy out s , 0 . 05 , 0 . 18 , spe , /NORMAL 
spe * ’Lab. Wave.’ 


xyouts ,0 .21 ,0. 18, spe, /NORMAL 
spe * ’f’ 

xyouts ,0.35,0. 18 , spe, /NORMAL 

spe « ’ Obs. Wave* 

xyouts , 0 . 43 , 0 . 18 , spe , /NORMAL 

spe • 1 EQW (A) 1st (km/s) 2nd (km/s/s)* 

xy out s , 0 . 57 , 0 . 1 8 , sp e , /N ORM AL 


«0,lcs-i do begin 
•data ■ mdata(k) 
yht ■ 0.16 - 0.04*k 
■pa ■ ■tring(adata.el) 

■pa " •trtri»(»pe,2)+' * 

■pa ■ atxmid (ape, 0,15) 

xyouta , 0 . 05 , yht , spa , /NORMAL 

■pa ■ string (format*’ (f8.3) ’ .adata.wl) 

xyouta, 0.21 ,yht , spa, /NORMAL 

■pa ■ »trimg(format«’ (f6.4) ’ ,adata.f ) 

xyouta , 0 . 33 , yht , apa , /NORMAL 

■pa « «tring(format“’ (f8.3) * .adata.owl) 

xyouta , 0 . 43 ,yht , ape , /SORMAL 

■p« ■ string(f ormat*’ (f 10.4) * ,sdata.eqw) 

xyouta, 0.55, yht, ape, /NORMAL 

■pa ■ *tring (format*’ (f 10.3) ' ,adata.fm) 

xyouta ,0.67, yht , ape , /NORMAL 

■pa ■ atring (format* * (f 10.3) ’ , sdata.sm) 

xyouta ,0 . 82 , yht , ape , /NORMAL 

yht « 0.14 - 0.04*k 

xyouta, 0.47, yht , ’Errors : ’ , /NORMAL 

■pa * string (format** (f 10.4) ’ .adata.me) 

xyouts ,0.55, yht , ape , /NORMAL 

■pa ■ atring(f ormat* ’ (f 10.3) * .adata.fme) 

xyouta ,0 .67 , yht , ape , /NDRMAL 

apa * atring (format*’ (f 10.3) ’ .adata.sme) 

xyouta ,0.82, yht , ape , /NORMAL 


For each Baaaurament: 

- Gat Measurements . 

- Ad j . print locat . 


; - Print momenta. 

; - Print momenta. 


- Build a string. 

; - Print momenta. 

; - Print momenta. 

; - Print moments. 

; - Adj . print locat . 

; values . 

; - Print errors. 

; - Print errors. 

; - Print errors . 


andfor 

if HC eq 1 then plotconf ig.-l , ’ ’,-l.kdev,” 

if HC aq 0 then begin 

HC ■ wmenu( [’ Copy? ’ , ’ Yes ’ , ’No ’] ,title=0, init=l) 
if HC aq 1 then goto, PLTSTRT 


; Send plot k -> term. 

; If not already making 
; a hardcopy, want to? 

; If yes , go do it . 


endif 

EXIT: CLOSE, 4 
■delete, 2 
RETURN 

END ; INTGRT 


; close it up 
; Delete extra window. 
; go back to DETAILS 


PRO EQUIVW , X , Y , C ,TAU , JE , XF , XI , kcnt , XLIMIT , sdat a ,mp .mcntrl ,up 

j ********************************** EQUIVW ******************************** 

; COMPUTES INTEGRATED FLUX - ABSORPTION OR EMISSION LINES 
« 

; TRAPOZOIDAL METHOD 

t 

; I wavelength vector 

; Y flux vector 

• C continuum vector 
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JE 

IF 

XI 

kcnt 

ILIHIT 

•data 


■P 

mcntrl 

up 


Index of cursor position 
right dsf . and stop point 
loft dsf. and stop point 
Busbar of points ussd in tbs intsgration 
cursor positions ussd to dsf ins continuum rsgions 
■tructurs containing tbs ID plus momsnts and aixorB 
•tructurs MSLAPparameter 
•tructurs for sons control ▼ariablos 
UssrParameter, a 30 point floating vector which is 
rsssrvsd sxclusivsly for tbs user . 


Translated to IDL by C. JOSEPH 7/1/89 from code by E. Jenkins 


********** 


J1 « 0 

while (I(J1) It XI) do Jl*Ji+l 
ttht * max(x,J2) 

•bile (X(J2) gt IF) do J2=J2-1 
kcnt * J2 - Ji + 1 
10 ■ 0.0 
BKG * mp.bg 
CF2 * 9.0 
C0HFAC ■ mp.cohfac 
EBKG - mp.bgerr 
V * fltarr(5,3) 

EM * fltarr(5) 

ER » fltarr(5) 

EE * EM 
F » EM 

A * fltarr(6,600) 

XSAV « X 

▼elf act * 2 . 9979o05/sdata . wl 
X * velfact*(X-sdata.wl) 
xje * X(JE) 

X - X - X(JE) 

QDELT « (X(J2) - X(J1))/FL0AT(J2 - Jl) 
for II ■ J1,J2 do begin 
D • (C(II) - Y(I«) 

B - C(II) - BKG 
1(0,11) * D/(B + EBKG) 

1(1,11) - D/B 
A(2,II) « D/(B - EBKG) 

XI ■ X(II) - X0 
12 * X1‘2 

for J * 0,2 do begin 

V(0,J) * V(0,J) + A(J,II) 
V(l.J) * V(l.J) + A(J,II)*X1 
V(2, J) - V(2,J) + A(J ,II)*X2 


Find left edge of integrat 
ing region. 

Find right edge of integr- 
ating region. 


Set Coherence Factor. 

Set background error est. 
create array: V 
create array: EM 
create array: ER 
create array: EE 
create array: F 
create array: A 


For each integration point: 
Calculate Continuum - Flux 
and Continuum - Background. 
Use these to find the tract 1 
ion of missing flux, taking 
into account BG errors . 
Calculate weighting for 
1st and 2nd moments. 

Update 0th, (equiv. width) 
1st, and 

2nd moments for the 


throe BG coses. 


endf or 


•ndfor 

for J » 1,3 do begin 

VCl.J-1) ■ V(1,J-1)/V(0,J-1) ; Formalize let t 2nd mom. 

V(2, J-l) - V(2,J-1)/V(0,J-1) - V(1,J-1)“2 
endf or 

EM(0) • FLOAT (J2 - J1 + O/COHFAC 

for j-1,4 do EMC j ) » EM(j-l)*EM(0) 

for J ■ 0,4 do F(J) « ABS(V(J,2) - V(J,0))/2. 

F(0) « F(0)*QDELT 
JCEHT - (J1 + J2)/2 
SI0C2 ■ l./(ap.SSR"2) 

SC0C2 » SI0C2/FL0AT(mcntrl.FNE) 

T1 - SI0C2*EH(0) 

T2 ■ SC0C2*(EM(0) - V(0,1)/C0HFAC)“2 
EE(0) « C0HFAC*SQRT(T1 + T2) 

VB2 - ((V(l.l) - I(JCENT) + X0)/QDELT)“2 
V3Q « V(2,1)/QDELT“2 

EE(1) « COHFAC*SQRT (SI0C2* (EM (2) *CF2/ 12 . + EM(0)*VB2) $ 

♦ SC0C2*EM(1)*VB2)/V(0, 1) 

Q * EM(2)*CF2/12. ♦ EM(0)*VB2 

EE(2) « COHFAC*SQRT(SIOC2*(EMC4)*CF2‘2/80. + EM(2)*CF2*VB2/2 . $ 

+ EM(0)*VB2"2 - 2.*V3Q*Q $ 

+ EM(0) *V3Q“2) + SC0C2* (Q - EM(0)*V3Q)“2)/V(0,1) 

EE(2) « EE(2)*QDELT 
for J ■ 0,4 do begin 
EE( J) « EE(J)*QDELT 
ER( J) ■ SQRT(EE(J)*2 + F(J)‘2) 


•ndfor 

V(0,1) « V(0,1)*QDELT 
•data.eqw * V(0, l)/velf act 
edata.me * ER(0)/velfact 

■data. owl » (V(l,l)+xje)/velfact + sdata.wl 

edata.fm ■ V(l,l)+xje 

edata.fme * ER(1) 

edata.sm * V(2,l) 

sdata.sme ■ ER(2) 

mcntrl.eO * F(0)/velfact 

mcntrl.el * F(l) 

mcntrl.e2 * F(2) 


Save Equivalent Width (A) 
Save its error. 

Save observed wavelength. 
Save first momement. 

Save its error. 

Save 2nd momement . 

Save its error. 

BG contrib. to EQW error. 
BG contrib. to 1st error. 
BG contrib. to 2nd error. 


Calculate the Optical Depths 


X « X + X(JE) I Change origin of vel. back. 

SCOC « SQRTCSC0C2) 

B « C(J1: J2) - BKG 
D « C(J1:J2) - Y(J1:J2) 

A(2,0: J2-J1) « D / B ; TAU’s for no BG error. 

CCOR * C*SC0C 
D * D - CCOR 
B * B - CCOR 

A(1,0: J2-J1) = D / (B + EBKG) 


; Minimum TAU’s. 



2 . *CCOR 
2 . *CCOR 


A(3,0: J2-J1) ■ D/(B - EBKG) 
A « 1.0 - A(*,0: J2-J1) 


A - -1.0*ALOG(A > 4. *-5) 


A » AL0G10CA > 0.005) 

X ■ XSAV 

A(0,0: J2-J1) - X(J1:J2) 
A(4,0: J2-J1) • Y(J1 : J2) 
A(5,0: J2-J1) - C(J1:J2) 

TAU • A 

RETURN 


END ; EQUIVW 


; Mar inram TAU* 8. 

; Limit the ring*. 

; Return X to wavelength. 
; Wavelengths for TAU ’ b . 

; Store Spectrum. 

; Store Fit Continuum. 

; All info. -> TAU array. 
; back to IWTGRT 


jhig routine sets plotting configuration throughout M.S.L.A.P. It sizes 
the plot windows, correctly adjusting for the environment. It also directs 
the plots to various terminals or hard copy devices. 


HC HardCopy flag -1 Close Hardcopy device A Bet to terminal. 

0 Terminal or console. 

1 HardCopy Unit - currently PostScript. 

tlO Title for window 0. 

tl2 Title for window 2. 

wtype open new window of various types . 


-2 «> No New window - reset to default plot parameters. 

- ALSO : Do not produce CURSORS . 

«> Ho New window - reset to default plot parameters. 

0 *> window 0 — the large, main plotting window. 

1 => windows 0*1 - MAIN plus ONE for the cursor coords. 

2 => window 2 the expanded results window. 

3 => window 3 small instructions window, then do 

an immediate return. 


kdev Terminal Device Type either a ’sun’ or ’xterm’ is returned, 

xmarg Sets !x.margin(0) , the left edge of the plot. 

marg * ’large* •*> more space to WMENU functions. 

Note: xmarg ■ ’large’ overrides the wtype variable. 


pro plotconf ig,HC,tlO,tl2, wtype , kdev .xmarg 


if wtype eq 3 then begin 

window, 3, title=” ,xpos=450,ypos=B30,xsize=450,ysize=40 
return 
endif 


Only open window 3 
; for intructions. 
RETURN. 


kdev = get env( ’TERM’ ) 

if kdev ne ’sun’ then kdev=’ xterm’ 


; Get terminal type. 

; Assume X window or SUN 


if HC It 0 then begin 
device, /CLOSE 


Finish HC plotting. 
Close device/PS-f ile . 


; Move the plot file. 

; Send to laser printer 
; Flag to terminal . 


•pawn.’mv -f idl.ps twmp.pB’ 
»pawn,*lpr tamp.ps’ 

HC « 0 
•ndif 


if (HC eq 0) AID (kdev eq ’xterm’) then begin ! If wcreen A X window. 

■•t_plot, , X’ 

if .q 0 thu b.gln : °=U 

window, 0,color»2,ipo#*200,ypo8»350 ( i«iz**620,yii2*»550,titl«*tl0 

•ndif 

if wtype eq 1 then begin » wi * dow8 Banted - 

window.O, color-2, xpos-200,ypos=350,xsize=820,ysize*550,title*tl0 

window, l.xpos-770,ypos»150,xsize-200,ysize*50. title- ’Cursor Coordinates’ 

B8et0 ; Ready for window 0. 

•ndif , . , 

if wtype eq 2 then begin J Expanded results wind. 

window, 2. color-2, xpos-250,ypos-50,xsize-770,ysize=820.title-tl2 

!y .margin(O) - 30 » Reserve text area ’ 

•ndif 

if wtype It 0 then !y.margin(0) * 4 » Default under graph, 

if wtype It 0 then !x. margin (0) - 10. I Default left margin. 

•ndif 

if (HC eq 0) AND (kdev eq ’sun’) then begin ; If screen t Sunviev, 

■et_plot , * SUN ’ , 

if wtype eq 0 then begin > 0nl y MAIN window - 

window.O, color-2, title=tl0,xpos=100,ypos=270,xsize=1030,ysize=600 

•ndif . 

. , 4 • _ ; Two windows wanted, 

if wtype eq 1 then begin * * 

window.O, color-2, title=tl0,xpos=100,ypos=270,xsize=1030,ysize=600 
window, 1, title- ’Cursor Coordinates’ ,xpos-820.ypos=180,xsize=200,ysize-50 

„ et ,o ; R « ad y for sindow °- 

•ndif , . , 

if wtype eq 2 then begin i Expanded results wind, 

window ,2 . color-2 ,title=tl2 ,xpos=250 ,ypos=50 .xsize-880 ,ysize=820 

!y .margin (0) - 21 ! Reserve text 

•ndif J . 

if wtype It 0 then !y.margin(0) « 4 * Default under grap . 

if wtype It 0 then !x.margin(0) - 10. ; Default left margin. 

•ndif 


if xmarg eq ’large’ then begin 
!x.margin(0) - 30. 

if kdev eq ’xterm’ then !x.margin(0) « 50. 
•ndif 


Extra room for HMENU. 
for Sunviev window. 
X window on SUN. 


if HC eq 1 then begin 
set.plot, *PS’ 


; If HardCopy Unit. 
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. /INCHES . Y0FFSET*2 . 8 , YSIZE-7 . 5 
“"^fy .margin (0) • 24. 

fx.margin(O) • 10. 

•ndif 


; Adjust plot area. 

; Reserve text ar«a. 

; Use all hor. space. 


if wtype aq -2 than raturn 


; Avoid m a kin g cursors. 


if BC aq 0 than device, /CURSOR.CROSSHAIR 
if HC aq 0 than TVCRS, 0.5, 0.5, /NORMAL 


; Graphical cursors A 
; placa on scraen. 


RETURN 

END ; plotconfig 

;**************************** COG. PRO ************************************* 

9 

; CURVE OF GROWTH - for workstations 

; - Input input files are from MSLAP 

I 

; By Charles L. Joseph: 5/22/79 

; Latest Modification: 10/19/90 


:****** 


*********************************************************************** 


PRO COG, STAR, date, libr 

nulld = { noda, el: bytarr(15), iaf: long(O), wl: 0., f: 0., vel: 0., aqw: 0., $ 
me: 0., fm: 0., fme: 0., sm: 0., sine: 0., com: bytarr(lO), up: fltarr(30) > 
dtl - replicate({ noda >, 200) ; Make array for storage 

kdev ■ getenv(’TERM’) 

FIREUP , IAF , DV , STAR , dtl ,n_opt , c.opt s . s_opt ,nns , iaf b4 , cch , knone 


wtd « wmenu(c_opts,title=0,init=l) 
if wtd aq 4 then stop 

if (wtd It 1) OR (wtd ge 7) then goto, EXIT 
if wtd gt 2 then wtd * 2 

gymnum = intarr(9) 
snum « intaxr (9) 
cdsav * fltarr(9) 
melsv ■ fltarr(9,2,20) 
cogi ■ string(80) 

STRT: ; Starting place for 1st species 

INUM * 0 

aqlsv » fltarr(9,20) 
wfsav * fltarr(9,20) 
if (wtd aq 1) OR (wtd aq 2) then begin 

which * wmenu(n_opt,title=0,init=l) 
if which aq knone then goto, EXIT 
snum(INUM) ■ which 
iaf * iaf b4 (which) 
tst ■ where (cch aq iaf) 
sz * size(tst) 

if (sz(0) aq 0) then cnum=3 else cnum=tst(0) 
openr , 1 , libr+ ’ /coginf o . tab ' 
for k*0,cnum do readf.l.cogi 


Input What-To-Do Flag 
Stop - Temporary Halt. 

Exit Options. 

Temp, void options. 

; Zero plot symbol table 
Zero species table. 

Column Density mstr tab 
Mean Error mstr table. 


Initialize Image Number 
log(EQW) master table. 
log(Wf) master table. 

Start with 1st data. 

; Get species choice. 
Choice: "No More". 

Save Species Number 

Computer’s Choice for 

the theo. CogNumber 
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^-^fose.l 

cwf ■ float(strmid(cogi,71,7)) 
print, cwf 

plotconf ig.O.cogi, ’ \0,kdev, ’large’ ; Open type 0 aindow. 

; fx.margin(O) -30. 1 Mai ® P lot square. 

I if kdev eq ’ item ’ then !x.margin(0)*50. > X window a little diff. 

lynnum(INUM) » wmenu(B_opt ,title=0,init«l) ; Get symbol choice. 

endif 

whichcogs ,nna ,IS ,NT ! C.O.G. ’s to use? 

T « 0. t B - 0. ft L - 12. ; Set plot for autoscale. 

START: ; ---------- Starting place for Additional Species ------- -------- 

iaf « long(O) : i6f code id ’ 8 8 P ecie8 ' 

print, iaf b4 (which) , long ( iaf b4 (which)) 

if (which gt 0) then iaf • long(iafb4(which)) 

GRAB, STAR, C, IAF, WF, EQL, MEL. WAF,F,W,EQW,NL, ML, ME.L.dtl.cnum 

wfsav(INUM.O:NL-l) ■ WF : Update master log(Wf ) . 

eqlsv(INUM,0:NL-l) « EQL • Update master log(EQW) . 

melsv(INUM ,0 : 1 ,0:ITL-1) * MEL 

IF NL LT 0 THEN GOTO, EXIT 


LOOP: 


Master LOOP to shift Obs. relative to theo. C.O.G. 


LINES, A, NT.NS.M.H.L.T.B, WO, EQL, libr, which ; DRAWS C.O.G. 

plotem , wf sav , eqlsv , symnum ,melsv ,n_opt , snum , INUM ; Overplot Obs. data. 
PRINT * 9 

print, ’Use left mouse to locate a starting data position for translation’ 
print, ’Use right mouse to bring up a MENU of options.’ 

CURSOR, XI, IY 

if kdev eq ’xterm’ then wait.l ; X window? Slow down. 


Master Branching Section inside LOOP 


if !ERR eq 4 then begin ’• If ri g ht mouse * then 
wtd * wmenu(c opts,title=0,init c 2) • call COG Menu, 

if wtd eq 4 then stop J Stop - Temporary Halt, 

if wtd le 0 then goto, LOOP 

if (wtd It 1) OR (wtd ge 7) then goto, EXIT ; Exit Options. 

if (wtd eq 1) then begin I Get lst data set ’ 

tst = j ; OK, if not much data. 

if INUM gt 1 then begin ! Check before data loss. 

print ,string(7B) 

print, ’Warning: Previous Column densities will be erased’ 
tst * wmenu( [’Confirm* , ’Yes ’, ’No’) ,title=0, init=2) 
endif 

if tst eq 1 then goto.STRT *. Confirmed, start over, 

goto. LOOP J Continue on. 

endif 

if (wtd eq 2) then begin '• Xdd more sets of data. 

which - wmenu(n_opt,title=0,init=l) ; Get species choice, 

if which eq knone then goto, LOOP * Choice: ’No More . 

INUM ■ INUM + 1 » * ncr - I® a £ e Number, 

snum (INUM) * which 


Exit Options. 

Get NEW 1st data set. 
OK, if not much data. 
Check before data loss. 
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^ ^ symnum(IWUM) • wmenu(s_opt,title«0,init«l) ; Get symbol choice. 

goto, START 5 Co K Bt E8W * et of ob8 - 

•ndif 

if wtd eq 3 then goto, HC 5 Co make hardcopy? 

if rtd eq 5 then $ * Adjuet plot Units? 

RElD,’Vhat are Iiain,Xmai,Ymin > YmaLX? > ,L,H,B,T 
if rtd eq 6 then whichcogs ,nns ,NS,NT ; Change # of curves? 

if (rtd eq 5) OR (rtd eq 6) then goto, LOOP 

endif ; - — — END Branching Section — 


print , * * 

print, 'Use left souse to indicate new location for the translation.’ 
CURSOR, IX, IY 

if kdev eq ’xtenn’ then wait.l ; X window? Slow down. 

VF-VF+IX-I1 

wfsav(INUM,0:NL-l) « WF 

LINES,A,NT,NS,M,H,L,T,B,WO,EQL,libr .which ; DRAWS C.O.G. 

plotem,wfsav,eqlsv,symnum,melsv,n_opt,snum,INUM ; Overplot Obs. data. 

print, ’Total SHIFT: ’ ,rf (O)-waf (0) 
cdl * rf(0) - waf(0) + cwf 
cdsav(INUM) • cdl 

print, ’Log Column Density: ’,cdl,’ for current set of observations’ 
goto, LOOP 


HC: PRINT,’ ’ 

; ! I. margin (0) =10. ; Set left margin back, 

plotconf ig,l, ’ ’.’ ’,-l.kdev,’ ’ ; Set up for hardcopy. 

LINES, A, NT.NS.M.H.L.T.B, WO, EQL.libr, which ; DRAWS C.O.G. 

if nt gt 5 then xyouts, 0.7,0. 66, ’b Values ’, /NORMAL $ ; Position A write 
else xyouts, 0.7, 0.6, ’b values ’, /NORMAL ; Doppler broadening vals. 

for k=NS .NS+NT-l do begin 

if nt gt 3 then yht * 0.63-0. 02* (k-NS) 
if nt le 3 then yht * 0.57-0. 02* (k-NS) 
xyouts ,0.69 ,yht ,DV (k) + ’ km/ s ’ , /NORMAL 


endfor 

plotem , wf sav , eqlsv , symnum ,melsv ,n_opt , snum , INUM 
!y.margin(0) * 4 

spe * *0bs. data from file: ’+STAR+’ .DTL’ 
xyouts ,0.1,0.28, spe , /NORMAL 
xyouts ,0.85,0. 28 , date, /NORMAL 
xyouts ,0.1,0.26, cogi , /NORMAL 
for k=0,INUM do begin 
yht * 0.22 - 0.02*k 
spe ■ strtrim(n_opt(snum(k))) 
spe ■ 'log N(’+spe+’)’ 
xyouts ,0 . 1 ,yht , spe , /NORMAL 
spe * ’ * ’+strtrim(string(cdsav(k))) 
xyouts ,0.35 ,yht , spe , /NORMAL 


; Overplot Obs. data. 
Enable text area. 

Create string. 

Print Data file name. 
Print date of analysis. 
Print Theo. COG info. 
For each species, write 
log N, column density. 
Get species name. 

Print species. 

; Print column density. 


endfor 

IF NL LE 7 THEN STP=NL ELSE STP=7 

plotconf ig,-l, ’ ’,’ ’ ,-l , kdev, ’large’ ! Send plot t -> terminal. 

; !x.margin(0) * 30. \ Make plot more square. 
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if kdev eq *xter®’ than !x.*argin(0)«50. 
GOTO. LOOP 


EXIT: wdelete.O 

; !x.margin(0) ■ 10. 

RETURN 


ETC) 


; X window a little diff. 
; Set back to dafault. 


PRO FIREUP , IAF.DV , STAR , dtl ,n_opt , c_opts , a_opt ,nns . iaf b4 . cch .knone 
; m*m**mm*m*****m* FIREUP. PRO ********************************** 

; TO DETERMINE WHICH C.O.G. IS TO BE USED. 

; AND TO SET VARIOUS PARAMETERS. 

; by Charles L. Joseph 5/19/79 

.*********************************************************************** 
! ******** OPEN THE INFO FILE 


close, 1 


openr, 1, STAR* ’ .DTL’ , ERROR * errtst 
close, 1 
app - ’ ’ 

if errtst eq 0 then begin 
CLOSE, 1 

OPENU , 1 , STAR* ’ . DTL ’ 

ML- 200 
readu, 1 ,dtl 

iafs « where(dtl.iaf ne 0.) 
sz * size(iafs) 

if (sz(O) eq 0) then ntl « 0 else 
close, 1 
endif 

plotconf ig,0, ’Curve of Growth’,’ ’ ,0 

i.opt « Btring(bytarr(26,99)) 

iafdat ■ dtl. iaf 

iafs ■ where (iafdat ne 0) 

iafb4 » long(intarr(99) ) 

iondat * dtl. el 

kk - 1 

for k * 1, ntl do begin 
nxt - iafs(k-l) 
iaftst « iafdat (nxt) 
tst - where (iaftst eq iafb4) 
sz - size(tst) 

if (sz(0) eq 0) then begin 
iafb4(kk) - iaftst 

i_opt(ki) « string (iondat (0:9, nxt)) 
kk « kk + 1 

end if 
endfor 


; Does file exists? 


; Get old file. 

; Get f of lines, 
ntl * sz(l) 

,kdev,’ ’ ; Set up for type 0 window. 

; Can handle 99 species. 


; Search through obs . data. 
Get address of next value. 

Put in test variable and 
see if there are other 
identical entries. 

; If not, add this species 
; to the list. 

; Incr. # of valid entries. 


i_opt(0) » ’C.O.G. Choices’ 
i.opt (kk) ■ ’>o more’ 

knone ■ kk 

n_opt • i.opt (0:kk) ; Take only sub-array. 

iafb4 - iafb4(0:kk) 

— Load Table of Available Tbeo. C.O.G. 

cch » intarr(5) 
cch(O) * 26020 
cch(l) • 12020 
cch (2) - 26020 
cch(3) ■ 8010 

cch(4) » 14020 

LOAD THE DOPP VEL VALUES 

DV -string (BYTARRC3, 11)) 

DV(1) - *1.0* 

DV(2) * ’1.5’ 

DV (3) « ’2.0’ 

DV (4) ■ ’3.0’ 

DV(5) » ’4.0’ 

DV(6) « '6.0’ 

DV (7) ■ *10.’ 

DV(8) « ’15.’ 

DV (9) * ’20. ’ 

DV(10)“ ’30.’ 

Load the COG Menu Options 

c.opts - string (bytarr (26 ,8)) 

c.opts(O) - 'Options:* 

c.opts (1) * ’Compare lBt Species to COG’ 

c.opts (2) * ’Add another Bpecies to COG’ 

c.opts (3) * 'Send Plot to Laser Printer’ 
c_opt8(4) * ’Stop - Temporary Halt’ 
c.opts (5) * ’Adjust Plotting Limits’ 
c.opts (6) « ’Adjust # of Theo. Curves’ 
c.opts (7) * ’Quit and Return’ 

— Load Plotting Symbol Options 

s.opt - string(bytarr(16,9) ) 
s_opt(0) * ’Plotting Symbols’ 
s.opt ( 1) * ’Plus sign’ 
s.opt (2) * ’Asterisk’ 
s.opt (3) * ’Filled Circle’ 
s.opt (4) - ’Diamond’ 
s.opt (5) * ’Triangle’ 
s_opt(6) * ’Square’ 
s.opt (7) * ’X’ 


Fe II COG 
Mg II COG 
Mn II COG 
0 I COG, default 
Si II COG 
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• 'Open Circle’ 


h&b 


;******** 

RETURN 

END 


— -- To number of curves options 

string (bytsrr ( 1 1 , 1 1 ) ) 


nns(l) 

m 

>< 

i > 

nns(2) 

m 

9 < 

2 > 

nns(3) 

m 

»< 

3 > 

nns(4) 

m 

»< 

4 > 

nns(5) 

m 

*< 

5 > 

nns(6) 

m 

»< 

6 > 

nns(7) 

m 

>< 

7 > 

nns(8) 

w 

’< 

8 > 

nns(9) 

m 


9 > 

nns(10) 

m 

'< 

10 > 


END OF FIREUP 


PRO LINES,A,NT,NS,H,H,L,T,B,X,EQL,libr,shich 

.**************♦**************** LINES. PRO ********************************** 


TO GET AND PLOT THE CURVES FOR C.O.G. 


by Charles L. Joseph 5/17/79 

****************************************************************************** 


if (T eq 0.) AND (B eq 0.) then begin 
tst * maz(EQL) 
if (tst It -5) then begin 
T * -3.0 


; If NOT Scaling Override 
; Test for only weak k 
; Set plot limits, 

; accordingly . 


B « -7.0 


L ■ 10. 

H - 16. 

end else begin : strong-lines-case 

j m -2. ; plotting limits. 


B * -6. 

L * 11.5 


H « 18. 
endelse 

endif 


ERASE 

A * libr+'/tabdata/feii2382.tab’ 
A * libr+’/tabdata/siiil260.tab* 


0PENR.1.A 

y = fltarr(lO) 

yy * fltarr(10,41) 
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fltarr(4l) 
f or k-0,40 do begin 

readf .l.foimaf (f 9. 3,10(113.4))’. x,y 


yy(0:9,k) - y 
rx(k) - x 
•ndfor 
close , 1 
IA-NS+NT-1 
FOR M-NS.NA DO BEGIN 
y • yy(M-l,0:40) 

IF M EQ IS THEN PLOT,xx,y,xrange»[L,H] ,yrange*[B,Tj , 
XTITLE* ’ ! 31og (Nf ! 4k ! 3) ’ , YTITLE* ’ ! 31og (W/!4k!3)’ 

IF M GT NS THEN OPLDT.xx.y 

END 

; XYOUT , 400 , 30 , * log N ( cm -2 ) ’ 

; YIOUT, 0,600,’ log eqw/ lambda’ 

CLOSE, 1 
; !PSYM=1 

M*0 


RETURN 

END 


$ 


PRO GRIB, STAR, C.IAF.WF.EQL, MEL, WAF.F ,W,EQW,NL,ML,ME,L,dtl,cnum 
j********************* GRAB. PRO ********************************** 
; TO GET THE EQW DATA FROM THE .DTL FILE 

; WRITTEN BY C. JOSEPH 1980 

; MODIFIED BY C. JOSEPH AND T. ARMITAGE ON 19 MAY 1983 

TO HANDLE A SINGLE OBSERVED EQUIVALENT WIDTH 


**************************************************************** 


OPENR, 1 ,STAR+’ .DTL’ 
readu, 1 ,dtl 
close, 1 

iafs ■ where (dtl.iaf eq iaf) 
NL * size(iafs) 

NL * NL(1) 


; Open data file and read 
; observed data. 

; Get addresses of desired data. 
; Get Number of Lines and put 
; value in NL. 


W * dtl.wl 
F - dtl.f 
EQW ■ dtl.eqw 
ME ■ dtl.me 


; Get those portions of the 
; .DTL file to be used. 


I 
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; Strip out only those measure 
; ments for the give species. 

; If there is only one measure 
; ment, make 2-elm . array of it. 


WAF ■ ALOG10(W*F) 

EQL - AL0G10CEQW/W) 

HEL«FLTARR(2.NL) 

FOR N«0,NL-1 DO BEGIN 

IF ME(N) NE -1000 THEN MEL(0,N) « (EQW(N)+ME(N))/W(N) 

IF ME(N) NE -1000 THEN MEL(l.N) * (EQW(N)-HE(N))/W(N) 

IF HE(N) EQ -1000 THEN MEL(0,N) ■ EQW(N)/W(N) 

IF ME(N) EQ -1000 THEN MEL(1,N) ■ (EQW(N)-17)/W(N) 

END 

KEL=AL0G10(KEL) 

SS * L + 1 - WAF*(F GT 0) \ Shift points so they are on 

WF=WAF+MIN(SS) ’• the C0G P lot • 

FIN: RETURN 
END 


****************************** PLOTEM .PRO ****************************** 
To overplot the various species equivalent widths plus error bars. 

By Charles L. Joseph S/15/90 

*********************************************************************** 


; Calculate log f -lambda’s. 

; Take log of equivalent widths 



1 then begin 

W(iafs) 

F(iafs) 

EQW(iafs) 

NE(iafs) 


W(0:1) * 
F(0:1) » 
EQU(0: 1) 
ME(0: 1) 

NL « 2 

endelse 


pro plot em , wf s av , eqlsv , Bymnum ,melsv , n_opt , snum , INUH 
asy * findgen(16)*(!PI*2/16.) 

DX * fltarr(2) 
for k*0,INUM do begin 
teql ■ eqlsv(k,0:19) 
tst ■ where (teql ne 0,NL) 
teql * teql (tst) 


; Used for plot symbol. 

; For each species, over 
; plot with correct symb. 
", Get addresses of valid 
; data and strip out. 
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"twf • wfaav(k,0:19) 
twf * twf(tst) 
tine ■ melsv(k, 0:1, 0:19) 
tme » tme(0,0: 1 ,tst) 


; Strip out log(fW) . 

I 

; Strip out Mean Errors . 


if symnum(k) eq 3 then begin 
UBersym,0.7*cos(asy) ,0.5*sin(asy) ,/FILL 
tsym ■ 8 

end else begin 

usersym,0.7*cos(asy) ,0.5*sin(asy) 
tsym ■ Bymnua(k) 

endelse 

OPLOT.twf ,teql ,PSYM=tsym 
for N*0,NL-1 do begin 
DX(0: 1) ■ twf (N) 

DY ■ tme(0,0:l,H) 

0PL0T.DX.DY 

endfor 

endfor 


; Make special symbols? 

; Convert dots to open 
; circles. 

; Option for open circles. 
; Most caseB, take IDL’s. 

; Over plot obs . DATA 
; Over plot Error Bars . 


; END Error Bars part. 


if INUM gt 0 then begin 
Y2 

ymg 
xrng 


!y .crange(O) 


If more than one species, 
make symbol key at left. 
Get range of plot so that 
the KEY can be positioned 


fltarr(2) 

!y.crange(l) 

ix.crange(l) - !x.crange(0) 
xoutl ■ O.iOxrng + !x.crange(0) 
xout2 ■ 0. 13*xrng + !x.crange(0) 
if INUM It 5 then tymx =0.75 else tymx=0.85 
DX(0: 1) * xoutl 
for k*0,INUM do begin 

if symnum(k) eq 3 then begin 

usersym,0.7*cos(asy) ,0.5*sin(asy) ,/FILL 
tsym » 8 
end else begin 

usersym,0.7*cos(asy) , 0 . 5*sin(asy) ; Option for open circles, 

tsym » symnum(k) J Most cases, take IDL's. 

endelse 

tdy ■ (tymx-0 . 05* (INUM-k) ) *yrng + !y.crange(0) ; Y-position for next 


Lots of obs.? Adjust TOP 
Stuff 2-elm. X array. 

Make special symbols? 
Convert dots to open 
circles. 


Y2(0 : 1) » tdy 
oplot ,DX, Y2 ,PSYM=tsym 
ape ■ n_opt(snum(k)) 
tdy = tdy - 0.01*ymg 
xyouts,xout2,tdy,spe 
endfor 
end if 
return 
end 


Stuff 2-elm. Y array. 
Over plot the symbol . 


Print species name. 

; END KEY -making part , 


********************* WHICHC0GS .PRO ***************************************** 

To determine starting and how many theoretical C.O.G. are to be used. 

By Charles L. Joseph 5/15/90 


pro whichcog8,nnB,IS,IT 

print , etring(7B) ! » ich C.0 .G.»b to u..? 

print, ’Indicate the total number of theoretical Curves of Growth 
print,’ that are to be plotted. (1 - 10)’ 

print, ’The b values are: 1.0, 1.5, 2, 3, 4, 6, 10, 15, 20, and 30 
print, ’ka/B, respectively’ 
nns(0) » ’# of Curves?’ 
wait , 1 

IT * wmenu(nns,title»0,init«10) 

if IT It 1 then IT « 1 
print , ’ ’ 

print, ’Indicate the curve number for the C.O.G. with the lowest’ 
print, ’b value to be plotted, (e.g. use 2 to get the 1.5 km/s.)’ 
nns(0) * ’Starting?’ 

IS ■ wmenu(nns,title=0,init=l) 

if IS It 1 then IS ■ 1 
if IT gt (11-IS) then IT » 11 - MS 
return 
end 


; Make sure range range 
; makes sense. 


f#**#*«##*****f##******#« END OF COG PACKAGE ##################»«**«*#* 


MSLAP SITE LICENSE AGREEMENT 


TERMS AND CONDITIONS 


1. The title and full ownership rights to each and every part of the Modular Spectral Line Analysis Pro- 
gram (hereinafter MSLAP) shall remain the sole property of Charles L. Joseph and Edward B. Jenkins 
(hereinafter the authors). You shall acquire no rights in the Program other than as expressly granted in 
this Agreement. Each portion of MSLAP constitutes valuable proprietary assets of the authors, em- 
bodying substantial creative efforts and significant expenditure of time. 

YOU MAY NOT USE, COPY, TRANSFER, OR MODIFY THE PROGRAM OR ANY COPY, 
MODIFICATION, OR MERGED PORTION THEREOF, IN WHOLE OR IN PART, EXCEPT AS EX- 
PRESSLY PROVIDED IN THIS AGREEMENT. YOUR LICENSE TO USE MSLAP IS AUTOMATI- 
CALLY TERMINATED IF YOU TRANSFER POSSESSION OF ANY COPY, MODIFICATION OR 
MERGED PORTION OF THE PROGRAM TO ANOTHER PARTY. 

Version 1.0 of MSLAP is licensed free of charge only to sites that are actively engaged in astronomical 
research for the sole purpose of pursuing astronomical research and that are in compliance with the 
terms of this Agreement. Possessing MSLAP in whole or in part indicates that the user has accepted all 
of the terms of this agreement. Violation of any terms of this Agreement shall make you liable for full 
payment for any damages incurred to the authors. 

The provisions of this section 1 shall survive any termination of this Agreement. 

2. The authors provide MSLAP and license its use on a single site bases. A site is defined for the pur- 
poses of this Agreement as a single Department, Center, or Institute that is located on a single campus 
of a university or government agency. Permission is granted to you to use or to copy MSLAP to any 
machine on your site that is controlled by your Department, Center, or Institute, provided such action is 
in fall compliance with all of the terms of this Agreement. 

3. MSLAP IS PROVIDED ’AS IS’ WITHOUT WARRANTY OF ANY KIND, EITHER EXPRESSED 
OR IMPLIED. NO WARRANTY OF ANY KIND IS IMPLIED OF MERCHANTABILITY OR FIT- 
NESS FOR A PARTICULAR APPLICATION. THE ENTIRE RISK AS TO THE QUALITY AND 
PERFORMANCE OF THE PROGRAM IS WITH YOU AND NOT THE AUTHORS NOR THEIR 
AGENTS. 

In no event will the authors or their agents be liable to you for any damages arising out of use or ina- 
bility to use MSLAP, even if the authors or their agents have been advised of the possibility of such 
damages, or for any claim by any other party. If any third-party claims arise against you, you shall 
bare aU such liabilities. 

4. Standard version 1.0 of MSLAP consists of the 34 routines licensed to you that initially are contained 
in the files: main.mslap, master^ux, compare.pro, posto.pro, edatdtl.pro, mantau.pro, inlgrl.pro, 
plotconfig.pro, and cog.pro as well as all of the associated tables found in files with the extension 
".tab". You are permitted to alter this standard source code as you desire for your personal or site s 
use, but you must maintain the original copyright notice in both the source file and in the initial start up 
of MSLAP, and you may NOT distribute the altered code off site without written authorization from 
Charles L. Joseph. If your site has 2 or more visitors in a given year who use MSLAP, you are con- 
sidered a "guest user facility", and you must comply with the conditions spelled out in section 8 of this 
Agreement. Any alterations to the standard source code must also meet the documentation require- 
ments spelled out in section 6 of this Agreement. 


^^a-ofroutines, contained in the file mslap.pro, are designed for individualized customization of 
ldSLAP. You are permitted to create new modules (subroutines) of the USERPROG s or to make 
adaptations of the routines originally found in the file mslap.pro, and to distribute these customized 
modules to other sites, provided the new code meets the requirement in section 6. 

6. Any and all modifications must comply with the following: 

L All modified code must describe all parameters in each procedure definition statement. 


ii. All modified code must indicate the original author, author’s affiliated institution, and a 
modification history including each and every major participant involved in the modification. 

iii. If the altered code is designed to override a particular behavior of standard MSLAP, a concise 
description of the former and new behavior must be output in a conspicuous manner and also in- 
cluded in the source code documentation. 

iv. The ALT INT routine has special rules. This routine is used to provide an alternate type of in- 
tegration. The following disclaimer must appear in a conspicuous manner whenever the routine is 
called. Standard MSLAP with a working instruction window at the lower left of the screen au- 
tomatically produces this disclaimer. The following lines from that disclaimer are listed below 
and may NOT be removed from the file intgrtpro. 

print,’ ’,string(7B) 

print,string(7B),’No Alternate Integration Routine is Provided’ 
print.’in standard M.S.L.A.P. - version 1.0’ 
print,’ ’ 

print.’Origin of the Alt. Integration Routine is as follows:’ 

7. A library of user-generated data-gcUing routines are made available to the astronomical community 
and users are encouraged to add to this library. MSLAP supports any combination of 5 or less of these 
routines at one time. (See the MSLAP Documentation Manuel for the general implementation of dgetl, 
dget2, dget3, dget4, and dget5 procedures and their appropriate protocols.) All "dget" routines are con- 
sidered to be formally outside of MSLAP and not subject to the terms of this Agreement However, it 
is strongly recommended that the documentation rules listed above be implemented to maximize the 
overall utility of any user-contributed routines for data retrieval. 

8. If you represent a guest user facility, defined as any site where more than 2 visitors use MSLAP in a 
given year, you are not permitted to incorporate modifications that remove the modularity of MSLAP, 
thus making it difficult for guest users to customize MSLAP. For example, you must avoid reserving 
more than 10 of the 30-element vector called "up" (UserParameter) since such action would severely 
inhibit the ability of a user to store and manipulate customized calculations. 

Some version of the original file mslap.pro must be retained and be made available to individual users 
for their personal customization of MSLAP. 

If you modify any portion of MSLAP including the original file mslap.pro to override a particular 
behavior of standard MSLAP, a concise description of the former and new behavior must be output in a 
conspicuous manner when MSLAP is running and also must be included in the source code documenta- 
tion. 

9. Permission is granted to incorporate MSLAP into a large driver program, provided the name MSLAP 
and the original copyright notice is maintained in any menu or other solicitation to activate MSLAP. 
Any incorporation into a driver program is also subject to sections 4, 5, 6, and 8. 
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>«r'Tour site is considered to be a Distribution Center and you may transfer copies of the STANDARD, 
version 1.0 of MSLAP to other astronomical sites only if this section of this agreement is complete with 
all necessary signatures. 


Site Name and Address of the Department, Center, or Institute that is to be a a Distribution Center, 


Date and Signature of authorized representative of the above site. 


Name and Position of the person signing above. (Please print or type) 


Date and Signature of Charles L. Joseph, the first author. 

11. The license granted under this Agreement is effective until terminated. You may terminate this 
license by destroying all copies, modifications and merged portions of MSLAP in your possession and 
notifying Charles L. Joseph in writing of such destruction and termination. 

The license granted under this Agreement will terminate if you violate any of the terms and conditions 
of this Agreement. The first author, Charles L. Joseph, reserves the right to terminate this Agreement, 
if he believes that any of the terms and conditions of this Agreement have been violated. Notice of 
such termination shall be made in writing to you. You agree upon such termination to destroy promptly 
all copies, modifications, and merged portions of MSLAP in your possession and to certify to the first 
author that such action has been taken. 

12. None of your rights, duties or obligations under this Agreement may be sold, sublicensed, or other- 
wise assigned without prior written consent of Charles L. Joseph. 

13. YOU ACKNOWLEDGE THAT (a) YOU HAVE READ THIS ENTIRE AGREEMENT AND 
AGREE TO BE BOUND BY ITS TERMS AND CONDITIONS; (b) THIS AGREEMENT IS THE 
COMPLETE AND AND EXCLUSIVE STATEMENT OF THE UNDERSTANDING AND CON- 
TRACT BETWEEN US AND SUPERSEDES ANY AND ALL PRIOR ORAL OR WRITTEN COM- 
MUNICATIONS RELATING TO THE SUBJECT MATTER HEREOF; AND (c) THIS AGREEMENT 
MAY NOT BE MODIFIED, AMENDED OR IN ANY WAY ALTERED EXCEPT IN WRITING AND 
SIGNED BY BOTH YOU AND THE FIRST AUTHOR. 


ORIGINAL PAGE IS 
Oh POOR QUALITY 
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APPENDIX B 


Data Structures in NSLAP 


There axe 3 structures used by MSLAP, "mp". "mcntrl". and dtl . HI 
global in nature, meaning they can be accessed from any subroutine. The firs , 
•W’ is the primary structure and the one most frequently used to modify a 
particular behavior of MSLAP. "mcntrl" is a second structure used primarily for 
record keeping and some control. Greater care should be exercised in make 
adjustments to "mcntrl" since MSLAP could become "confused" and may lead to 
SSo- results that may not be obvious. This structure ("mcntrl") is intended 
for the advanced user of MSLAP. The final structure is "dtl", which contains 
the measurements that are written to the output data file. This "dtl" structure 
can be edited using option 4 in the beginning Main MSLAP Menu. Below is a list 
of th» parsjn.ters .long with a description of each. 


The primary (non-data etorag.) etmctur. used b, MS UP 1. "m p" . Hs lapParameter. 

mp”dget - a _ flag”used to indicate the presence of a dget routine. During the 
Betup stage, MSLAP calls each dget routine (i.e. dgetl, dget2, .... 
dget5). At this time mp.cntrl equals -1, a flag to indicate the 
initial call. If the user-installed dget routine is to be 
recognized by MSLAP, it must return mp.dget not equal to zero on 
this initial call. The dummy dget routines leave mp.dget=0. 

mp.cntrl - The control flag for MSLAP operations. Some care should be 
exercised when changing the value of mp.cntrl in any USERPROG, but 
generally speaking mp.cntrl can be safely set in any DGET. 

mp.cntrl ■ -1 «> This is the first time the dget routine is being 

called. Hardware graphics are configured. 

mp.cntrl « 1 «> No problems have been encountered in the sub- 

routine. OK flag for storing the measurement. 

mp.cntrl * -99 *> Abort Condition, MSLAP will exit. 

mp.cntrl * -10 «> Make HardCopy Flag. 

mp.cntrl * 0 *> Subroutine Returned as if NEVER CALLED. MSLAP 

will just replot spectra and start over . 

mp.cntrl - 2 «> Get NEXT Spectra. 

mp. order - is used for echelle spectra. It can be used to select a portion 
of the spectra when more than one record is placed in a given data 
file. (Standard MSLAP does not make use of this parameter; it is 
supplied solely for the user's convenience [e.g. it can be displayed 
using plotlabl or 2.]) 

mp.CAM - Specifies a camera number. Many satellites have more than one 
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and it is of tan usaful to racord this information. 

(Standard MSL1P doas not make usa of this parameter; it is 
supplied solely for the user's convenience [e.g. it can be displayed 
using plotlabl or 2.3) 


■p.SHO - a flag that causes the graphics to smooth the data by a 3 point 
running box car. (This feature is becoming obsolete and may not 
available in future versions of MSLAP.) It is best to perform 
any smooth ing inside the user’s DGET or USERPROG routine and then 
adjust the mp.cohfac value accordingly. 


■p.cohfac - the coherence factor of the noise in the data. (See MSLAP the 
Manuel for details.) Basically, this is a measure of the linear 
independence of the one pixel to other nearby pixels. For example, 
data that have been smoothed by a 5-point ru nnin g box car would 
have a coherence factor of 5. 


mp.bg - the background level, if not zero. 

mp.bgerr - the 1-sigma uncertainty in the background level. This value is 
VERY IMPORTANT since in many applications this error is the 
dom inan t source. (See the MSLAP Manuel for details.) 

mp. window - this parameter is used to define a fixed-sized integration window. 
In data with poor signal-to-noise ratios or in cases where the 
intrinsic strength of a spectral feature is expected to be weak 
compared to the noise, this feature of MSLAP allow the user to 
make unbiased measurements of the signal strength based on info, 
obtain from other spectral features. This parameter is set and/or 
used in response to the MENU option which selects the order of the 
polynomial to fit the continuum. (See the MSLAP Manuel for 
information regarding the "Predetermined Window" option.) 

mp.DTY - a parameter with integer values in the range [1, ...,5] which 
selects which [DGET1 , .... DGET 5] routine will be called. 

mp.STAR - Contains the Output File name. This string plus the " .DTL" 
ext ens ion form the primary output data file that MSLAP creates. 

sq>.FNAM - Contains the Input File name from which the spectral data will be 
retrieved. This string may include a directory path as well. 

mp.poly - Holds the order number of the polynomial that was used to fit the 
continuum. 

mp.SNR - is the APPARENT Signal-to-Noise Ratio. It is the RMS value of the 
residuals of the read data minus the polynomial fit. The Real SNR 
is obtained by dividing mp.SNR by the square-root of the the 
coherence factor [i.e. Real SNR * mp . SNR/ (mp . cohf ac 0.5)]. This 


RMS Talus is added in quadrature with the mp.bgerr to estimate the 
total uncertainties. 
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A structure called "mcntrl" is primarily used for record keeping while KSLAP is 
running. It is recommended that only advance users of MSLAP adjust these para- 
meters. 


mcntrl . I - is the pointer of the number of measurements that have been made 
If a "output. DTL" file exists at start up. this file is opened and 
mcntrl. I is set to the number of non-zero data entries. The user 
is then prompted to accept or reinitialize this pointer. 


mcntrl. date a string holding the date obtained from the system clock 


mcntrl . int opt an integer flag 
Normally, mcntrl . intopt= 
is called. 


indicating which integrating routine is to be 
1. If mcntrl. intopt=0, the ALT_INT routine 


used. 


mcntrl. WE the number of continuum points used in the polynomial fit and in 

the uncertainty calculations. 


mcntrl. ESAV the RMS uncertainty, calculated from the residual differences 
between the real continuum and the polynomial fit. 


mcntrl. eO Background's contribution to the zeroth moment (Equivalent Width). 
If multiple measurements are made simultaneously, only the first 
measurement is stored in this parameter. The same is true for 
mcntrl. ei and mcntrl.e2 

mcntrl. el same as mcntrl. eO, except it is for the first moment. 

mcntrl. e2 same as mcntrl. eO, except it is for the second moment. 

mcntrl. wtol wavelength tolerance used in testing for a match between the 
selected observed wavelength and comparison to the USER Loofc Up 
Table of rest wavelengths. Normally, this is taken to be 10 
spectral elements wide or 0.5 Angstroms, which ever is largest. 

mcntrl. libr a string holding the main "library" directory of routines and 
tabled data. 

mcntrl. mtot the largest polynomial that will be fit by MSLAP. Standard MSLAP 
uses mcntrl .mtot«5. A larger value will automatically fit corresp- 
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Hold* *11 of the First Momsnts of the profiles. 

dtl.fme Hold* all of the Error* associated with the First Moments. 

dtl.*m Hold* all of the Second Moments of the profiles. 

dtl.sme Hold* all of the Errors associated with the Second Moments. 

dtl.com Holds all of the byte arrays containing the information that was 

entered a* comments. This part is similar to dtl.wl above. 

dtl.up Holds all of the UserParameter Arrays. Each dtl. up holds a 30 

element floating-point vector that is reserved for the user’s 
exclusive use* 


MSLAP SITE LICENSE AGREEMENT 


TERMS AND CONDITIONS 

1. The title and full ownership rights to each and every part of die Modular Spectral Line Analysis Pro- 
gram (hereinafter MSLAP) shall remain the sole property of Charles L. Joseph and Edward B. Jenkins 
(hereinafter the authors). You shall acquire no rights in the Program other than as expressly granted in 
this Agreement. Each portion of MSLAP constitutes valuable proprietary assets of the authors, em- 
bodying substantia] creative efforts and significant expenditure of time. 

YOU MAY NOT USE, COPY, TRANSFER, OR MODIFY THE PROGRAM OR ANY COPY, 
MODIFICATION, OR MERGED PORTION THEREOF, IN WHOLE OR IN PART, EXCEPT AS EX- 
PRESSLY PROVIDED IN THIS AGREEMENT. YOUR LICENSE TO USE MSLAP IS AUTOMATI- 
CALLY TERMINATED IF YOU TRANSFER POSSESSION OF ANY COPY, MODIFICATION OR 
MERGED PORTION OF THE PROGRAM TO ANOTHER PARTY. 

Version 1.0 of MSLAP is licensed free of charge only to sites that are actively engaged in astronomical 
research for the sole purpose of pursuing astronomical research and that are in compliance with the 
terms of this Agreement Possessing MSLAP in whole or in part indicates that the user has accepted all 
of the terms of this agreement Violation of any terms of this Agreement shall make you liable for full 
payment for any damages incurred to the authors. 

The provisions of this section 1 shall survive any termination of this Agreement. 

2. The authors provide MSLAP and license its use on a single site bases. A site is defined for the pur- 
poses of this Agreement as a single Department Center, or Institute that is located on a single campus 
of a university or government agency. Permission is granted to you to use or to copy MSLAP to any 
machine on your site that is controlled by your Department, Center, or Institute, provided such action is 
in full compliance with all of the terms of this Agreement 

3. MSLAP IS PROVIDED ’AS IS’ WITHOUT WARRANTY OF ANY KIND, EITHER EXPRESSED 
OR IMPLIED. NO WARRANTY OF ANY KIND IS IMPLIED OF MERCHANTABILITY OR FIT- 
NESS FOR A PARTICULAR APPLICATION. THE ENTIRE RISK AS TO THE QUALITY AND 
PERFORMANCE OF THE PROGRAM IS WITH YOU AND NOT THE AUTHORS NOR THEIR 
AGENTS. 

In no event will the authors or their agents be liable to you for any damages arising out of use or ina- 
bility to use MSLAP, even if the authors or their agents have been advised of the possibility of such 
damages, or for any claim by any other party. If any third-party claims arise against you, you shall 
bare all such liabilities. 

4. Standard version 1.0 of MSLAP consists of the 34 routines licensed to you that initially are contained 
in the files: main .ms lap, mas ter .a ux, compare.pro, posto.pro, edatdtl.pro, man tail, pro, intgrtpro, 
plotconfig.pro, and cog.pro as well as all of the associated tables found in files with the extension 
".tab". You are permitted to alter this standard source code as you desire for your personal or site’s 
use, but you must maintain the original copyright notice in both the source file and in the initial start up 
of MSLAP, and you may NOT distribute the altered code off site without written authorization from 
Charles L. Joseph. If your site has 2 or more visitors in a given year who use MSLAP, you are con- 
sidered a "guest user facility", and you must comply with the conditions spelled out in section 8 of this 
Agreement Any alterations to the standard source code must also meet the documentation require- 
ments spelled out in section 6 of this Agreement 


5. A set of routines, contained in the file m slap .pro, are designed for individualized customization of 
MS LAP. You are permitted to create new modules (subroutines) of the USERPROG’s or to make 
adaptations of the routines originally found in the file mslap.pro, and to distribute these customized 
modules to other sites, provided the new code meets the requirement in section 6. 

6. Any and all modifications must comply with the following: 

L All modified code must describe all parameters in each procedure definition statement. 

ii. All modified code must indicate the original author, author’s affiliated institution, and a 
m odificati on history including each and every major participant involved in the modification. 

iii. If the altered code is designed to override a particular behavior of standard MSLAP, a concise 
description of the former and new behavior must be output in a conspicuous manner and also in- 
cluded in the source code documentation. 

iv. The ALT_INT routine has special rules. This routine is used to provide an alternate type of in- 
tegration. The following disclaimer must appear in a conspicuous manner whenever the routine is 
called. Standard MSLAP with a working instruction window at the lower left of the screen au- 
tomatically produces this disclaimer. The following lines from that disclaimer are listed below 
and may NOT be removed from the file intgrtpro. 

print,’ ’,string(7B) 

print^tring(7B),’No Alternate Integration Routine is Provided’ 
print,’in standard M.S1.AP. - version 1.0’ 
print,’ ’ 

print/Origin of the Alt. Integration Routine is as follows:’ 

7. A library of user-generated data-getting routines are made available to the astronomical community 
and users are encouraged to add to this library. MSLAP supports any combination of 5 or less of these 
routines at one time. (See the MSLAP Documentation Manuel for the general implementation of dgetl , 
dgel2, dget3, dget4, and dget5 procedures and their appropriate protocols.) All "dget" routines are con- 
sidered to be formally outside of MSLAP and not subject to die terms of this Agreement However, it 
is strongly recommended that the documentation rules listed above be implemented to maximize the 
overall utility of any user-contributed routines for data retrieval. 

8. If you represent a guest user facility, defined as any site where more than 2 visitors use MSLAP in a 
given year, you arc not permitted to incorporate modifications that remove the modularity of MSLAP, 
thus making it difficult for guest users to customize MSLAP. For example, you must avoid reserving 
more than 10 of the 30-element vector called "up" (UserParameter) since such action would severely 
inhibit the ability of a user to store and manipulate customized calculations. 

Some version of the original file mslap.pro must be retained and be made available to individual users 
for their personal customization of MSLAP. 

If you modify any portion of MSLAP including the original file mslap.pro to override a particular 
behavior of standard MSLAP, a concise description of the former and new behavior must be output in a 
conspicuous manner when MSLAP is running and also must be included in the source code documenta- 
tion. 

9. Permission is granted to incorporate MSLAP into a large driver program, provided the name MSLAP 
and the original copyright notice is maintained in any menu or other solicitation to activate MSLAP. 
Any incorporation into a driver program is also subject to sections 4, 5, 6, and 8. 


-3- 


10. Your site is considered to be a Distribution Center and you may transfer copies of the STANDARD, 
version 1.0 of MSLAP to other astronomical sites only if this section of this agreement is complete with 
all necessary signatures. 


Site Name and Address of the Department, Center, or Institute that is to be a a Distribution Center. 


Date and Signature of authorized representative of the above site. 


Name and Position of the person signing above. (Please print or type) 


Date and Signature of Charles L. Joseph, the first author. 

11. The license granted under this Agreement is effective until terminated. You may terminate this 
license by destroying all copies, modifications and merged portions of MSLAP in your possession and 
notifying Charles L. Joseph in writing of such destruction and termination. 

The license granted under this Agreement will terminate if you violate any of the terms and conditions 
of this Agreement The first author, Charles L. Joseph, reserves the right to terminate this Agreement 
if he believes that any of the terms and conditions of this Agreement have been violated. Notice of 
such termination shall be made in writing to you. You agree upon such termination to destroy promptly 
all copies, modifications, and merged portions of MSLAP in your possession and to certify to the first 
author that such action has been taken. 

12. None of your rights, duties or obligations under this Agreement may be sold, sublicensed, or other- 
wise assigned without prior written consent of Charles L. Joseph. 

13. YOU ACKNOWLEDGE THAT (a) YOU HAVE READ THIS ENTIRE AGREEMENT AND 
AGREE TO BE BOUND BY ITS TERMS AND CONDITIONS; (b) THIS AGREEMENT IS THE 
COMPLETE AND AND EXCLUSIVE STATEMENT OF THE UNDERSTANDING AND CON- 
TRACT BETWEEN US AND SUPERSEDES ANY AND ALL PRIOR ORAL OR WRITTEN COM- 
MUNICATIONS RELATING TO THE SUBJECT MATTER HEREOF; AND (c) THIS AGREEMENT 
MAY NOT BE MODIFIED, AMENDED OR IN ANY WAY ALTERED EXCEPT IN WRITING AND 
SIGNED BY BOTH YOU AND THE FIRST AUTHOR. 


